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1  INTRODUCTION 


1 . 1  GENERAL 

The  conventional  approach  for  implementing  adaptive 
signal  cancellers  to  reject  radar  interference  is  to  build  a 
separate  electronic  loop  for  each  required  degree  of  freedom. 
The  process  consists  of  cross-correlating  the  radar  signal 
with  interference  samples  received  by  wide-beamwidth  auxilia¬ 
ry  antennae.  These  correlation  values  are  then  used  to  ad¬ 
just  the  interference  waveforms  in  amplitude  and  phase  so 
that  a  precise  subtraction  of  the  offending  interference  can 
be  accomplished.  [1],  [2],  [3]. 

The  interference  components  can  arrive  reflected  from 
various  objects  causing  "multi-path"  delay.  Each  resolvable 
time  delay  requires  a  new  degree  of  freedom  in  the  proces¬ 
sing,  and  in  many  cases  an  even  denser  time  sampling  is  re¬ 
quired.  This  requirement  may  be  thought  of  as  the  need  of 
the  processor  to  have  a  sample  of  the  interference  waveform 
of  the  correct  structure  and  timing  to  cancel  such  a  waveform 
in  the  radar  signal. 

In  addition,  if  there  are  independent  interference 
sources  in  the  radar  space,  a  new  set  of  adaptive  loops  with 
multi-path  capability  must  be  provided  for  each  of  these 
sources.  The  result  is  often  an  expensively  large  number  of 
electronic  loops. 

An  optical  implementation  of  such  a  processor  is  the 
subject  of  the  effort  described  in  this  report.  The  optical 
AOAP  concept  was  first  suggested  in  1975  by  Dr.  Prank  R. 
Dickey,  Jr.  of  the  General  Electric  Co.,  Syracuse,  N.Y.  [4] 
In  the  period  of  1976  -  1980,  the  General  Electric  Co.  was 
diligent  in  developing  this  idea,  culminating  in  several 
internal  publications  [5],  [6],  and  the  effort  for  the  Rome 
Air  Development  Center  which  is  described  in  this  report. 

In  the  optical  configuration  which  has  been  developed, 
the  required  operations  of  multiplication  and  time  delay  are 
provided  by  acousto-optical  (AO)  delay  lines.  The  required 
time  integration  is  provided  by  an  image  converter  having  a 
suitable  delay  time.  In  the  optical  realization,  each  re¬ 
solvable  optical  element  along  the  AO  delay  lines  used  as 
input  modulators  represents  an  additional  delay  time  degree 
of  freedom.  Thus  we  find  a  "continuum  of  correlation 
weights"  developed  spatially  along  the  delay  line.  The  opti¬ 
cal  dimension  transverse  to  the  direction  of  acoustic  propa¬ 
gation  can  be  used  to  provide  another  array  dimension. 


Each  additional  interference  source  can  be  properly 
cancelled  by  providing  another  AO  column  in  the  optical  aper¬ 
ture.  Thus  with  an  array  of  delay  lines  in  the  aperture,  a 
single  optical  adaptive  system  can  be  equivalent  to  many 
hundreds  of  electronic  control  loops. 

In  the  following  report,  the  details  of  the  configura¬ 
tion  which  has  been  developed  are  described.  An  analysis  of 
the  system  is  presented  which  investigates  the  effects  of 
system  parameters  and  the  loop  gain  on  cancellation  perfor¬ 
mance.  The  experimental  feasibility  model  of  the  system 
which  has  been  built  is  described,  and  the  results  obtained 
with  this  system  presented.  Limitations  in  performance  which 
were  experienced  in  the  effort  are  analytically  evaluated. 
Recommendations  for  improved  performance  are  made.  In  addi¬ 
tion,  a  theoretical  study  is  included  which  evaluates  an 
adaptive  canceller  when  operating  against  a  single  interfer¬ 
ence  source,  with  a  single  multi-path  environment,-  i.e.  the 
possibility  of  direct  path  reception  of  the  interference,  and 
a  single  alternative  delayed  path.  The  results  of  this 
study,  which  are  contained  in  an  extensive  set  of  computer 
runs,  will  serve  as  a  guide  to  evaluating  the  performance  of 
the  eventual  closed  prototype  systems  which  will  be  tested 
first  with  laboratory  signals,  and  then  in  a  radar  environ¬ 
ment. 

1.2  ADAPTIVE  CANCELLER  LOOPS 

A  variety  of  wide-band  electronic  adaptive  systems  have 
been  built  to  perform  interference  cancellation  in  radar  and 
sonar  systems.  Interference  is  a  propagating  component  in 
the  radar  or  sonar  field,  which  can  be  sensed  by  an  auxiliary 
antenna  or  array.  This  interference  component  is  detected 
and  is  cross-correlated  with  the  output  from  the  primary,  or 
main  antenna  or  array.  Any  correlation  thus  found  indicates 
the  presence  of  an  unwanted  component  in  the  signal  being 
generated  by  the  system.  This  interference  can  be  generated 
from  either  natural  or  hostile  sources. 

When  a  correlation  is  established,  the  precise  complex 
value  of  this  correlation  is  produced  by  the  signal  proces¬ 
sing  system.  This  value  is  then  imposed  (multiplied)  onto 
the  auxiliary  (interference)  signal  to  derive  a  version  of 
the  interference  adjusted  in  amplitude  and  phase,  that  when 
subtracted  from  the  main  signal  will  cause  interference  can¬ 
cellation.  The  residue  of  this  subtraction  is  fed  back  into 
the  adaptive  loop  thus  formed,  in  the  manner  of  a  closed  loop 
control  system,  which  causes  adjustment  of  the  correlation 
coefficient  to  improve  the  cancellation.  With  proper  inte¬ 
gration  in  the  control  loop,  a  null  is  developed  in  regard  to 
the  interference,  and  the  desired  signal  is  passed  without 
appreciable  degradation. 


through  a  tapped  delay  line  to  account  for  all  the  expected 
multi-path  time  delays  which  might  be  experienced  in  the  main 
channel  for  the  given  interference  component.  Each  of  these 
taps  is  provided  with  a  separate  electronic  loop.  The  appro¬ 
priate  one  of  this  subset  of  loops  matches  the  actual  time 
delay  of  the  noise  path.  This  particular  loop  will  be  char¬ 
acterized  by  a  high  correlation  between  its  delayed  auxiliary 
signal  and  the  multi-path  delayed  interference  component 
found  in  the  main  signal  m(t).  This  correlation  is  developed 
by  a  multiplication  between  the  main  and  auxiliary  signals 
followed  by  a  time  integration.  This  develops  a  correlation 
"weight"  W,  which  then  multiplies  the  auxiliary  signal,  thus 
adjusting  it  in  amplitude  and  phase  so  that  a  cancellation 
occurs  with  subtraction. 

A  whole  set  of  loops  spaced  in  time  delay  is  made  avail¬ 
able  for  each  auxiliary  input  in  the  system,  which  is  a  de¬ 
sign  response  to  the  expectation  of  a  multiplicity  of  noise 
sources.  Thus  a  degree  of  freedom,  i.e.  a  separate  adaptive 
loop,  is  required  in  the  most  complete  design  for  each  combi¬ 
nation  of  auxiliary  input  and  time  delay  tap.  The  totality 
of  signal  outputs  of  the  loops,  each  seeking  its  own  null  for 
its  particular  combination  of  input  waveform  and  time  delay, 
is  summed  (summing  bus  E )  and  the  result  is  subtracted  from 
the  main  loop.  Analysis  shows  that  the  required  spacing  in 
time  delay  for  good  cancellation  is  the  reciprocal  bandwidth 
of  the  system,  assuming  an  unlimited  offset  time  window 
(tapped  delay  line  length)  and  quadrature  sampling.  For  a 
finite  window,  i.e.  a  finite  number  of  taps,  the  tap  density 
must  be  even  finer,  according  to  interpolation  theory.  Fur¬ 
ther,  enough  delay  and  corresponding  number  of  taps  must  be 
provided  to  accommodate  multipath  components  in  the  auxiliary 
signal,  and  the  combination  may  exceed  the  primary  signal 
delay  by  a  large  factor.  [2]  As  a  result,  in  many  practical 
cases,  many  hundreds  of  loops  are  required.  Obviously,  re¬ 
duction  of  the  tap  density  has  a  direct  effect  on  the  amount 
of  equipment  required.  The  resulting  compromise  of  perfor¬ 
mance  is  often  a  difficult  issue. 


2  MATHEMATICAL  DESCRIPTION  OF  OPTICAL  CANCELLER 


The  functions  of  the  two  multiplications,  the  time  inte¬ 
gration,  and  the  delay  line  can  be  provided  by  an  optical 
system  which  includes  acousto-optic  (AO)  Bragg  lines.  The 
time  integration  can  be  supplied  for  example,  by  an  image 
converter  for  which  a  phosphor  screen  provides  an  optical 
output  with  a  particular  decay  time.  Other  devices  are  pos¬ 
sible  and  will  be  discussed. 


2 . 1  SIGNAL  CORRELATION 


2.1.1  Time  Integrating  Correlator 


One  possible  arrangement  of  such  a  system  is  shown  in 
Figure  2-1.  The  light  source  is  modulated  with  the  cancelled 
main  signal.  An  optical  modulator  such  as  an  electro-optical 
(EO)  Pockel's  cell  can  be  used  with  a  cw  laser,  or  a  semicon¬ 
ductor  laser  diode  can  be  modulated  directly.  Thus  all  of 
the  light  entering  the  system  varies  in  amplitude  according 
to  the  main  signal  after  it  has  been  cancelled,  i.e.  r(t)  = 
m(t)  -  m(t)  in  Figure  2-1.  The  light  then  passes  through  a 
Bragg  cell  delay  line,  DL1 .  The  acoustic  wave  launched  in 
this  cell  is  obtained  from  the  auxiliary  channel.  This  func¬ 
tion  is  denoted  as  iT(x-vat),  where  va  is  the  acoustic 
velocity.  This  form  of  the  argument  expresses  the  fact  that 
it  is  a  traveling  wave.  At  each  point  a  multiplication  is 
achieved  between  a  and  ?.  The  time  offset  between  the  two 
functions  varies  with  position  along  the  delay  line.  This 
product  function  is  imaged  point  by  point  on  the  photocathode 
of  the  image  converter.  The  spatial  pattern  is  transferred 
to  the  output  phosphor  screen  of  the  image  converter,  which 
performs  the  desired  time  integration.  It  is  presumed  that 
the  integration  is  characterized  by  an  exponentially  decaying 
time  weighting  function.  For  this  application  a  Pi  phosphor 
was  selected,  providing  a  time  constant  on  the  order  of  10 
milliseconds.  The  spatial  pattern  thus__developed  represents 
the  correlation  function  between  a  and  r,  over  a  time  history 
corresponding  to  the  exponential  decay  of  the  optical  device. 
The  correlation  function  is  presented  in  the  spatial  coordi¬ 
nate  as  a  continuous  function  of  signal  offset.  This  corre¬ 
sponds  to  the  "correlation  weight"  of  the  conventional  time- 
sampled  architecture  composed  oj.  electronic  loops.  If  we 
denote  this  weight  function  as  W(x,t)  we  have: 
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where  t  =  current  time 

t  =  time  history  variable 

T0  =  1/e  decay  time  of  optical  integrator 

x  =  position  along  delay  line  or  delay  line 

image,  representing  correlation 
offset  time 

The  relation  between  correlation  offset  time  t'  and 
delay  line  position  x  is  given  in  terms  of  the  acoustic  ve¬ 
locity  va  as 


=  X/V; 


(2-2) 


Both  a(x,t)  and  r(t)  are  carrier  borne  signals.  The 
carrier  is  the  key  to  handling  signals  in  both  amplitude  and 
phase  which  are  thus  mathematically  representable  as  complex 
quantities.  The  i.f.  carrier  therefore  assumes  the  role  of 
the  original  radar  or  system  propagation  carrier.  Thus  if 
am  and  rm  are  the  complex  modulation  signals,  a  and  r  may 
be  represented  in  the  following  way: 
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r  ( t )  =  r  0  ®  +  r  +  r  e  ^ 
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where  rm  =  complex  modulation  signal  (over-bars 

indicate  complex  quantities) 

Fm*  =  conjugate  of  rJn 

wa  =  acoustic  carrier  radian  frequency. 
Similarly, 
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(2-4) 


This  form  includes  the  case  of  optical  phase  modulation 
in  the  delay  line.  Phase  and  amplitude  modulation  are  dif¬ 
ferentiated  by  the  relative  phase  of  the  am  terms  and  aQ. 

The  light  beam,  bearing  signal  r,  passes  through  the  delay 
line  which  contains  a,  and  thus  a  product  is  formed.  This 
product  is  imaged  on  the  image  converter  and  integrated 
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there  with  time  constant  T0,  according  to  equation  (2-1). 
The  product  a(x,t)7(t)  from  (2-3)  and  (2-4)  forms  9  terms. 
These  terms  all  have  a  factor 


2  . 


For  n  =  1,2  the  term  is  characterized  by  an  acoustic 
carrier  frequency,  and  hence  will  not  survive  the  integra¬ 
tion.  Only  the  three  n  =  0  terms  are  thus  of  interest,  and 
these  lead  to: 


W(x,t)  =  WQ  + 
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where  WQ  =  r0a0  =  optical  bias  term 


(2-5) 
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This_  is  recognized  as  the  complex  correlation  function 
between  am  and  rm,  involving  exponentially  weighted  inte¬ 
gration.  Equation  (2-5)  is  equivalent  to 


W(x,t)  =  WQ  +  2  | Wm|  cos  (u»a  —  +  4>w  )  (2-7) 

va 

where  <)>w  -  phase  angle  of  correlation  weight. 

This  may  produce  optical  amplitude  or  phase  modulation 
depending  on  the  relative  phase  of  WD.  Since  the  delay 
line  is  an  optical  phase  modulator,  and  W0  is  the  product 
of  rp  and  a^,  equation  (2-4)  and  therefore  equation  (2-5) 
implies  optical  phase  modulation.  This  condition  may  be 
converted  to  an  equal  combination  of  amplitude  and  phase 
modulation  by  eliminating  one  sideband,  e.g.  the  third  term 
in  equation  (2-4).  The  spectrum  of  the  product  function 
(corresponding  to  the  acoustic  carrier)  is  spatially  dis¬ 
played  in  the  focal  plane  of  lens  L2 .  The  extra  term  can 
thus  be  eliminated  in  the  optical  spatial  frequency  plane, 
between  lenses  L2  and  L3  with  a  simple  optical  stop.  This 
third  term  must  be  optically  removed  only  when  the  array  of 
optical  orders  are  being  generated  in  a  symmetrical  manner, 
such  as  by  a  Raman-Nath  delay  line.  This  is  a  line  which  is 
thin  enough  in  the  direction  of  optical  travel,  so  that  it 
effectively  acts  as  a  phase  modulator  in  a  single  plane.  In 
the  case  of  a  Bragg  line,  or  a  "thick"  line,  when  adjusted  at 
the  Bragg  angle,  only  one  strong  diffracted  optical  order  is 
generated.  This  order  is  reconstructively  related  to  the 
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input  beam  over  the  volume  of  the  delay  line.  Thus,  essen¬ 
tially  only  two  orders  emanate  from  the  delay  line,  and  the 
desired  amplitude  modulation  is  automatic,  without  the  need 
of  a  stop. 

2.1.2  Space  Integrating  Correlator 

The  T7(x,t)  pattern  developed  by  the  integrating  device 
(image  converter)  is  now  imaged  to  a  second  optical  delay 
line,  DL2,  similar  to  the  first  delay  line  (see  Figure  2-1). 
Again  a  product  is  formed.  In  the  case  of  an  optically  co¬ 
herent  system,  the  required  amplitude  modulation  can  again  be 
provided  by  the  elimination  of  an  optical  order.  However,  if 
the  integrating  device  provides  a  light  pattern  which  is 
optically  incoherent,  some  different  scheme  must  be  used  to 
provide  the  required  amplitude  modulation.  An  approach  to 
the  incoherent  case  is  to  use  polarizers,  and  rely  on  the 
acoustically  induced  birefringence  in  the  delay  line  to  alter 
the  polarization.  It  can  be  shown  that  the  response  of  the 
output  light  intensity  to  a  birefringent  phase  shift  <t>  (phase 
difference  between  principal  axes)  is  given  by 


—  =  sin  2  (-) 
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(2-8) 


In  this  case,  a  birefringent  bias,  such  as  is  provided 
by  a  phase  plate,  must  be  used  to  linearize  the  output  ex¬ 
pressed  by  (2-8).  The  product  function  between  T7(x,t)  and 
5(x,t)  is  spatially  integrated  across  DL2  by  lens  L5,  which 
is  focussed  on  the  output  detector.  The  function  launched 
into  DL2  is  identical  to  the  signal  in  DL1 .  This  can  be 
understood  by  referring  to  Figure  1-1,  in  which  the  generic 
architecture  of  an  adaptive  loop  is  shown,  with  the  two  suc¬ 
cessive  multiplications  of  the  auxiliary  signal.  Thus  we 
have  for  the  detector  signal: 
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where 


xA  is  the  interval  of  spatial  integration  across 
DL2. 


This  again  produces  9  terms.  Now  all  but  three  of  the 
terms  have  a  spatial  frequency  corresponding  to  the  acoustic 
delay  line  carrier  frequency.  The  integration  expressed  by 
(2-9)  will  reject  these  spatial  carrier-borne  terms,  since 
this  is  a  spatial  integration.  The  result  is: 


where 
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(2-11) 

The  result  makes  clear  that  a  coherent  phase-sensitive 
operation  is  accomplished.  This  is  to  say  that  the  correct 
form  of  the  complex  correlation  function  between  a  and  r  is 
indicated  by  the  integral,  and  the  complex  "weight"  is  used 
to  adjust  the  amplitude  and  phase  of  a^  by  multiplication 
for  subsequent  subtraction  and  cancellation  of  the  input 
signal . 


2.1.3  Space-time  Transfer  of  Phase  Coherent  Signals 

It  is  noteworthy  to  observe  the  manner  in  which  the 
complex  value  of  the  signal  is  transferred  through  the  sys¬ 
tem,  with  phase  preserved.  This  happens  because  the  signals 
are  borne  on  either  a  time  or  space  carrier  throughout,  in 
the  manner  of  a  narrow  band  signal.  The  first  signal  applied 
to  the  system  is  the  modulation  on  the  light  source.  The 
phase  of  this  signal  is  with  respect  to  an  assumed  reference 
carrier  frequency,  which  is  also  the  acoustic  carrier  in  the 
delay  lines.  The  correlation  operation  between  the  signal 
applied  to  the  light  source,  r,  and  the  signal  in  DL1 ,  5, 
senses  the  phase  between  r  and  a.  This  phase  information  is 
imparted  to  the  spatial  signal  developed  in  the  time  integra¬ 
tor  (image  converter).  Here  the  physical  position  of  the 
spatial  carrier  cycles  with  respect  to  an  assumed  position 
denote  the  phase  of  the  correlation  pattern  (W).  As  the 
phase  difference  between  r  and  a  changes,  the  position  of  the 
carrier  cycles  in  W  will  move  accordingly.  The  dual  of  this 
sequence  occurs  in  the  second  optical  section.  Here  the 
spatial  W-pattern  interferes  optically  with  the  traveling 
acoustic  wave  in  DL2 ,  which  upon  spatial  integration  by  the 
final  lens,  causes  a  time  carrier  to  appear  in  the  detector. 
The  time  phase  of  this  output  depends  on  the  spatial  position 
of  the  cycles  in  the  W- function,  and  the  phase  of  the  acous¬ 
tic  wave  in  DL2.  Thus  a  transfer  of  signal  is  made  from  the 
time  domain  to  the  space  domain  and  back  again,  always  with 
the  phase  being  coherently  transferred  at  each  point. 

Inspection  of  (2-11)  reveals  that  the  proper  adaptive 
operation  is  realized.  Signal  r  is  inspected  for  a  partial 
containment  of  a  by  a  complex  cross-correlation  of  r  with  a 
at  an  arbitrary  offset  time  of  x/va.  This  correlation 
weight  is  applied  to  a,  again  offset  by  x/va.  The  result 
is  the  correct  subtraction  for  an  interference  component 
delayed  by  x/va.  Finally  this  result  is  integrated  over 


all  values  of  x  that  are  in  an  expected  range  to  develop  the 
complete  subtraction  signal,  (equation  2-10). 

2.2  CLOSED  LOOP  CONTROL  SYSTEM 

When  the  subtraction  is  made,  the  system  becomes  a 
closed  control  system.  Some  insight  into  the  performance  of 
this  closed  feedback  configuration  can  be  gained  by  assuming 
that  a  unity  amplitude,  zero  phase  signal  is  introduced  into 
the  system  as  the  auxiliary  input,  and  a  constant  signal,  as 
represented  by  a  step  function  is  applied  as  the  main  signal 
input.  Hopefully,  the  system  will  adapt  to  this  constant 
signal,  and,  since  the  main  signal  now  resembles  the  auxilia 
ry  signal  within  a  proportional  constant,  the  entire  signal 
will  be  adaptively  cancelled  in  the  limit,  and  the  output 
will  approach  zero.  The  optical  aperture  is  assumed  to  be 
reasonably  small  in  this  simple  case,  so  that  only  one  effec 
tive  degree  of  freedom  need  be  assumed  for  this  particular 
signal.  The  simplified  diagram  for  this  assumed  case  is 
shown  in  Figure  2-2.  With  the  unity  amplitude,  zero  phase 
auxiliary  signal,  the  two  multipliers  have  the  effect  of 
unity  multiplication.  The  well-known  closed  loop  response 
for  a  unity  forward-gain  control  loop  is 
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where  H(u>)  is  the  open  loop  frequency  response. 
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Figure  2-2.  Simplified  Control  Loop 

As  shown  in  Figure  2-2,  the  transfer  function  e^9 
is  placed  in  the  loop  to  account  for  any  incidental  phase 
shift  in  the  loop.  A  single-break  integrator  is  assumed  to 
model  the  optical  storage  device,  such  as  an  image  converter 
Collecting  the  appropriate  factors  and  substituting  these 
into  (2-12)  produces  the  closed  loop  transfer  function: 


1  +  H 


(2-13) 


1  +  jwT~ 


This  function  has  a  pole  at  the  complex  frequency 
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As  the  spurious  phase  shift  in  the  loop,  0 ,  varies,  the 
location  of  this  resonant  pole  will  move  along  the  circular 
trajectory  shown  in  Figure  2-3.  Values  in  the  right  half 
plane  lead  to  loop  instability.  T.t  is  significant  that  the 
response  of  this  system  is  characterized  by  a  single  complex 
pole,  rather  than  the  usual  case  of  a  conjugate  pair.  It 
will  be  remembered  that  the  variable  which  is  mathematically 
represented  throughout  the  loop  is  a  carrier-borne  signal, 
and  the  phase  angle  of  the  variable  represents  the  instanta¬ 
neous  phase  of  the  signal  with  respect  to  an  assumed  stable 
carrier  frequency.  In  the  usual  control  analysis,  a  complex 
pair  of  poles  represents  a  sinusoidal  oscillation  in  the 
amplitude  of  the  controlled  variable  with  an  exponential 
growth  or  decay.  Here  a  single  pole  represents  a  decay  or 
growth  of  the  signal,  but  the  oscillation  now  manifests  it¬ 
self  as  a  forward  or  backward  progression  of  phase,  i.e.  a 
frequency  shift  for  a  time  variable,  and  a  spatial  motion  for 
a  space  variable.  Consider,  for  example,  the  spatial  sine- 
wave  pattern  representing  the  continuum  of  correlation 
weights,  on  the  image  intensifier.  If  the  input  is  suddenly 
reduced  to  zero,  implying  that  the  correlation  weights  will 
also  then  approach  zero,  a  moving  sine-wave  pattern  will  be 
observed,  with  the  a.c.  amplitude  performing  an  exponential 
decay.  If  the  modulation  function  for  this  behavior  is  plot¬ 
ted  in  the  complex  time  domain  plane,  the  trajectory  is  an 
inward  spiral,  approaching  zero  as  a  center  point.  This 
behavior  has  been  described  by  Morgan  (7] .  The  speed  of  this 
motion  will  depend  on  the  amount  of  spurious  phase  shift. 


Figure  2-3.  Trajectory  of  Closed  Loop  Pole 
in  Complex  Plane 


3  DESCRIPTION  OF  EXPERIMENTAL  SYSTEM 


The  experimental  system  was  arranged  essentially  as 
shown  in  Figure  2-1.  It  is  composed  of  two  parts.  The  first 
section,  from  the  light  source  to  the  optical  integrator 
{image  intensifier)  is  a  time-integrating  correlator.  The 
second  section,  from  the  integrator  to  the  output  detector  is 
a  space-integrating  correlator. 

3.1  TIME  INTEGRATING  CORRELATOR 

3.1.1  General 


The  operation  of  the  time  integrating  correlator  depends 
on  signal  modulation  of  the  light  source,  and  the  first  delay 
line.  The  main  residue  signal  is  imposed  on  the  light 
source,  and  the  auxiliary  signal  is  passed  through  the  first 
delay  line  in  the  form  of  a  sound  wave.  Thus  at  each  point 
along  the  delay  line  an  optical  multiplication  is  achieved, 
and  at  each  such  point  a  different  correlation  offset  time  is 
represented.  This  product  function  is  then  imaged  on  the 
image  intensifier  or  any  alternate  integrating  device  that 
might  be  used  to  develop  a  spatial  display  of  the  continuous 
correlation  function  between  the  two  signals. 

3.1.2  Light  Source 

One  possibility  for  a  satisfactory  light  source,  which 
would  have  several  desirable  characteristics,  is  a  visible 
gas  laser,  such  as  an  argon  ion  laser.  This  requires  the  use 
of  an  external  optical  time  modulator,  such  as  a  Pockel's 
cell.  One  advantage  of  this  approach  is  the  abundance  of 
input  coherent  light  that  can  be  thus  obtained.  On  the  order 
of  1  watt  of  usable  spatially  filtered  light  can  be  made 
available  in  the  illumination  of  the  first  delay  line.  To 
evaluate  this  approach  a  Lasermetrics  3079  Pockel's  Cell 
modulator  was  obtained.  This  is  a  KD*P  (KD2P04  or  potassium 
deuterium  phosphate)  crystal,  with  a  clear  aperture  of  2  mm 
diameter.  The  half-wave  voltage,  is  that  voltage 

required  to  produce  a  wave  retardation  difference  between  the 
two  optical  axes  of  it  radians.  For  the  3079  cell  VH/W 
is  162  v.  Experiments  were  conducted  at  60  MHz  using  a  tuned 
r.f.  driver.  This  experience  shows  that  it  is  very  difficult 
to  shield  the  high  r.f.  power  from  the  subsequent  electronics 
of  the  system.  In  particular  r.f.  feedthrough  can  be  experi¬ 
enced  which  feeds  directly  into  the  optical  detector  circuit. 
For  this  reason,  and  the  general  inconvenience  of  providing 
high  power  r.f.  at  a  substantial  bandwidth,  the  alternative 
of  using  a  semiconductor  laser  diode  was  elected. 


The  use  of  a  semi-conductor  laser  diode  offers  a  much 
more  convenient  signal  modulation  mechanism,  i.e.,  the  diode 
can  be  electrically  modulated  directly.  Modulation  band- 
widths  up  to  1  GHz  are  available,  as  far  as  the  intrinsic 
response,  of  the  device  is  concerned. 

A  Hitachi  1600  Ga  A1  As  diode  is  used  in  our  experiment. 
This  device  operates  nominally  at  830  nm,  i.e.  the  near 
infra-red  region.  The  diode  is  mounted  on  a  heat  sink  to 
help  avoid  mode-hopping  due  to  temperature  shift.  The  tem¬ 
perature  of  this  heat-sink  is  regulated  by  a  thermal  feedback 
control  loop  to  keep  the  temperature  within  approximately  a 
1 *C  range. 

Figure  3-1  shows  the  vulnerability  of  this  device  to 
mode-hopping  with  temperature  change.  Approximately  4  fre¬ 
quency  jumps  occur  over  a  temperature  range  of  6*C,  covering 
a  wavelength  shift  of  about  1  nm  according  to  Hitachi 
data.  [8] 


Temperature,  Th_(°C) 


Figure  3-1.  Peak  Wavelength  vs.  Heat  Sink  Temperature  of  an 
HLP1000  Diode  at  a  Constant  DC  Current* 


*  Taken  from  Reference  [8] 


Lasing  threshold  current  for  this  device  is  rated  to  be 
typically  70  ma.  The  device  was  operated  in  the  experiment 
at  a  bias  current  of  80  to  90  ma.  It  is  believed  that  ap¬ 
proximately  2  mw  of  optical  power  was  thus  obtained,  with  a 
maximum  rating  of  10  mw  available.  Added  to  this  bias  cur¬ 
rent  is  the  r.f.  modulation  developed  by  a  current  gain  am¬ 
plifier,  fed  by  a  signal  generator  with  an  output  of  approxi¬ 
mately  0  dBm.  By  crossing  the  laser  threshold  on  the  swings 
of  the  carrier  cycle,  100%  light  modulation  was  observed. 

3.1.3  Delay  Line 

The  delay  line  used  for  the  initial  experiment  is  an 
Isomet  Type  1201  AO  modulator.  This  is  a  glass  unit  operat¬ 
ing  at  a  center  frequency  of  37  MHz.  It  has  an  aperture  of  2 
mm  x  5  mm.  Our  tests  have  been  with  cw  excitation  only.  The 
unit  is  approximately  35  mm  thick,  which  puts  it  fairly  far 
into  the  Bragg  mode.  With  an  input  drive  of  1  watt  a  dif¬ 
fraction  efficiency  of  50%  is  easily  achieved.  The  same 
signal  generator  output  used  to  drive  the  laser  diode  light 
source  was  used  as  the  signal  for  this  delay  line.  The  sig¬ 
nal  generator  output  is  amplified  in  an  ENI  3-watt  r.f.  am¬ 
plifier.  This  is  enough  to  provide  50%  of  diffraction  effi¬ 
ciency,  which  means  that  50%  of  the  input  light  energy  in  the 
optical  zero-order  is  deflected  into  the  first  order  dif¬ 
fracted  order.  These  two  orders  are  then  used  to  provide  an 
amplitude  modulated  space  carrier  on  the  image  intensifier. 
Since  both  orders  are  equal,  100%  amplitude  modulation  is 
obtained. 

3.1.4  Optics 

The  imaging  optics  used  for  the  AOAP  experiment  consist 
of  telecentric  systems.  A  telecentric  system  uses  two  lenses 
with  a  common  front/back  focal  plane,  as  shown  in  Figure  3-3. 
The  imaging  then  occurs  between  the  front  focal  plane  of  the 
first  lens  (object  plane),  to  the  back  fecal  p-lene  of  the 
second  lens.  The  advantages  of  this  configuration  are  (1) 
exceptional  freedom  from  aberrations,  since  each  lens  is  used 
at  infinite  focus,  for  which  it  is  usually  designed,  (2) 
three-dimensional  undistorted  imaging  when  the  lenses  have 
identical  focal  lengths,  and  (3)  the  provision  between  the 
lenses  of  a  true  Fourier  spatial  frequency  plane,  when  coher¬ 
ent  light  is  used.  In  the  present  configuration,  two  Bausch 
and  Lomb  152  mm  Baltar  f/2.8  lenses  were  used,  as  well  as  a 
third  such  lens  to  accomplish  collimation  of  the  laser  beam. 
This  provides  precise  1:1  imaging  between  the  delay  line  and 
the  optical  integrator/image  converter. 
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Figure  3-3.  Telecentric  Imaging  Optics 


A  frequency  plane  is  thus  provided,  again  as  shown  in 
Figure  3-3,  in  which  the  spatial  dimension  xF,  is  related 
to  the  spatial  frequency  in  the  delay  line,  or  the  output 
image  at  the  image  intensifier  by  the  relation: 


XF  =  Fo 


(3-1  ) 


1  -  (X  fx)‘ 


where 


F0  =  focal  length 
fx  =  spatial  frequency 

For  paraxial  systems  such  as  this,  (meaning  all  light 
rays  are  close  to  being  parallel  to  the  optical  axis),  the 
relation  reduces  to 


xp  =  Fq  X  fx 


(3-2) 


For  the  acoustic  carrier  frequency  of  fQ  =  60  MHz,  an> 
an  acoustic  velocity  in  quartz  of  5.95  mm/u sec  we  have 

fXol  =  spatial  frequency  of  sound  carrier  in  delay 
line 


Lo  60  MHz 


v_  5.95  mm/w sec 

a 

10.7  cycles/mm 

=  acoustic  carrier  frequency 


At  37.5  MHz,  and  va  =  3.7  mm/ysec  for  the  Isomet  glass 
line  which  was  used,  we  have  the  same  spatial  frequency. 

From  (3-2),  the  spatial  separation  between  the  zero 
order,  and  the  diffracted  first  order,  corresponding  to  the 
acoustic  carrier  frequency  is: 

^o  X 

XF1  =  Fo  X  f xol  =  Fo  X  ~  =  Fo  ”T~  (3-4) 

va 

where 

F0  =  Lens  focal  length 

Fxol  =  Spatial  frequency  of  acoustic  wave  in  delay 
line 

fQ  =  Acoustic  carrier  frequency 

A  =  Acoustic  wavelength 

For  the  case  of  10.7  cycles/mm,  and  X  =  830  nm,  one 
obtains : 

xpi  =  152  x  .830  x  10-3  x  10.7 

=  1 . 35  mm 

On  this  scale,  it  is  quite  convenient  to  sort  orders  if 
necessary. 

3.1.5  Image  Converter 

Our  original  concept  has  been  to  use  an  image  intensi- 
f ier/converter  as  the  optical  temporal  integrator  required 
for  AOAP  operation.  Such  an  image  converter  will  accept  an 
input  light  pattern  which  is  composed  of  either  coherent  or 
incoherent  light  on  its  photocathode.  Electrons  generated 
from  the  photo-cathode  then  are  accelerated  to  a  phosphor 
output  screen,  where  the  output  image  is  found.  In  this  AOAP 
configuration  it  is  being  used  as  a  coherent  to  incoherent 
converter.  The  optical  wavelength,  at  the  same  time,  is 
converted  from  830  nm  to  the  mid-visible  band,  centered  in 
the  region  of  550  nm. 

The  desirable  range  of  integration  time  required  for  an 
AOAP  system  is  estimated  to  be  from  3  to  100  milliseconds. 

The  integration  time  is  defined  as  the  time  required  to  reach 
the  1/e  point  in  the  natural  exponential  decay  characteristic 
of  the  device.  Gain,  while  not  vital,  is  desirable.  The 


resolution  of  the  device  should  be  adequate  to  handle  the  10 
cycles/nun  spatial  frequency  generated  by  the  delay  lines. 

Thus  the  MTF  curve  should  be  well  above  50%  at  10  cycles/mm. 

Several  image  intensifier  configurations  could  be  se¬ 
lected  for  the  AOAP  system.  The  following  discussion  covers 
three  major  classes  of  this  type  of  device. 

Magnetically  focussed.  Figure  3-4  shows  the  three 
image  intensifier  configurations  presently  in  use.  The  first 
is  the  magnetically  focussed  intensifier  shown  in  Figure  3- 
4a.  This  configuration  consists  of  a  cylindrical  envelope, 
usually  made  of  glass,  with  windows  at  each  end.  The  inter¬ 
nal  surface  of  the  input  window  is  coated  with  a  thin  layer 
of  material  that  emits  free  electrons,  when  illuminated  with 
light  energy  in  the  acceptance  band.  Most  photoemitters  have 
a  response  that  includes  the  visible  spectrum  with  several 
extending  a  short  distance  to  adjacent  bands.  These  photo¬ 
emitters  are  the  same  surfaces  used  for  photocathodes  in 
photomultiplier  tubes.  The  opposite  end  of  this  intensifer 
envelope  contains  a  window  with  a  phosphor  layer.  Each  win¬ 
dow  has  a  transmission  reponse  matched  to  the  wavelength 
range  of  its  respective  surface.  These  devices  are  sometimes 
called  image  converters  because  the  input  radiation  can  be  in 
a  different  wavelength  from  the  output  radiation. 

Operation  of  this  device  proceeds  with  the  input  image 
interacting  with  the  photoemitter,  or  photocathode  to  produce 
an  electronic  image  of  low-energy  emitted  electrons.  A  large 
DC  potential  is  maintained  between  the  photocathode  and  the 
phosphor  coated  anode  at  the  other  end  of  the  envelope.  Thus 
the  low  energy  electronics  are  accelerated  to  provide  a  high 
energy  input  to  the  phosphor.  During  this  transit  from  the 
photocathode  to  the  phosphor  plate  the  electrons  are  focussed 
under  the  action  of  the  magnetic  field  from  the  focus  coil. 

In  the  arrangement  as  shown,  the  magnetic  focus  coil  will 
operate  at  a  magnification  of  nearly  unity  and  provides  an 
exceptionally  high  resolution  image.  It  is  possible  to  ob¬ 
tain  resolutions  as  high  as  100  cycles/mm  with  this  magnetic 
focus  arrangement. 

Proximity  focussed.  Figure  3-4b  shows  the  proximity 
focussed  diode  image  intensifier.  This  is  a  very  simple,  but 
effective  device.  Here  the  focus  is  effected  by  the  same 
electrostatic  field  as  is  used  for  acceleration.  This  is  of 
course,  a  compromise.  The  resolution  cannot  be  as  high  as 
the  magnetically  focussed  intensifier  because  here  there 
really  is  no  focus  action,  and  the  high  axial  field  merely 
holds  off  the  normal  defocussing  action  due  to  space  charge 
spreading  of  the  electronic  image.  Likewise  since  the  plates 
must  be  kept  very  close  one  cannot  use  as  high  an  acceler- 


ating  voltage,  and  therefore  the  gain  is  lower.  In  spite  of 
these  differences  the  proximity  focussed  diode  intensifier 
has  substantial  advantages.  It  is  physically  small,  light¬ 
weight,  has  few  parts,  has  only  two  electrodes  and  one  power 
supply,  gives  good  focussed  performance  at  any  voltage  in  its 
range.  In  short,  it  is  simple  and  reliable  and  is  the  lowest 
cost  option. 
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Figure  3-4.  Image  Intensifier  Configurations 
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MicroChannel  Plate.  The  third  dev-ice  is-  called  a 
microchannel  plate  because  of  its  construction.  This  inten¬ 
sifier  uses  a  large  collection  of  parallel  microchannel  chan¬ 
nel  to  provide  an  image  gain.  Each  microchannel  is  a  tiny 
secondary  emission  amplifier  made  by  coating  the  thin  semi¬ 
conductor  film  on  the  inside  of  a  micro  glass  cylinder. 

Under  the  axial  potential  generated  by  contacts  at  each  end 
of  this  tube,  one  entering  electron  impacts  the  film  on  the 
wall,  generating  more  than  one  secondary  electron.  This 
process  repeats  many  times  over  as  the  number  of  avalanching 
electrons  grow.  This  type  of  secondary  emission  amplifier  is 
a  relatively  low-noise  amplifier  capable  of  gains  as  high  as 
10“  or  105. 
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As  shown  by  the  sketch  of  Figure  3-4c,  the  microchannel 
plate  image  intensif'ier  is  arranged  to  be  operated  with  three 
power  supplies  (different  voltages).  Voltage  VI  places  a 
potential  between  the  photocathode  and  the  microchannel  input 
face.  This  provides  proximity  focus,  as  accomplished  in 
Figure  3-4b.  The  same  approach  is  used  to  focus  the  output 
electronic  image  to  the  phosphor  plate.  At  the  same  time 
these  potentials  increase  the  energy  of  the  initially  low- 
energy  electronic  image. 

Table  I  shows  a  device  performance  comparison  for  the 
three  configurations.  This  table  is  an  approximate  summary 
of  presently  available  devices.  One  can  quickly  get  a  gener¬ 
al  idea  of  the  capability  of  the  different  configurations. 

It  is  to  be  noted  that  in  cost  the  magnetic  focussed  is  the 
most  expensive  with  the  MCP  about  3/4  and  the  proximity  fo¬ 
cussed  about  1/4  this  cost. 


TABLE  I.  PERFORMANCE  COMPARISON  FOR  IMAGE  INTENSI FIER 

CONFIGURATIONS 


40  MM  APERTURE  (INPUT/OUTPUT) 
S-20  Photocathode,  P-20  Phosphor 


Device  Type 


1  )  Magnetically 
Focussed 

2)  Proximity 
Focussed 

3)  MicroChannel 
Plate  (MCP) 


Limiting 
2-5%  MTF 

Resolution 

(cycles/mm) 

Gain 

Cost 

(Approximate ) 

~ 50%  MTF 

85-90 

45 

60-120 

20$K-$25K 

45 

18 

20-50 

5$  K 

20-25 

10 

10- 

1  5$K 

When  comparing  gain,  the  proximity  focussed  is  the  low¬ 
est  gain,  with  the  magnetically  focussed  about  twice  this 
gain.  The  MCP  has  the  capability  of  exceedingly  high  gain 
and  in  fact,  this  is  the  major  advantage  of  the  microchannel 
plate.  The  gain  of  the  MCP  device  is  so  high  that  one  must 
exercise  care  not  to  exceed  some  rather  low  level  like  10-4 
foot-candles  or  the  output  is  saturated. 

Spatial  resolution  in  the  MCP  is  limited  by  the  spatial 
structure  of  the  microchannel  array  to  about  20  to  25  cy¬ 
cles/mm.  In  the  proximity  device,  the  resolution  is  limited 
to  nearly  45  cycles/mm  by  the  charge  spreading  in  the  proxim¬ 
ity  focus  region  between  the  photocathode  and  the  phosphor 
plate.  Resolution  in  the  magnetically  focussed  device  is 
limited  by  the  electron  optical  aberrations  of  the  magnetic 
focus  coil,  and  this  limit  is  about  100  cycles/mm. 

After  careful  consideration  and  comparison  of  the  vari¬ 
ous  capabilities  of  the  three  intensifier  configurations  we 
selected  the  proximity  focussed  device.  It  offers  the  best 
tradeoff  in  performance  specifications  for  the  intended  AOAP 
application. 

Most  of  the  present-day  image  intensifiers  are  intended 
to  intensify  an  image  for  direct  visual  observation  or  as  an 
input  to  a  TV  imaging  device  such  as  a  vidicon.  Consequently 
these  devices  are  usually  sold  with  the  S-20/P-20  combina¬ 
tion.  For  the  AOAP  system  we  want  a  special  phosphor  to 
obtain  the  proper  integration  time.  Table  II  shows  a  list  of 
phosphor  types  and  certain  charactistic  data.  Of  specific 
interest  are  those  phosphors  with  a  medium  to  medium  long 
persistence.  We  wanted  to  select  a  suitable  phosphor  with  a 
decay  to  10%  brightness  in  the  range  of  5  to  no  more  than  200 
milliseconds  to  yield  a  decay  to  50%  in  the  range  of  1  to 
about  50  milliseconds.  As  one  can  see,  there  are  at  least  6 
choices  open.  The  specific  selection  was  made  on  more  de¬ 
tailed  performance  features. 

The  resolution  of  10  and  20  cycles/mm  is  the  assumed 
resolution  at  a  50%  MTF  for  the  MCP  and  the  proximity  fo¬ 
cussed  devices  respectively.  We  see  that  the  MCP  does  not 
give  sufficient  resolution  to  readily  yield  300  or  more  cy¬ 
cles  per  side  dimension  except  in  the  largest  apertures.  The 
proximity  focussed  intensifier  of  40  mm  is  our  choice  since 
it  can  yield  the  desired  resolution  in  a  40  mm  diameter  aper¬ 
ture.  Larger  apertures  produce  inconveniences  in  the  optical 
system,  as  well  as  requiring  more  expensive  optics. 

Figure  3-5  shows  the  MTF  curve  for  the  ITT  40  mm  proxim¬ 
ity  focussed  (diode)  intensifier.  This  curve  shows  the  reso¬ 
lution  at  50%  MTF  is  about  17-18  cycles/mm. 


TABLE  II.  SUMMARY  OF  CHARACTERISTICS  OF  COMMONLY  USED 
"STANDARD"  PHOSPHOR  SCREENS  AS  PUBLISHED  BY  El A 
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Figure  3-5.  Modulation  Transfer  Function 


Two  ITT  proximity  focussed  intensifiers  were  obtained 
for  use  in  the  present  experiment. 


The  first  is  an  ITT  F4108  diode.  This  has  a  usable 
aperture  diameter  of  25  mm.  This  also  has  an  "extended  red" 
S-20  photoemitter.  The  spectral  response  of  this  emitter  is 
discussed  in  Section  4.3.1.  The  design  accelerating  voltage 
is  6  kilovolts.  At  this  level  the  tube  has  a  luminous  gain 
of  100. 


The  phosphor  screen  of  this  device  is  type  P-20.  The 
center  color  of  the  emission  is  yellow-green,  peaking  at 
550  nm.  The  response  time  is  in  the  region  of  100  usee. 

Thus  this  device  does  not  provide  the  integration  time  needed 
for  the  AOAP  system,  but  it  has  proved  very  useful  in  testing 
the  time-integrating  correlator.  It  was  used  to  achieve  the 
cw  correlation  pattern  described  in  the  next  section. 


The  second  device  obtained  is  an  ITT  F4109  unit.  This 
device,  a  non-standard  specification,  differs  from  the  F4108 
unit  in  three  ways: 


1)  It  has  a  usable  aperture  of  40  mm,  which  provides  a 
space-bandwidth  product  of  70  for  the  expanded 
AOAP,  corresponding  to  7  usee  at  10  MHz  bandwidth. 
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2)  It  has  a  P-1  phosphor,  which  is  characterized  by  a 
10  millisecond  response,  quite  adequate  for  AOAP 
integration. 

3)  It  does  not  have  a  fiber-optic  transfer  face-plate, 
used  on  standard  devices  for  ease  of  image  trans¬ 
fer.  This  eliminates  a  slight  reduction  of  MTF, 
and  slight  spatial  phase  distortion  encountered  in 
fiber  optic  faceplates. 

3.1.6  Time-integrating  Correlator  Experimental  Results 

Using  the  configuration  described  in  the  preceding  dis¬ 
cussion  the  auto-correlation  of  a  cw  carrrier  was  developed 
on  the  output  screen  of  the  F4108  intensifier.  A  37  MHz 
carrier  obtained  from  the  signal  generator  was  applied  to 
both  the  delay  line,  an  Isomet  1201  acousto-optic  modulator, 
and  the  laser  diode  light  source.  The  electrical  input  power 
to  the  delay  line  was  1  watt.  The  current  modulation  on  the 
laser  diode  was  approximately  50%. 

The  moving  phase  disturbances  in  the  delay  line  are 
optically  projected  on  the  photocathode  of  the  image  intensi¬ 
fier  as  a  moving  image.  The  time  modulation  favors  a  partic¬ 
ular  spatial  phase  position  of  this  image  depending  on  the 
phase  relation  between  the  laser  output  and  delay  line  exci¬ 
tation.  The  phase  modulation  is  converted  to  amplitude  modu¬ 
lation  by  virtue  of  the  delay  line  operating  in  a  Bragg  mode. 
When  the  delay  line  is  carefully  adjusted  in  angle  with  re¬ 
spect  to  the  input  optical  beam  at  the  Bragg  angle,  only  one 
important  diffraction  beam  is  developed  (a  first  order). 
Essentially  only  this  one  diffracted  order  and  the  zero  order 
exist  in  the  output  space.  The  interference  between  these 
two  orders  develops  the  desired  carrier  amplitude  modulation. 
If  the  delay  line  were  in  the  Raman-Nath  regime,  with  normal 
input  incidence  of  the  illuminating  light  beam,  a  symmetrical 
pattern  of  orders  would  be  produced,  and,  most  importantly, 
these  would  be  the  -1,  0,  and  +1  orders.  The  phases  of  these 
orders  would  be  arranged  in  accordance  with  optical  phase 
modulation.  One  way  to  produce  carrier  amplitude  modulation 
from  this  situation  is  the  elimination  of  one  of  the  first 
orders,  leaving  only  the  other  first  order  and  the  zero  or¬ 
der.  This  is  then  equivalent  to  the  Bragg  effect. 

The  resulting  spatial  pattern  obtained  on  the  output 
screen  of  the  intensifier  for  both  inputs  to  the  delay  line 
and  the  laser  modulation  coming  from  a  37  MHz  sine-wave  gen¬ 
erator  is  shown  in  Figure  3-6.  This  pattern  represents  the 
auto-correlation  of  the  37  MHz  input.  The  non-uniformity 
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Figure  3-6.  Spatial  Correlation  Function  for  a  Simple  R.F. 

Carrier  (Autocorrelation) 


of  the  light  pattern  is  due  to  the  modal  structure  of  the 
laser  beam.  This  is  correctable  with  a  pinhole-spatial  fil¬ 
ter  in  front  of  the  laser  output. 

The  grating  pattern  disappears,  as  it  should,  when  ei¬ 
ther  the  laser  modulation  drive,  or  the  input  to  the  acoustic 
transducer  on  the  delay  line  is  removed.  If  the  delay  line 
is  moved,  the  pattern  moves  a  corresponding  distance. 

It  is  clear  that  if  a  more  complex  waveform  were  applied 
to  the  laser  and  delay  line  the  appropriate  carrier-borne 
auto-correlation  function  would  appear  as  a  spatial  pattern 
on  the  image  intensifier. 

3.2  SPACE- INTEGRATING  CORRELATOR 

3.2.1  General 

The  second  optical  section  of  the  AOAP,  from  the  image 
intensifier  to  the  output  detector  is  the  space  integrating 
correlator,  in  which  the  input  to  this  section  is  an  optical 
spatial  pattern,  and  this  is  imaged  to  a  delay  line.  As  the 
acoustic-borne  function  travels  by  the  pattern  image,  a 
point-by-point  spatial  multiplication  is  performed  at  the 
delay  line.  This  "continuum"  of  products  is  collected  by  the 
final  lens  and  the  spatial  integral  of  this  product  is  then 
delivered  to  the  final  detector  by  action  of  the  final  lens. 
Thus  the  cross-correlation  between  the  spatial  pattern  on  the 
intensifier,  and  the  function  in  the  delay  line  is  imple¬ 
mented. 


Our  original  concept  was  to  use  the  incoherent  optical 
field  generated  by  the  image  intensifier  as  the  source  for 
the  space-integrating  correlator.  The  optical  radiation  from 
the  image  intensifier,  which  is  self-luminous  radiation  from 
a  phosphor,  is  totally  incoherent,  and  radiates  into  an  angu¬ 
lar  hemisphere.  Of  this,  only  a  small  sector  of  ray  direc¬ 
tion  can  be  used.  Rays  in  this  sector  are  sufficiently  well 
collimated  so  as  to  essentially  register  in  angle  with  the 
acoustic  wavefront  in  the  delay  line.  Rays  which  cross 
acoustic  wavefronts  receive  very  little  net  modulation  in 
optical  phase.  Rays  which  are  perfectly  aligned  with  the 
acoustic  wavefronts  receive  consistent  modulation  during 
their  entire  travel  through  the  delay  line.  This  will  be 
quantitatively  discussed  in  a  later  section. 

If  the  acoustic  wavefronts  have  significant  curvature, 
it  may  be  that  the  crossing  of  wavefronts  is  unavoidable,  and 
poor  modulation  will  result  for  all  rays. 

3.2.3  Delay  Line 

The  incoherent  mode  was  attempted  with  a  delay  line 
built  by  Andersen  Laboratories.  This  is  a  fused  quartz  line 
with  a  velocity  of  approximately  6mm/ysec.  Thus  at  60  MHz,  a 
spatial  frequency  is  generated  which  is  compatible  with  10 
cycles/mm  produced  on  the  image  intensifier.  The  delay  line 
was  tested  and  shown  by  Andersen  to  be  capable  of  43%  dif¬ 
fraction  efficiency  in  the  primary  Bragg  1st  order.  This 
line  is  15  mm  thick,  which  classifies  it  as  weakly  Bragg. 

The  43%  efficiency  implies  total  optical  phase  modulation  of 
approximately  1.43  radians,  which  is  to  say  that  the  peak 
optical  phase  deviation  caused  by  a  negative  or  positive 
swing  of  the  acoustic  carrier  corresponds  to  1.43  radians. 
This  result  assumes  strictly  Bragg  operation. 

3.2.4  Optics 

The  optics  used  for  the  space  integrating  correlator  is 
similar  to  those  used  for  the  time  integrating  correlator.  A 
telecentric  system  using  two  Bausch  and  Lomb  Baltar  f/2.8 
152mm  lenses  was  used,  although  this  is  not  explicitly  shown 
in  Figure  2-1.  This  system  images  the  input  pattern  to  the 
delay  line.  The  telecentric  system  provides  the  equivalent 
of  a  spatial  frequency  plane  between  the  two  lenses,  although 
a  spatial  frequency  is  not  displayed  here  in  an  incoherent 
system.  Position  in  this  plane  corresponds  to  ray  angle  in 
either  the  object  (grating)  or  image  (delay  line)  planes. 

Thus  the  correct  sector  of  angles  can  be  selected  here  to 
properly  register  with  the  acoustic  wavefronts  in  the  delay 


line. 


The  optical  field  produced  by  the  intensifier  may  be 
adequately  simulated  by  a  photographic  grating  transparency 
with  a  10  cycle/mm  pattern  on  it.  With  incoherent  illumina¬ 
tion  through  a  diffuse  screen,  the  radiation  obtained  is 
dispersed  into  approximately  a  hemisphere,  as  it  is  in  the 
image  intensifier.  The  light  level  may  be  easily  matched  to 
simulate  that  of  the  intensifier.  By  using  a  spectral 
source,  such  as  a  mercury  vapor  lamp,  a  spectral  optical 
width  roughly  equivalent  to  the  intensifier  output  may  also 
be  obtained. 

3.2.5  Detector 

An  avalanche  photodetector  (APD)  diode  was  used  as  the 
output  detector  in  the  space-integrating  correlator.  This 
photodetector  was  used  with  a  broad-band  three-stage  transis¬ 
tor  amplifier  which  had  current  and  voltage  limiting  circuits 
to  protect  the  diode.  The  APD  has  voltage,  current  and  power 
limits  that  must  be  strictly  adhered  to  in  order  to  avoid 
damaging  it. 

A  constant-current  mode  was  selected  as  an  alternative 
to  temperature  compensating  the  applied  diode  voltage.  The 
diode  itself  is  an  RCA  C30817  device,  with  a  rated  voltage 
dependence  with  temperature  of  1/2  volt  per  °C.  This  temper¬ 
ature  dependence  is  fairly  linear  between  -40 °C  and  +60 °C. 

The  circuitry  thus  allows  one  to  select  the  average  diode 
current  over  fairly  wide  limits  of  temperature  and  average 
input  light  intensity. 

3.2.6  Space  Integrating  Correlator  Experimental  Results 
3.2.6. 1  Incoherent  mode 

In  the  incoherent  mode,  the  frequency  plane  elimination 
of  all  but  two  diffracted  orders,  which  was  the  scheme  used 
in  the  time  integrating  correlator,  cannot  be  used.  Such 
orders  are  not  separated  in  the  frequency  plane  since  the 
input  radiation  is  directed  in  a  wide  range  of  angles  due  to 
its  incoherent  nature.  An  alternative  scheme  is  to  convert 
the  optical  phase  modulation  to  amplitude  modulation  in  the 
object  plane.  This  is  accomplished  by  the  use  of  crossed 
polarizers  around  the  delay  line.  Then  in  the  absence  of 
acoustic  stress,  the  polarizers  cause  optical  ccinction  (in 
the  case  of  an  extreme  adjustment).  In  the  presence  of 
stress,  the  optical  material  becomes  biref ringent ,  which 
alters  the  polarization  state  and  results  in  the  passage  of 
light  through  the  second  polarizer.  The  result  is  an  ampli¬ 
tude  pattern  of  light  which  has  as  a  factor  the  acoustic 
amplitude  function  in  the  delay  line. 


A  calculation  shows  that  a  positive  signal-to-quantum 
noise  performance  should  be  obtainable  with  this  configura¬ 
tion.  In  the  experiment,  no  signal  could  be  recovered,  even 
with  the  achievement  of  considerable  integration  gain,  using 
a  spectrum  analyzer  on  the  output  waveform. 

It  is  felt  that  this  lack  of  output  is  due  to  a  combina¬ 
tion  of  non-optimum  modulation  of  the  delay  line,  and  poor 
registration  between  acoustic  wavefronts  and  light  rays. 

This  misregistration  may  be  attributable  to  curvature  in  the 
acoustic  wavefronts. 

3. 2. 6. 2  Coherent  Experiment 

As  an  alternative  to  the  foregoing  experiment  a  coherent 
space-integrating  correlator  was  implemented.  This  merely 
requires  the  substitution  of  coherent  (laser)  illumination 
for  the  incoherent  source  behind  the  input  photographic  grat¬ 
ing  transparency.  This  now  causes  separation  of  diffraction 
orders  in  the  spatial  frequency  plane  as  in  the  first  time- 
integrating  section,  so  that  the  use  of  crossed  polarizers  is 
not  required. 

A  Spectra-Physics  Model  125A  60  mw  laser  was  used.  Its 
1  mm  diameter  beam  was  focussed  into  a  25  micron  pinhole  and 
the  resulting  expanded  beam  was  collimated  with  a  162  mm 
lens.  This  beam  was  used  to  illuminate  the  grating  transpar¬ 
ency. 


The  resulting  spatial  pattern  was  imaged  on  a  second 
delay  line,  which  was  a  60  MHz  fused  quartz  line.  With  a 
15  mm  interaction  length,  this  line  is  again  a  Bragg  line. 

The  light  from  this  delay  line  was  collected  by  a  final  lens 
and  delivered  to  an  avalanche  silicon  photodiode.  In  Figure 
3-7  is  shown  a  series  of  photographs  of  the  r.f.  sine  wave 
recovered  from  the  detector.  The  oscilloscope  was  synchro¬ 
nized  to  the  r.f.  wave  applied  to  the  delay  line.  Each  of 
the  traces  shown  represents  a  different  positional  setting  of 
the  input  grating,  in  which  the  grating  was  moved  a  physical 
distance  corresponding  to  1/3  spatial  cycle.  The  proper 
corresponding  change  in  temporal  phase  is  demonstrated  in 
Figure  3-7. 

3. 2. 6. 3  Summary 

Experimental  results  have  been  shown  which  demonstrate 
the  successful  operation  of  a  time-integrating  correlator,  at 
a  wavelength  of  830  nm.  The  operation  of  a  coherent  second 
section  correlator,  which  is  required  for  an  adaptive  system 
is  demonstrated  at  633  nm.  The  required  marriage  of  the  two 
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Bottom  Traces:  Reference  oscillator  waveform. 
Top  Traces:  Detector  signals. 


Figure  3-7.  Time  Signal  Recoveries  from  Coherent  Output 
Correlator  for  Simple  R.F.  Carrier 


correlators  has  not  as  yet  been  satisfactorily  accomplished, 
due  to  the  poor  performance  of  the  incoherent  second  section. 
It  is  expected  that  with  the  addition  of  a  coherent  optical 
modulator,  which  will  allow  the  coherent  mode  in  the  second 
correlator  as  well  as  the  first,  a  complete  adaptive  loop  can 
then  be  demonstrated. 


4  ANALYSIS  OF  COMPONENT  PERFORMANCE 


DELAY  LINE  PERFORMANCE 


Introduction  to  Bragg  Diffraction 


We  now  take  up  the  theory  of  modulation  of  light  with  an 
AO  delay  line.  One  needs  to  know  the  amount  of  optical  modu¬ 
lation  achieved  as  a  function  of  the  total  phase  modulation 
provided  by  the  AO  line.  Further,  the  relation  between 
acoustic  power  in  the  AO  line  and  the  resulting  phase  modula¬ 
tion  is  also  needed. 

In  evaluating  the  performance  of  diffraction  of  light  by 
a  delay  line  in  the  Bragg  mode,  a  fundamental  variable  is  the 
total  phase  shift  imposed  on  a  light  ray  assuming  that  it 
travels  entirely  through  the  line,  perpendicular  to  the  di¬ 
rection  of  acoustic  propagation,  as  shown  in  Fig  re  4-1. 
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Figure  4-1 


A-0  Delay  Line  as  Phase  Modulator  of 
Perpendicular  Light  Beam 
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We  define  this  total  phase  shift  according  to  the  value 
of  the  peak  amplitude  on  the  positive  and  negative  swings  of 
the  sinusoidal  variation: 


where  An  is  the  peak  magnitude  of  the  optical  index  change 
caused  by  the  acoustic  wave,  and  L  is  the  AO  interaction 
length  as  shown  in  Figure  4-1.  The  variation  of  4>  at  any 
point  along  the  delay  line  is  given  by 

x  ±  vat 

4>  =  4>0  +  4>msin2n  ( - - - )  (4-2) 

A 

where  A  =  acoustic  wavelength. 

In  Bragg  diffraction,  incoming  light  rays  do  not  follow 
an  acoustic  wavefront,  but  rather  are  correctly  oriented  at 
the  "Bragg  angle",  which  is  given  by 

—  1  ^ 

eR  =  sin  1  -  (4-3) 

B  2A 

In  this  position,  the  diffracted  first  order  is  oriented 
at  -00,  and  the  geometry  is  then  such  that  there  is  a  con¬ 
tinuous  coherent  reinforcement  of  one  wave  by  the  other. 

This  is  caused  by  the  volume  grating  pattern  corresponding  to 
the  index  variation  induced  by  the  acoustic  wave.  This  geom¬ 
etry  is  illustrated  in  Figure  4-6. 

4.1.2  Photo-Elastic  Effect. 


To  understand  this  type  of  acoustic  modulation  one  must 
consider  the  relation  between  mechanical  stress,  strain,  and 
the  resulting  optical  photo-elastic  effect.  The  optical 
behavior  of  a  stressed  solid  is  very  similar  to  the  behavior 
of  the  optical  behavior  experienced  in  a  crystal  [9],  [10], 

In  the  general  case,  i.e.  anisotropic  crystals,  the 
optical  behavior  varies  between  the  principal  axes  of  such  a 
crystal.  The  most  general  crystal  is  the  biaxial  crystal  in 
which  the  behavior  is  different  between  three  mutually  per¬ 
pendicular  principal  axes.  A  somewhat  simpler  case  is  the 
uniaxial  crystal  in  which  one  axis  behaves  differently  than 
the  other  two,  which  are  isotropic  between  themselves.  The 
simplest  <"ase  is  the  isotropic  crystal  in  which  behavior  is 
uniform  in  any  direction. 

The  uniaxial  case  is  of  the  most  interest  to  us,  since 
when  an  isotropic  material  is  stressed  with  a  simple  compres¬ 
sive  (or  tensile)  stress  in  some  direction,  the  isotropic 
material  becomes  uniaxial  in  its  optical  behavior,  with  its 
principle  uniquely  altered  axis  (sometimes  referred  to  as  the 
extraordinary  index  axis)  coincident  with  the  stress  vector. 


In  a  uniaxial  crystal/  it  is  convenient  to  define  an 
ellipsoid  of  wave  normals/  or  "optical  indicatrix".  This  is 
given  by  the  equation:  [9] 


=  E  •  D  =  Electric  Energy  Density 


(4-4) 

where  Dx,  Dy,  and  Dz  are  electric  displacement/  or  "D" 
vectors.  This  is  the  equation  of  an  ellipsoid,  and,  it  can 
be  shown  from  the  Maxwell  theory,  that  the  principal  direc¬ 
tions  of  the  D-vector  associated  with  any  wave  normal  can  be 
obtained  from  the  construction  in  Figure  4-2. 
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Figure  4-2.  Principal  Axes  of  D-Vector  as  Obtained  from 

Ellipsoid  of  Wave  Normals 


For  any  optical  wave  normal  S  (not  necessarily  the  ray 
direction)  the  origin  plane  perpendicular  to  S  intersects  the 
ellipsoid  of  wave  normals  in  a  plane  ellipse,  shown  in  dotted 
lines  in  Figure  4-2.  The  minimum  and  maximum  radii  of  this 
ellipse,  i.e.  the  principal  axes,  define^the  two  perpendicu¬ 
lar  principal  D-vectors  associated  with  S. 

In  particular,  if  If  is  coincident  with  one  of  the  prin¬ 
cipal  axes,  say  _the  z-axis  for  example,  the  two  D-vectors 
associated  with  s'  are  the  other  two  principal  axes,  the  x- 
and  y-axes  in  this  case.  Thus  in  this  case  of  principal 
stress-optical  alignment,  the  D  and  E  vectors  are  colinear, 
and  are  related  by  e  x  in  the  x-directior.  and  ey  in  the  y- 
direction.  It  follows  that  the  phase  velocities  Vp  of 
the  two  waves  propagating  in  the  z-direction  are  1  // u e x’  for 
eT  in  the  x-direction,  and  1  // y e y k  for  E  in  the  y-direction. 


The  appropriate  axis  is  selected  from  the  ellipsoid 
according  to  the  E-vector,  since  in  a  dielectric  =  yQ, 

B  =  y qH  in  all  directions.  Only  the  relation  of  if  and  If 
are  affected  by  the  changes  in  the  dielectric  constant  e . 


Thus,  referring  to  Figure  4-1  if  an  incident  ray  travel¬ 
ing  in  the  z-direction  is  incident  on  a  medium  stressed  (or 
strained)  in  the  x-direction,  two  results  are  possible. 

1 )  the  electric  vector  is  oriented  along  the  y-axis, 
in  which  case 

1 

vp  =  . . .  (4“5> 

V  uOE0RD 

where  e0RD  =  ey  =  ez 

and  the  ray  is  called  the  "ordinary  ray". 

2)  the  electric  vector  is  oriented  along  the  x-axis  in 
which  case 


VP 


/  u  e 
M  o  x 


(4-6) 


and  the  ray  is  called  the  "extraordinary  ray”. 


The  extraordinary  ray  thus  has  a  different  velocity  than 
with  rays  oriented  with  E-vectors  perpendicular  to  the  x- 
axis. 


In  the  case  of  a  general  ray  direction  with  a  general  E- 
orientation,  a  uniaxial  medium  refracts  two  rays,  according 
to  the  theory  previously  described.  The  travel  of  these  two 
rays  are  in  slightly  different  directions  in  this  general 
case.  The  important  consequence  to  us  is  that  the  phase 
delay  of  the  two  rays  differ,  that  of  the  ordinary  ray  and 
that  of  the  extraordinary  ray. 

The  equation  of  the  ellipsoid  of  wave  normals,  equation 
(4-4),  is  written  for  the  principal  optical  axes  and  the 
x,y,z  coordinate  system  being  co-aligned.  The  coefficients 
of  this  equation  are  ax  =  1  /ex,  ay  =  1/£y,  and  az  = 

1/e?.  In  the  case  of  a  biaxial  crystal  in  a  general  orien¬ 
tation,  with  a  general  mechanical  stress  applied,  the  change 
in  the  a^j's  =  1/eij's  is  given  by 


Aaij  "  <*ijkl  Tkl 


(The  RH  side  is  summed 
over  k,l) 


(4-7) 


where 

Aa^j  =  dielectric  impermeability  tensor 

(2nd  rank) 

Tkl  =  stress  tensor  (2nd  rank) 

9ijkl  =  stress-optical  tensor  (4th  rank) 

Each  index  covers  three  spatial  dimensions.  Thus  there 
are  nine  terms  in  Aaij  and  T^i,  and  81  terms  in 
‘Hjkl* 

However,  the  cross  terms  of  A  a  and  T  are  symmetric,  and 
so  the  9  terms  of  each  reduces  to  six.  The  terms  of 
Tki,  f°r  example,  are  given  by 

fp  rp  m  fn  rp  rp 

xx'  yy'  zz'  Ayz'  xz'  xy* 

qijki  now  reduces  to  36  independent  terms. 

In  the  case  of  a  material  which  is  isotropic  when  un¬ 
stressed,  it  can  be  shown  that  this  reduces  to  2  independent 
coefficients.  We  then  have,  by  orienting  the  coordinates 
along  the  principal  axes  of  the  stress  tensor,  and  thereby 
eliminating  the  cross  (shear)  terms:  [ 9 ] 
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Coefficients  can  also  be  defined  which  relate  the  change  of 
ai  =  1/e  j_  to  the  component  of  strain  instead  of  stress. 

Thus  we  have: 

Sii  =  Strain  in  itn  coordinate 


Aax 

=  Pi  1 

sxx 

+ 

Pi  2 

>i 

>i 

to 
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Aay 
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Referring  to  Figure  4-3  we  now  consider  the  special 
simple  case  which  is  of  primary  interest  to  us  here.  If  the 
optical  wave  is  propagating  in  the  z-direction,  in  an  isotro¬ 
pic  medium  which  is  stressed  in  the  x-direction,  then  the 
medium  will  act  like  a  uniaxial  crystal  with  the  optical 
axis,  (direction  of  extraordinary  index  of  refraction)  in  the 
x-direction.  y 


ACOUSTIC 

AXIS 

Figure  4-3.  Geometry  of  A-0  Delay  Line  and 
Longitudinal  Stress,  Txx,  and  Corresponding 
Ellipsoid  of  Wave  Normals 

If  a  longitudinal  acoustic  wave  is  traveling  in  the  x- 
direction,  this  will  cause  x-oriented  stress,  alternating 
between  compression  and  dilation  according  to  the  acoustic 
wavefronts. 

It  can  be  shown  that  in  most  cases  when  an  isotropic 
material  is  under  dilatory  or  tensile  stress,  it  behaves  like 
a  "positive"  uniaxial  crystal,  such  as  crystalline  quartz, 
and  when  it  is  in  compressive  stress,  it  behaves  like  a 
"negative"  uniaxial  crystal  such  as  felspar. 

A  positive  crystal  has  an  ellipsoid  of  wave  normals  or 
indicatrix  similar  in  shape  to  that  shown  in  Figure  4-3,  in 
which  the  index  along  the  optical  x-axis  ne  is  larger  than 
along  the  other  axes  (y  or  z)  which  is  n0.  Hence  a  wave  in 
which  the  electric  vector  is  aligned  along  the  x-axis  will 
undergo  more  phase  delay  than  a  wave  in  which  the  E-vector  is 
pointed  in  the  y-direction,  in  the  case  of  tensile  stress. 
Thus  in  tension  the  Ex  wave  will  be  retarded  with  respect 
to  the  Ey  wave. 

Referring  to  the  first  two  of  equations  (4-9),  and  re¬ 
membering  that  the  a's  refer  to  1/e  =  1/n2 ,  we  see  that  in 
tension  (Aax-Aay)  is  negative  (i.e.  Anx>Any)  for 
tensile  stress  in  the  x-direction.  For  only  x-strain  these 
equations  reduce  to: 


Now,  the  general  effect  of  tension  is  to  cause  a  lower 
molecular  density,  which  lowers  the  index  of  refraction, 
raising  the  corresponding  a  coefficient.  Thus  as  defined  in 
(4-9),  the  p  coefficients  are  positive  if  sxx  represents 
tensile,  rather  than  compressive  strain.  If  sxx  is  thus 
positive  for  tension,  A ax  -  Aay  is  negative  for  tension, 
and  the  p-coef f icients  are  positive,  then  we  conclude  from 
(4-10)  that  the  value  of  p,  2  is  larger  than  pj j  .  The  p12 
coefficient  relates  optical  behavior  of  a  wave  in  which  the 
E-vector  is  parallel  to  the  y-axis,  with  strain  in  the  x- 
direction.  This  is  an  "ordinary"  wave  (the  extraordinary 
index  is  coaligned  with  the  stress).  The  Pj  1  coefficient,  by 
like  reasoning,  treats  the  "extraordinary  ray".  Thus  we 
conclude  that  the  phase  advance  (retardation)  for  a  tensile 
(compressive)  stress  is  greater  for  the  ordinary  ray  than  for 
the  extraordinary  ray. 

This  conclusion  is  verified  by  the  example  of  "stress- 
optical  coefficients"  associated  with  fused  quartz,  which 
are:  [10]  [11] 

q  =  -.65 

C2  =  -4.2 


The  stress-optical  coefficients  relate  index  change 
rather  than  the  a-coeff icients  directly  to  stress.  They  are 
like  the  q's  of  equations  (4-8)  in  the  sense  that  they  are 
based  on  stress  rather  than  strain.  Thus: 


A  ne  =  C1 T 


A  nQ  =  C2  T 


(4-11  ) 


The  coefficient  q  -  C2  is  positive,  and  so  we  see  that 
in  the  condition  of  tensile  stress,  in  an  isotropic  medium, 
ne  is  larger  than  nQ,  a  condition  of  a  positive  uniaxial 
crystal. 


The  quantity  (C1  -  C2 )  is  called  the  relative  stress- 
optical  coefficient.  By  multiplying  any  of  the  C’s  by  the 
modulus  of  elasticity  for  pure  extension/  Ea,  the  strain- 
optical  coefficients/  S1  and  S2  which  relate  index  change  to 
strain  can  be  obtained. 

The  value  of  Ea  which  is  used  in  these  calculations 
pertains  to  ultrasonic  wavefronts,  in  which  the  lateral  di¬ 
mension  of  the  wavefront  is  many  wavelengths.  Under  these 
conditions  a  plane  compression  is  not  accompanied  by  lateral 
expansion  (governed  by  Poisson's  ratio).  This  is  called  the 
bulk  modulus,  and  is  sometimes  denoted  as  k*.  In  terms  of 
the  Lame*  elastic  moduli,  y  (shearing  modulus)  and  X,  Ea 
(or  k)  is  given  as  X  +•  2y  whereas  Young’s  modulus  for  uncon¬ 
strained  lateral  dimensions  is  given  by  [12] 


.  3X  +  2y  . 
v  l — — - J 

X  +  y 


In  elasticity  theory,  X  accounts  for  volumetric  change, 
and  y  for  shearing  effects.  In  fact,  y  is  the  modulus  of 
shear  commonly  denoted  by  engineers  as  G.  In  the  case  of 
pure  extension,  with  no  lateral  strain,  the  extension  is 
equivalent  to  2  shears  equal  to  the  longitudinal  strain 
s  x  at  45°,  one  in  the  x-y  plane,  and  one  in  the  x-z 
plane,  and  a  volume  change  equivalent  to  the  longitudinal 
strain  sx.  Thus  the  total  modulus  is  given  by  X  +  2y . 


From  equation  (4-4),  and  the  definition  of  index  of 
refraction  n  =  /T',  e  is  the  relative  dielectric  constant,  the 
Aa^’s  can  be  related  to  the  An’s  as  follows: 


Since 


1  1 


we  have: 


(4-12) 


(4-13) 
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Where  the  change  A  a  is  related  to  strain  by  p,  and  A  a  is 
related  to  stress  by  q,  we  have  defined  coefficients  S  and  C 
which  relate  the  change  in  index  of  refraction  An  to  strain, 
s,  and  stress,  T  respectively. 

Hence: 


Si  =  - 


T  Pi 


(4-14) 


ci  =  -  2  qi  i 

The  interrelations  between  the  p's,  q's,  C's,  and  S's 
are  thus  found  to  be  as  follows: 


TABLE  III.  CONVERSIONS  BETWEEN  PHOTO-ELASTIC  PARAMETERS 
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4.1.3  Acousto-Optic  Interaction 

The  preceding  theory  relates  the  changes  in  the  dielec¬ 
tric  parameters,  and  specifically  the  index  of  refraction,  in 
a  solid  material  to  the  stress  and  strain  pattern  induced  in 
that  solid  by  an  acoustic  wave.  The  next  issue  is  to  summar¬ 
ize  the  resulting  optical  diffraction  behavior  in  the  Raman- 
Nath  and  Bragg  modes  which  is  caused  by  the  spatially  period¬ 
ic  dielectric  variations  that  accompany  such  an  acoustic 
wave. 


4. 1.3.1  Raman  Nath  Diffraction 

Assume  that  an  acoustic  line  is  thin  enough  in  the  di¬ 
rection  of  optical  propagation  so  that  a  light  ray  essential¬ 
ly  stays  registered  with  an  acoustic  wavefront  during  its 
entire  transit  through  the  delay  line.  The  physical  transit 
length  which  is  involved  has  been  denoted  by  L  in  the  pre¬ 
vious  discussion.  The  optical  behavior  of  this  configuration 
is  commonly  referred  to  as  Raman-Nath  diffraction.  In  this 
case,  a  differential  phase  shift  is  imposed  on  the  optical 
ray*  ♦in/  which  is  a  coherent  accumulation  of  phase  shift 
caused  by  the  index  of  refraction  change  over  the  interaction 
length  L.  Thus  the  value  of  the  index  change  An  which  is 
experienced  by  the  ray  is  constant  over  L.  The  resulting 
phase  shift  is  then  given  by  <t>m  =  An  k  L 

♦  m  ~  A n  k  L  (4-15) 


where 


k  =  free  space  wave  number  =  2n/X0. 

If  the  acoustic  propagation  is  along  the  x-axis,  the 
spatial  index  variation  is  then  in  terms  of  x.  If  we  assume 
that  a  simple  cw  sinusoidal  carrier  exists  of  frequency 
<*>a  =  2*  va/A ,  resulting  in  a  spatial  variation  of 
n  =  nQ  +  An  cos  2"  x/A ,  then  this  spatial  variation 
will  be  imposed  as  a  spatial  phase  modulation  on  a 
collimated  input  beam  passing  through  the  delay  line. 

If  this  input  light  beam  is  given  functionally  as 
E0  e3 (kz-w0t ) ,  where  the  E-vector  is  assumed 
to  be  aligned  in  the  x-y  plane,  and  u>0  is  the  optical 
frequency,  the  light  emerging  from  the  delay  line  will  be 
given  by: 


E  =  Eq  e 


[j(kz-o>  t)  +  jkL  An  cos  2 it 


(x-vat) 


(4-16) 
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We  now  note  that  e^  008  can  be 
mathematically  expanded  into  the  following  series: 

e3*"  CO.  •  .  f  jn  Jn  ,*m)  ejn6 


(4-17) 


Using  this  expansion,  we  find  that  the  optical  wave 
emerging  from  the  delay  line  can  then  be  described  by  the 
following  expression:  [13] 


[j(kz-ut)  +  j2*n  (■ 


x-vat 


(4-18) 


Each  term  in  the  summation  corresponds  to  an  optical 
diffracted  order.  In  particular,  the  first  order  of 
diffraction  has  a  coherent  amplitude  which  is  given  by: 


E1  *  (kLAn) 


(4-19) 


The  behavior  described  by  these  expressions  is  Raman- 
Nath  diffraction,  and  is  for  relatively  thin  lines.  For 
Raman-Nath  theory  to  apply,  the  diffracted  ray  must  stay 
essentially  registered  with  the  acoustic  wavefront,  and  not 
cross  over  to  adjacent  wavefronts.  The  situation  is  shown  in 
Figure  4-4. 
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Figure  4-4.  Registration  of  Light  Rays  with  Acoustic 

Wavefronts 


The  angle  between  the  zero  order  and  the  first 
diffracted  order  is  given  by  X/nA.  Hence  the  linear  spread 
between  the  two  orders  after  passing  through  L  is  XL/nA. 
This  should  not  exceed  1/4  acoustic  wavelength.  Hence  for 
Raman-Nath  operation: 


(4-20) 
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other  words ,  the  interference  between  successive  diffraction 
contributions  as  one  proceeds  through  the  interaction  length 
is  coherently  reconstructive.  The  interference  effect  for 
the  other  orders  is  destructive,  and  hence  all  other  orders 
are  greatly  attenuated.  For  a  strong  enough  signal  and  re¬ 
sulting  strong  phase  modulation  essentially  100%  of  the  input 
beam  can  be  diffracted  into  the  first  order.  The  situation 
is  pictured  in  Figure  4-6. 
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Figure  4-6.  Geometry  of  Bragg  Diffraction  Showing  Alignment 
of  Acoustic  and  Optical  Wavefronts 
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In  Figure  4-6  we  indicate  a  thin  slice  of  delay  line 
medium  with  width  dL.  When  the  strong  zero  order  impinges  on 
this  thin  slice,  Raman-Nath  diffraction  will  occur  according 
to  the  preceding  theory.  The  resulting  array  of  orders  will 
appear  symmetrically  around  the  zero  order.  Let  the  ampli¬ 
tude  of  the  zero  order  be  denoted  by  c  (for  carrier),  and  the 
first  order  denoted  by  £  (for  sidebancT).  The  following  argu¬ 
ment  is  taken  from  reference  [14]. 


When  the  zero  order  with  amplitude  c  passes  through  the 
thin  phase  grating,  the  complex  amplitude  of  the  emerging 
wave  is  given  by 
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(4-22) 


This  expression  accounts  for  Raman-Nath  diffraction  from 
the  thin  grating,  in  which  the  phase  modulation  involves  a 
small  angle  d<|>  =  k  dL  An.  The  last  two  terms  in  the  above 
expression  represent  the  positive  and  negative  first  optical 
orders.  Since  we  assume  that  only  the  positive  first  order 
has  reconstructive  interference,  this  order  will  build  up 
through  successive  thicknesses  dL,  and  the  other  order  will 
have  destructive  interference  and  will  not  appear  to  any 
significant  degree  in  the  emerging  wave.  This  constitutes 
the  Bragg  effect.  The  differential  buildup  of  the  diffracted 
order  s  can  be  expressed  in  terms  of  the  differential  accumu¬ 
lation  of  phase  shift  d$  from  equation  (4-22)  as 


c  ds  c 

ds  =  -  d<j>  =*>  -  =  —  (4-23) 

2  d*  2 

In  a  like  manner,  when  the  diffracted  order  s  builds  up 
to  a  significant  amplitude,  it  will  act  like  a  new  zero  order 
and  will,  by  Bragg  diffraction,  create  a  sideband  which  is 
oriented  exactly  with  the  original  zero  order,  and  will  co¬ 
herently  combine  with  the  original  zero  order.  This  argument 
does  not  answer  the  question  of  the  phase  relation  between 
the  original  zero  order  c,  and  the  sideband  of  the  diffracted 
order  which  is  combining  with  c.  This  problem  is  resolved 
indirectly  by  appealing  to  the  assumption  that  all  the  dif¬ 
fraction  in  the  Bragg  mode  goes  into  the  one  first  order,  and 
then  back  again  into  the  zero  order,  these  two  orders  being 
symmetrically  oriented  with  respect  to  the  acoustic  beam  at 
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i 6 b ,  the  Bragg  angle  (equation  4-3).  Under  this  assump¬ 
tion,  all  the  energy  associated  with  the  zero  order  entering 
the  interaction  region  at  the  left  hand  side  of  the  interac¬ 
tion  region  must  equal  the  sum  of  the  energies  of  the  two 
orders  at  any  arbitrary  point  within  the  line,  or  at  the 
plane  of  emergence  on  the  right  hand  side  as  shown  in  Figure 
4-6.  Thus  we  have: 

c2  +  s2  =  cQ2  (4-24) 


where 


cQ  =  amplitude  of  incident  zero  order  light  beam. 
Differentiating  with  respect  to  $ ,  we  have: 
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(4-25) 


Equations  (4-23)  and  (4-25)  can  be  separated  into  two 
homogeneous  equations  in  c  and  £  by  the  dif ferentat ion  by 
d/d$  producing 
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(4-26) 


The  solutions  to  these  equation  are: 
♦ 


C  =  Cq  cos 


.  ♦ 

s  =  c„  sin  — 


(4-27) 


The  relative  intensity  of  the  diffracted  order  to  the 
incident  intensity  is: 


I 
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=  sin 


2  1 


(4-28) 


4.1.4  Imaging  with  Bragg  Diffraction 


The  two  orders  developed  by  the  Bragg  diffraction  in  the 
delay  line  appear  collimated  and  at  a  relative  mutual  angle 
corresponding  to  the  acoustic  spatial  frequency  in  the  image 
plane.  Thus  the  relative  spatial  phase  of  the  two  orders 


varies  at  a  spatial  rate  which  is  the  acoustic  spatial  fre¬ 
quency.  The  coherent  sum  of  the  two  orders  can  then  be  ex¬ 
pressed  as: 


sum 
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(4-29) 
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8  =  acoustic  signal  phase 


(4-30) 


The  light  waves  are  detected  by  the  image  intensifier 
(or  other  detecting  device)  which  can  be  adequately  modeled 
as  a  square  law  detector.  The  light  intensity  is  given  as: 


1  ,  .  jV  ,  “ je  a  , 

-  (c  +  }se  J  (c*  +  (3s)*e  ) 


(4-31  ) 


For  this  argument,  we  assume  that  c  and  s  are  the  real 
functions  derived  previously,  equation  (4-27).  Thus  we  have 


1  .  3ea 

— —  (  cc*  +  ss*  +  }c*s  e  -  }cs*  e  J 


co  ,  *m  ^m 

-  (.1  +  2  cos  -  sin  -  sin  8  _J 

2  2  2  a 


Hence: 


I 

1  +  sin  4> m  sin  8g 

c  2 
co 

2 

(4-32) 


In  the  above  expression,  the  presence  of  the  factor  sin 
9.  tells  us  that  the  intensity  varies  as  a  linear  function 
of  the  acoustic  carrier,  and  this  variation  or  carrier  is 
multiplied  by  sin  4ro.  The  function  <|>m  is  here  regarded 
as  the  signal  amplitude  (total  phase  shift  across  the  Bragg 
line  due  to  the  photoelastic  effect).  If  <f>m  is  appreciably 
less  than  1  radian,  the  modulation  amplitude  becomes  linear 
with  signal  amplitude  as  desired.  This  variation  appears 
around  a  bias  term  equal  to  1/2. 


The  maximum  and  minimum  values  of  I  occur  at  values  of 
9a  equal  to  (2n  +  1)  w/2  where  sin  6a  =  ±1.  Thus: 


1  +  sin  <f>, 


1  -  sin 


(4-33) 


The  usual  definition  of  intensity  contrast  ratio  is  thus 
given  as: 


Contrast  Ratio  = 


rmax  “  Imin 
*max  +  ^-min 


=  sm  <j>. 


(4-34) 


so  that  the  contrast  ratio  is  also  a  linear  function  of  <f>m 
for  small  values  of  <j>m. 

4.1.5  Acoustic  Wave  Propagation 

We  now  address  the  relation  between  acoustic  power  and 
the  resulting  phase  modulation.  For  the  simplest  case ,  let 
us  assume  that  a  longitudinal  acoustic  wave  exists  in  the 
delay  line,  traveling  in  the  x-direction,  as  shown  in  Figure 
4-7.  A  thin  slice  of  medium  is  indicated  with  width  dx.  The 
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Figure  4-7.  Acoustic  Longitudinal  Wave  Propagation 
Characterized  by  Stress,  T,  and  Particle  Velocity,  v  on  Faces 

of  Thin  Slice  of  Delay  Medium 
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acoustic  wave  is  now  characterized  by  the  acoustic  sinusoidal 
pressure  on  the  plane  boundaries  of  the  slice  dx,  and  the 
particle  velocity  varying  sinusoidally  in  the  x  direction. 

The  differential  force  on  the  slice  will  accelerate  the  par¬ 
ticles  lying  within.  Thus  we  can  write: 


T2  -  T1 

3Txx  3  v 

■  A  dx  -  =  p  A  dx  - 

9  x  9  t 

3  TXX  9  V 

9  x  P  3  t 

(4-35) 

where 

Txx 

=  acoustic  pressure  (stress) 

Tl  ,T2 

=  two  values  of  pressure  on  the  boundaries 
of  dx 

P 

=  density  of  medium 

V 

=  particle  velocity 

A 

=  face  area  of  slice. 

Similarly,  we  argue  that  a  space  gradient 
will  cause  a  time  rate  of  change  of  strain,  s. 

in  velocity 
We  have: 

v2-vi  = 

3  v  9  S  9  , T„„  . 

-  dx  =  -  dx  =  —  ( -**■  )  dx 

9  x  9 1  9t  Ea 

3  v  1  3  Txx 

(4-36) 

9  x  ~  Ea  9  t 

where  Ea  is  the  appropriate  acoustic  modulus  of  elasticity, 
given  by  the  linear  combination  of  Lame'  coefficients  X  +  2u 
as  explained  in  Section  4.1.2. 


Combining  equations  (4-35)  and  (4-36)  we  have  the  second 
order  equation: 


The  most  general  solution  to  this  equation  is 
Txx  5  Txx  (x  1  vafc) 


where 


pva2  =  Ea 


4.1.6  Acoustic  Energy  and  Power 


(4-38) 


(4-39) 


The  mechanical  energy  density  stored  in  any  microscopic 
region  of  the  medium  is  given  by  1/2  the  product  of  acoustic 
pressure  and  strain.  Thus: 


T  S 
Am  m 


w  =  Energy  Density  = 


T™  =  Peak  value  of  T, 


sm  =  Peak  value  of  strain  (sxx) 


E  _  s_/ 
a  m 


(4-40) 


The  acoustic  power  is  the  time  rate  of  flow  of  acoustic 
energy.  Thus  the  acoustic  intensity  (power  density)  is  the 
acoustic  velocity  times  the  energy  density.  The  result  is 
expanded  using  (4-39)  and  (4-40): 


Ia  =  wa  va  *  J  sm  p  va 

From  equation  (4-14)  and  Table  III  we  have: 


(4-41  ) 


An-  ■  S  •  Sm  ■ 
i  i  m 


Pi  nr 


(4-42) 


Referring  to  equation  (4-1)  we  have  a  relation  between 
the  peak  phase  shift  and  the  peak  shift  in  index  of  refrac¬ 
tion.  Combining  this  with  the  above  produces: 


2*L  pi  ni 


i. 

p  vaJ 


(4-43) 


(4-44) 


All  of  the  quantities  contained  within  the  brackets  of 
the  above  are  parameters  of  the  medium  material.  These  are 
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popularly  combined  into  a  material  "figure  of  merit."  This 
particular  mix  of  parameter  values  is  known  as  M2.  If*  in 
the  above  derivation,  instead  of  solving  for  <j>m,  we  solve 
for  the  acoustic  intensity  Ia,  there  results:  [14] 


Pa  1  2  p  va 

Ia  .  —  -  -  X  4(Ant)2  X 


_  2  „  6 
Pi  ni 


PA  *  2  X2 

—  -  2  C-&)  - 

*  '  _  '  a 


2ir  ‘ 


L^Mr 


M2  = 


Pi2  ni6 


p  v. 


It  will  be  recalled  that  from  (4-34), 


(4-45) 


»  sin”1  (contrast  ratio)  (4-46) 

so  that  for  imaging  with  Bragg  diffraction,  <j>m  may  be  taken 
as  the  contrast  ratio  for  values  of  this  quantity  substan¬ 
tially  less  than  unity. 


4.1.7  Power  Calculations  for  Various  Optimized  Geometries 

To  achieve  minimum  power,  for  a  given  diffraction  effi¬ 
ciency,  we  minimize  the  expression  for  P^.  This  is 
achieved  by  maximizing  L  and  minimizing  a.  There  is  no  real 
limitation  for  a  single  acoustic  frequency  and  an  unlimited 
incident  light  power.  However,  a  more  sensible  constraint 
results,  if  a  definite  acoustic  bandwidth  a f  is  considered. 
Such  a  bandwidth  implies  a  spatial  element  length  along  the 
acoustic  wave,  (along  the  y-direction  in  Figure  4-3)  given  as 


b  ■  -  Af  =  acoustic  bandwidth 

A  f 


(4-47) 


Assume  the  acoustic  beam  is  roughly  distributed  over  a  rect¬ 
angular  cross  section  (corresponding  to  a  rectangular 
acoustic  transducer),  of  dimensions  a  in  the  y-direction,  and 
L  in  the  z-direction  (interaction  length).  This  results  in 
an  acoustic  cross-section  area  of  aL,  and  equation  (4-45)  can 
be  written: 


x  2a 


2  (lm_) 

2*  LM 


(4-48  ) 
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The  next  question  involves  the  permissible  geometry  of 
the  Bragg  diffractor.  First,  there  is  a  limit  to  the  length 
of  the  interaction  region,  i.e.  the  dimension  L  in  the  z- 
direction.  This  is  because  the  incoming  light  and  the 
acoustic  beam  have  to  be  related  by  the  Bragg  angle,  which 
from  equation  (4-3)  is  seen  to  be  sensitive  to  the  acoustic 
frequency.  If  a  wide-band  acoustic  frequency  is  being  used 
for  diffraction,  this  implies  an  angular  spread  in  the  orien 
tation  angle.  This  spread  of  angle  must  be  within  the  range 
of  Bragg  compatibility.  A  convenient  way  of  looking  at  this 
problem  is  to  consider  the  far-field  divergence  of  both  the 
optical  angular  spectrum  diffracted  by  the  acoustic  signal, 
and  the  far-field  acoustic  divergence  of  the  acoustic  beam. 
For  valid  Bragg  diffraction  of  all  the  frequency  components 
of  the  signal  in  the  acoustic  beam,  the  far  field  divergence 
of  the  acoustic  beam  must  be  at  least  as  large  as  the  diver¬ 
gence  of  the  optical  diffracted  beam.  Thus  we  have: 

X 

A0B  *  Divergence  of  the  diffracted  optical  beam  =  — 

nb 

(4-49) 


where  b  *  spatial  dimension  (in  x)  of  acoustic  resolution 
element  as  limited  by  the  acoustic  bandwidth 
=  va/Af 

Af  *  acoustic  bandwidth. 

Let  the  far  field  divergence  of  the  acoustic  beam,  as 
limited  by  the  transmitting  aperture  dimension  L,  which  is 
also  the  interaction  length,  be  denoted  as  A0A  *  A/L. 

Hence: 


A 0  g  <  A0a  =» 


X 

-  < 

nb 


nbA 

L  <  - 

X 


n 


X  Af  fQ 


(4-50) 


The  next  geometrical  constraint  to  consider  is  the  ver¬ 
tical  extent  of  the  acoustic  beam,  a.  This  must  be  large 
enough  so  that  no  appreciable  acoustic  diffraction  spreading 
of  the  beam  occurs  in  this  direction,  which  would  cause 
curved  wavefronts,  and  would  cause  the  acoustic  beam  to  in¬ 
tercept  the  side-walls  of  the  delay  line.  Thus,  in  essence, 
we  want  the  spreading  to  be  less  than  the  original  size  of 
the  acoustic  beam.  Hence: 
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(4-51  ) 


Vertical  acoustic  beam  spread  =  -  <  a 

a 


=*>  a  >  /AW  =  /A  V-T  =  V_  /  — 

a  a  c 

la 


(4-52) 


where  W  is  the  length  of  the  delay  line  in  the  direction  of 
acoustic  propagation,  and 

A  =  acoustic  wavelength. 

Now  consider  the  case  of  a  simple  acoustic  modulator,  in 
which  no  specific  delay  time  is  required.  In  this  case  the 
optical  aperture  need  only  subtend  an  acoustic  resolution 
element,  b.  Now  equation  (4-51)  becomes 


vw 


<  a  =»  a  >  /bA  = 


(4-53) 


'Af  f. 


The  relaxation  of  the  acoustic  spreading  constraint  now 
means  that  the  optical  spreading  associated  with  a  becomes 
the  important  limitation.  If  we  demand  that  the  optical 
spreading  is  restricted  to  a  we  have: 


and  from  (4-50)  a  >  /  bA 


which  is  the  same  result  as  (4-53). 

The  limitation  on  the  ratio  a/L  is  now  given  as: 


a  /Af  fn 

a  o 

—  >  - 

L  tr  2 


\  /Af  f. 


(4-54) 
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A  Af  f. 


where  6  =  A f /fQ  *  fractional  bandwidth. 
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Substitution  of  this  limiting  result  in  equation  (4-48 } 
produces  the  minimum  acoustic  power  required  for  optimum 
delay  line  geometry. 


p 


2  m 


2* 


X  StiTtl 


n  v 


a 


(4-55) 


In  this  expression  we  see  that  a  new  grouping  of  the 
material  parameters  appears:  [14] 


n  va  M2  =  M3  (4-56) 

The  required  acoustic  power  for  a  given  amount  of  opti¬ 
cal  modulation  is  related  to  this  new  figure  of  merit,  M3, 
for  the  case  of  optimized  delay  line  geometry.  Using  this 
notation,  equation  (4-55)  can  be  written 


P 


a 


Af 


>  2 


M3  ST 


(4-57) 


where 


6  = 


Af 


i.e.,  the  fractional  acoustic  bandwidth. 


It  is  clear  that  minimum  drive  power  is  achieved,  for 
the  single  element  acoustic  modulator,  when  the  fractional 
bandwidth  is  made  its  maximum  value,  1.  We  will  take  this 
maximum  value  to  be  1 ,  although  theoretically  the  value  is  2 
for  a  base-band  situation.  If  6  =  1,  we  have  for  the  re¬ 
quired  power: 


>  2 


4  2 

fa 


(4-58) 


Now  consider  a  delay  line  modulator  in  which  a  total 
delay  time  T  is  required.  Equations  (4-50)  and  (4-52)  are 
used  to  evaluate  the  minimum  value  of  a/L  as  follows: 


53 


(4-59) 


Substituting  this  into  equation  (4-48)  we  have: 


—  >  2  ( —  •  — - 

A  f  2it  Mj  n  va 


/T.  2  il  /T 

1  2ir  ;  M- 


(4-60) 


where  6  *  fractional  bandwidth,  as  before,  and 

N  *  time-bandwidth  product  of  delay  line 
*  A  f  •  T 


As  before,  the  minimum  power  is  achieved  with  the  maxi¬ 
mum  fractional  bandwidth.  Thus  one  should  select  the  frac¬ 
tional  bandwidth  in  due  regard  for  the  avoidance  of  spectral 
folding  or  aliasing  in  the  active  band.  For  the  present 
experiment  a  fractional  bandwidth  of  1/6  was  chosen,  corre¬ 
sponding  to  a  required  bandwidth  of  10  MHz,  and  a  selected 
center  acoustic  frequency  of  60  MHz,  when  using  a  quartz  line 
with  an  acoustic  velocity  of  5.95  mm/sec.  Alternatively,  a 
frequency  of  37  MHz  was  used  with  a  flint  glass  line  with  a 
velocity  of  3.7  mm/sec.  Thus  in  either  case  a  spatial  fre¬ 
quency  of  10  cycles/mm  was  achieved. 

4.1.8  Delay  Line  Medium  Comparisons 


We  will  now  use  these  results  to  calculate  the  required 
acoustic  power  level  for  various  delay  line  materials.  For 
the  purposes  of  comparison,  the  maximum  modulation  level  will 
be  assumed,  providing  a  contrast  ratio  of  unity  from  the 
image  intensifier,  according  to  equation  (4-46).  As  seen 
from  this  relation  this  implies  a  peak  phase  deviation  of 
n/2.  A  reinforcement  of  this  conclusion  may  be  gained  from 
equation  (4-28).  From  this  relation,  it  is  seen  that  if  the 
peak  phase  deviation  is  i  /2,  half  the  intensity  (power)  of 
the  incoming  light  beam,  IQ,  is  transferred  to  the 


diffraction  order  I.  With  equal  amplitude  1st  order  and  zero 
order  light  beams  being  used  for  imaging  of  the  acoustic 
carrier,  total  nulls  will  be  developed  at  points  of  maximum 
interference,  and  amplitude  doubling  will  be  obtained  at 
points  of  opposite  phase  relation.  Thus  the  image  varies  in 
intensity  from  total  blackness  to  maximum  brightness,  corre¬ 
sponding  to  100%  modulation. 

Substituting  =  n/2  in  the  expression  for  a  simple 
modulator,  with  no  total  delay  requirement,  equation  (4-57), 
we  have  ( m/2ir ) 2  =  1/16,  and 


8M3  nr 


which  for  a  fractional  bandwidth  of  100%  becomes 


(4-61  ) 


For  the  delay  line  with  a  total  delay  time  requirement 
of  T,  corresponding  to  a  time-bandwidth  of  T  Af  =  N,  we  have, 
from  equation  (4-60) 


In  Table  IV  we  show  the  values  of  the  important  physical 
parameters  of  candidate  materials  for  the  delay  lines,  as 
pertains  to  acoustic-optical  performance,  index  of  refrac¬ 
tion,  acoustic  velocity,  and  the  figures  of  merit  M2  and 
M3,  in  MRS  units. 

We  repeat  the  definitions  of  the  figures  of  merit,  M2 
and  M3 


M3  =  n  va  M2 


where  p  *  photo-elastic  constant 
p  =  density. 


TABLE  IV.  PHYSICAL  PARAMETERS  OF  CANDIDATE  MATERIALS  FOR  AO 

DELAY  LINES ,  MRS  UNITS 


n 

va  meters/sec 

M2  sec3  /kg 

M3  sec2  -meter/kg 

Quartz 

1.46 

5950 

1.51  x 

10"15 

1.31  x  io-11 

Dense 

Flint 

Glass 

1 .92 

3100 

18.9  x 

10-15 

• 

to 

X 

0 

1 

LiNb03 

2.2 

6570 

6.95  x 

m 

7 

0 

1 

0 

X 

0 

• 

0 

PbMo04 

2.39 

3660 

35.8  x 

1  o-1 5 

31 .3  x  I0“n 

Te02 

2.27 

617 

793  x 

10“15 

111.0  x  10-11 

Water 

1.33 

1500 

158  x 

10“15 

31  .4  x  IQ”11 

In  Table  v  we  show  the  acoustic  attenuation  for  several 
candidate  materials.  These  are  given  for  the  proposed  param¬ 
eters  of  the  AOAP  system,  i.e.  a  center  frequency  of  60  MHz 
and  a  total  delay  of  7  nsec.  Also  tabulated  is  the  required 
length  for  7  nsec.,  and  the  acoustic  wavelength  at  60  MHz. 


TABLE  V.  PROPAGATION  PARAMETERS  AND  ATTENU AT I ON  OF  DELAY 

LINE  CANDIDATE  MATERIALS 


Acoustic 

Velocity 

meters/sec 

Acoustic  A 
for  60  MHz 

Length  for 
7ns 

Attenuation  at 

60  MHz,  7  ns 

Fused 

Quartz 

5950 

100  microns 

41  .6  mm 

.015  dB 

LiNb04 

6570 

109.5 

46 

2.45  x  10"3  dB 

PbMo04 

3660 

61.7 

25.6 

.140  dB 

TeOz 

617 

10.3 

4.32 

.450  dB 

Equation  (4-62)  is  used  to  generate  Table  VI,  which 
tabulates  the  minimum  acoustic  power  required  in  geometrical¬ 
ly  optimized  delay  lines  of  various  materials,  for  100%  Bragg 
modulation,  using  100%  fractional  bandwidth  for  Af  =  10  MHz. 
Values  are  shown  for  the  two  optical  wavelengths  of  633  nm 
(HeNe  laser)  and  830  nm  (GaAlAs  semiconductor  laser).  Also 
the  optimum  dimensions  a  and  L  are  shown  as  calculated  from 
equations  (4-53)  and  (4-50). 


TABLE  VI.  ACOUSTIC  POWER  AND  GEOMETRY  FOR  100%  BRAGG 
MODULATION,  100%  FRACTIONAL  BANDWIDTH,  Af  =  10  MHz,  AND  NO 
TIME  DELAY  REQUIREMENT,  FOR  OPTIMIZED  GEOMETRY 


X  = 

633  mm 

X 

=  830  mm 

Power 

a 

L 

Power 

a 

L 

Quartz 

24.2  mw 

.  595  mm 

816  mm 

54.6  mw 

.595  mm 

623mm 

Dense 

Flint 

Glass 

2.83 

.3100 

291 

6.38 

.310 

222 

LiNb03 

3.17 

.657 

1500 

7.15 

.657 

1144 

PbMo04 

1  .01 

.366 

506 

2.28 

.366 

386 

Slow 

Shear 

Te02 

.286 

.0617 

13.7 

.644 

.0617 

10.4 

Thus  Table  VI  shows  the  theoretically  minimum  acoustic 
power  (for  100%  modulation,  Af  =10  MHz)  with  the  delay  line 
geometry  totally  optimized.  As  can  be  seen,  this  implies 
somewhat  bizarre  geometry.  For  example,  in  the  case  of  fused 
quartz,  the  interaction  length  is  almost  1  meter,  and  the 
vertical  dimension  (perpendicular  to  both  acoustic  and  opti¬ 
cal  propagation)  is  only  1/2  millimeter. 

The  values  of  the  required  power  shown  in  this  table, 
which  vary  from  1/4  to  50  milliwatts,  are  so  low,  that  it 
becomes  apparent  there  is  little  reason  to  deal  with  these 


extreme  geometries  since  it  is  easy  to  achieve  substantially 
larger  power  levels. 


In  a  more  realistic  vein,  we  now  tabulate  the  required 
electrical  power  under  the  constraints  of  providing  7  ms 
delay,  10  MHz  bandwidth,  60  MHz  center  frequency,  and  under 
the  assumption  that  the  electrical  to  acoustic  efficiency  is 
10%,  which  is  a  fairly  good  number  for  most  cases.  Equation 
( 4-63 )  now  becomes 


Pa  >  Af  • 


1 


PE 


where  pg  *  electrical  to  acoustic  efficiency. 


Placing  the  parameters: 


(4-64) 


p  E  =  *1 
Af  =  10  MHz 

=■>>  6 

fQ  =  60  MHz 

T  =  7  ms  N 

into  the  above  equation  we  obtain: 


1/6 

70 


PE( coherent) 


2.05  x 


(4-65) 


with  all  parameters  in  the  MKS  system. 

The  required  value  of  the  vertical  dimension  a  is  given 
by  equation  (4-52),  and  the  interaction  length  L  is  given  by 
equation  (4-50).  With  the  above  parameters  these  become 

/  T  /  7  x  icr6 

a  =  v,  / -  =  v_  / - 

f o  60  x  io6 

,  meters.  , 

a  (millimeters)  =  v,  ( - - )  x  .3416  x  10“J 

sec 


(4-66) 


L  < 


X  Af  fQ 


n  v. 


12 


X  • 


10  x  60  x  10 


(4-67) 


Using  equations  (4-65),  (4-66),  and  (4-67)  we  generate 
Table  VII. 


TABLE  VII.  ELECTRICAL  POWER  AND  GEOMETRY  FOR  100%  BRAGG 
MODULATION,  1/6  FRACTIONAL  BANDWIDTH,  f=10MHz,  T=7  SEC, 
ELECTRICAL  TO  ACOUSTIC  EFFICIENCY  =0.1,  FOR  OPTIMIZED 

GEOMETRY 


A  =  633  nm 


Electrical 

Power 


The  entries  in  Table  VII  still  imply  some  rather  expen¬ 
sive  crystals.  For  example,  the  required  interaction  length 
for  quartz  is  now  136  mm,  and  for  LiNbO,  is  now  250  nm.  With 
the  exception  of  quartz,  the  required  electrical  drive  powers 
are  less  than  1  watt  for  633  nm,  and  less  than  2  watts  for 
830  nm.  It  may  be  that  to  achieve  less  expensive  crystal 
dimensions  that  the  interaction  length  L  should  be  reduced 
from  these  optimum  values,  with  an  attendant  increase  in 
drive  power. 

Table  VII,  however,  provides  an  excellent  basis  of  com¬ 
parison  of  the  candidate  materials.  Clearly,  Te02  is  the 
best  choice  for  low  drive  power.  The  extremely  small  value 
of  acoustic  wavelength  (10  microns)  and  the  small  aperture 
(4  mm)  may  cause  some  optical  problems  associated  with  the 
required  magnification.  PbMo04  runs  a  close  second  in  drive 
power,  with  more  reasonable  optical  size,  and  this  is  also  a 
strong  contender,  especialy  if  a  non-optimum  geometry  is 
considered. 

4.2  INCOHERENT  OPTICAL  OPERATION 

4.2.1  Optical  Response  of  Birefringent  Modulator 

Our  original  concept  was  to  use  the  optical  output  pat¬ 
tern  developed  by  the  image  converter  directly  in  the  second 
optical  section.  With  this  approach  the  image  converter 
output  is  imaged  on  the  second  delay  line,  and  the  light  is 
collected  by  a  lens-detector  combination.  Since  the  optical 
output  from  the  image  converter  is  radiated  from  a  phosphor 
screen,  it  is  totally  incoherent.  In  an  optically  incoherent 
system  no  spatial  frequency  plane  is  available,  and  thus  the 
phase  modulation  imparted  by  the  delay  line  cannot  be  con¬ 
verted  to  amplitude  modulation  by  manipulation  of  the  optical 
frequency  plane.  Instead,  the  birefringent  effect  in  the 
delay  line  is  utilized,  which,  with  the  correct  orientation, 
causes  a  shift  in  phase  difference  between  the  orthogonally 
polarized  components  of  the  optical  beam.  When  used  in  con¬ 
junction  with  a  polarization  analyzer,  this  causes  a  varia¬ 
tion  in  optical  amplitude. 

In  a  region  of  birefringence,  the  optical  phase  is 
shifted  by  the  active  material,  according  to  the  direction  of 
polarization  of  the  light  with  respect  to  the  crystal  axes  in 
the  direction  transverse  to  the  light  beam.  If  the  light  is 
resolved  conceptually  into  components  in  which  the  polariza¬ 
tions  are  aligned  with  the  crystal  axes,  the  two  components 
become  the  ordinary  and  extraordinary  rays  of  the  resulting 
birefringent  behavior.  Linearly  polarized  light,  then,  will 
be  equally  shared  between  these  two  orientations  if  the 


direction  of  polarization  of  the  incident  light  is  45*  to 
each  of  the  orthogonal  axes  of  the  crystal,  in  the  manner 
shown  in  Figure  4-8. 


INPUT 

POLARIZATION 


OPTICAL  PROPAGATION 
INTO  PAPER 


Figure  4-8.  Relation  of  Input  Optical  Polarization  to 

Crystal  Axes 


In  order  to  convert  the  optical  phase  shifts  induced  by 
birefringence  into  intensity  variations,  an  optical  polarizer 
is  inserted  after  the  active  crystal,  oriented  in  a  direction 
perpendicular  to  the  unmodulated  direction  of  optical  polari¬ 
zation  of  the  input  beam.  Thus,  the  output  polarizer  is 
crossed,  i.e.  it  is  perpendicular  to  any  existing  input  po¬ 
larizer.  The  input  polarizer  can  be  omitted  if  the  input 
light  is  already  polarized,  such  as  from  a  laser.  With  this 
arrangement  no  light  will  be  translated  by  the  second  polar¬ 
izer  in  the  absence  of  any  birefringence  in  the  crystal. 

When  the  material  becomes  birefringent  the  light  becomes 
elliptically  polarized,  so  that  there  is  then  a  component  in 
the  crossed  direction  that  will  pass  through  the  second  po¬ 
larizer. 

This  effect  can  be  analyzed  in  the  following  .way.  Con¬ 
sider  Figure  4-9,  which  illustrates  th/s  relations  between 
crystal  orientation  and  input  light  polarization.  The  light, 
of  amplitude  A0,  is  shown  polarized  at  45*  to  the  two  crys¬ 
tal  axes,  the  ordinary  axis  "0"  and  the  extraordinary  axis 
"E".  Conceptually,  we  may  resolve  A0  into  two  mutually 
perpendicular  components,  aligned  with  axis  E  and  0  as  shown. 
It  is  obvious  that  the  amplitude  of  each  of  these  components 
will  then  be  h0//^T~. 


Figure  4-9.  Resolution  of  Input  Polarization  into  Parallel 
and  Perpendicular  Components 


It  is  assumed  that  the  output  polarizer  is  oriented 
perpendicular  to  A0,  or  horizontally  in  the  figure.  This 
polarizer  will  pass  an  "0"  component  labeled  x  and  an  "E" 
component  labeled  y.  These  components  are  analytically  cre¬ 
ated  by  performing  a  second  resolution  of  the  E  and  0  vectors 
into  orthogonal  components,  this  time  breaking  each  of  these 
into  two  components,  one  aligned  with  the  original  A0  vec¬ 
tor  and  one  perpendicular  to  this  vector.  These  latter  per¬ 
pendicular  components  then  become  the  x  and  y  vectors.  It 
can  be  appreciated  (from  the  figure)  that  the  magnitude  of 
the  output  light  amplitude  will  be  the  magnitude  of  the  dif¬ 
ference  between  x  and  y, 

Ac  ■  |x  -  y  |  .  (4-68) 


We  now  account  for  the  birefringent  effect  as  an  optical 
phase  shift  between  the  components  aligned  in  the  E  and  0 
directions,  say  $  radians.  From  inspection  of  Figure  4-9,  we 
see  that  the  x-component  is  given  as  the  45*  component  of  the 
"0"  component,  and  the  y-component  is  given  as  the  45®  compo¬ 
nent  of  the  "E"  component.  Also  (from  the  figure)  we  see 
that  both  the  E  and  0  vectors  will  generate  45®  components 
along  the  original  A0  axis,  and  these  two  components  will 
add  with,  in  general,  the  same  phase  difference  of  4.  How¬ 
ever,  this  resultant,  whatever  its  time  phase,  is  polarized 
at  right  angles  to  the  output  polarizer  and,  hence,  will  be 
rejected  by  the  system. 


Prom  the  preceding  discussion,  we  conclude  that  the  E 
and  0  components  of  the  A0  input  can  be  written  as  the 
following  phasors,  where  we  assume  that  no  phase  shift  is 
imparted  by  the  ordinary  polarization  (reference  phase)  and 
that  all  the  differential  phase  shift  is  vested  in  the  E 
polarization. 
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where  u»0  =  optical  radian  frequency. 
Hence  x  and  y  are  evaluated  as: 
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The  desired  cross-polarized  output  component,  which  is 
that  passed  by  the  output  polarizer,  is  then. 


The  output  light  intensity  is  the  magnitude  squared  of 
the  preceding  result. 


or 

Io 

I  =  -  (1  -  cos  <J>  )  .  (4-73) 

2 


Equation  (4-72)  will  be  recognized  as  identical  to  equa¬ 
tion  (4-28).  Thus  imaging  by  Bragg  diffraction  produces  the 
same  modulation  level  as  birefringent  polarized  readout. 

It  will  be  appreciated  that  the  above  function  ap¬ 
proaches  zero  as  the  birefringent  phase  shift,  <J> ,  approaches 
zero.  For  small  values  of  $  this  behavior  is  approximated  by 
a  square-law  variation,  with  the  normal  light  intensity  equal 
to  zero  (dark  field). 

In  the  present  application  a  linear  variation  in  inten¬ 
sity  is  desired,  and  this  can  be  provided  (approximately)  by 
biasing  the  phase  shift  with  a  fixed  amount  of  birefringence 
in  the  system.  Thus,  for  example,  if  a  quarter-wave  bire¬ 
fringent  plate  is  sandwiched  with  the  active  material,  the 
value  of  <|>  in  equation  (4-73)  is  altered  by  90*,  and  (4-73) 
becomes 


o 

1  *  ~2~  (1  +  Sin  * )  '  (4-74 ) 


The  result  is  approximately  linear  with  with  the 
approximation  that  sin  $  =  <t» . 

This  result  essentially  squares  with  the  Bragg  analysis, 
equations  (4-32)  and  (4-33).  The  conclusion  is  that  with  the 
proper  birefringent  bias,  polarized  readout  of  a  sine-wave 
induced  phase  modulation  of  amplitude  <t»ra  produces  the  same 
intensity  variation  as  Bragg  imaging,  in  which  the  two  orders 
are  interfered  to  produce  the  image.  In  either  case  <f>m 
represents  the  peak  value  of  phase  variation  in  the  delay 
line  modulator.  As  seen  from  equation  (4-33)  for  Bragg  imag¬ 
ing,  if  4> m  is  tt /2  the  image  intensity  varies  from  black  to 
maximum  brightness  (unity),  which  is  also  true  of  polarized 
readout  in  equation  (4-74),  as  the  value  of  varies  from  - 
n  /2  to  +n  /2  . 


The  only  difference  between  the  two  results  is  a  slight 
non-linearity  in  the  polarized  readout.  If  4  varies  sinusoi¬ 
dally  along  the  delay  line  as  induced  by  the  acoustic  wave, 
and  we  denote  the  phase  of  the  acoustic  wave  by  6 a  as  be¬ 
fore,  (4-74)  becomes 

—  *  -  (1  +  sm  (4m  sm  ea))  (4-75) 

Io  2 


The  non-linearity  caused  by  the  outer  sine  function 
disappear?  for  small  4>m,  of  course.  At  such  small  phase 
deviations,  equations  (4-32)  and  (4-75)  become  identical. 

4.2.2  Power  Calculation  for  Incoherent  Case 

To  evaluate  the  acoustic  power  requirements  for  the 
incoherent  case,  using  polarized  readout,  the  question  of  a 
figure  of  merit  must  be  addressed.  In  Sections  4.1.6  and 
4.1.7,  the  figure  of  merit  is  based  on  the  appropriate  photo¬ 
elastic  constant,  i.e.  one  of  the  p's  of  the  photo-elastic 
tensor  described  in  4.1.2.  The  single  appropriate  p  is  used 
in  the  calculation  of  the  figure  of  merit,  since  only  the 
optical  phase  modulation  along  a  single  axis  is  of  interest 
there.  In  the  present  incoherent  case,  the  relative  phase 
between  two  perpendicular  stress-induced  optical  axes  deter¬ 
mine  the  amount  of  light  that  passes  the  crossed  polarizer. 
Thus  with  stress  along  the  x-axis  in  the  delay  line,  the 
phase  difference  is  desired  between  light  passing  through 
with  polarization  in  the  y-direction  (ordinary  ray)  and  light 
polarized  in  the  x-direction  (extraordinary  ray). 

In  other  words,  under  tensile  or  compressive  stress,  the 
delay  line  medium  acts  as  a  positive  or  negative  uniaxial 
crystal  respectively,  and  the  relative  phase  modulation  cor¬ 
responding  to  the  difference  of  the  indices  of  refraction 
nQ  and  ne  is  the  pertinent  quantity  determining  the  opti¬ 
cal  response.  Thus  instead  of  a  single  photoelastic  coeffi¬ 
cient,  the  difference  between  two  p's  is  used.  This  is  the 
quantity  Pj x  -  p12,  in  the  manner  of  equation  (4-10).  The 
desired  acousto-optical  figures  of  merit  can  be  defined  just 
as  before,  including  both  M2  and  M3 ,  except  that  pn  -  p12  is 
used  as  the  appropriate  photo-elastic  constant. 

As  is  pointed  out  in  the  discussion  following  equation 
(4-10),  the  amount  of  phase  modulation  experienced  in  the 
ordinary  wave  is  substantially  larger  than  for  the  extraordi¬ 
nary  ray.  For  fused  quartz  the  respective  stress-optical 
coefficients  are  4.2  compared  to  0.65.  Thus  one  concludes 
that  a  reasonably  accurate  approximate  result  can  be  obtained 


by  using  the  single  photo-elastic  constant  associated  with 
the  ordinary  ray  to  represent  the  relative  photo-elastic 
constarit  required  here.  This  single  photo-elastic  constant 
is  the  one  most  widely  tabulated,  and  is  the  basis  of  the 
published  values  of  Mj  and  M3  which  have  been  used  in  the 
previous  section. 

In  the  incoherent  case,  the  evaluation  of  geometrical 
constraints  is  similar  to  the  previous  calculations.  As 
before,  the  spreading  of  the  acoustic  beam  in  the  y-direction 
must  be  less  than  the  near-field  height  of  the  acoustic  beam 
leaving  the  transducer  in  this  direction.  Thus  the  con¬ 
straint  as  expressed  by  equation  (4-52)  stands  for  this 
case. 


The  determination  of  the  limiting  value  of  the  interac¬ 
tion  length  L,  however,  is  somewhat  different.  The  general 
notion  of  the  limitation  in  this  case  is  that  the  optical 
spreading  caused  by  diffraction  by  the  acoustic  wave,  as  the 
optical  propagation  proceeds  through  the  interaction  length  L 
of  the  delay  line,  must  not  become  large  enough  so  that  adja¬ 
cent  acoustic  periods  are  involved.  Such  mixing  of  the  opti¬ 
cal  ray  into  a  range  of  acoustic  phases  means  that  the  phase 
modulation  of  the  optical  ray  will  be  degraded.  The  spread¬ 
ing  between  orders,  inside  the  medium  by  the  acoustic  wave 
is: 
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(4-76) 


The  requirement  is  that  this  be  less  than  order  A,  the 
acoustic  wavelength.  Hence 


AL  , 

L  •  A0  =  -  <  A 

0  n  A 


Hence  L  must  be  less  than  order 


n  A 


(4-77) 


A  more  precise  conclusion  can  be  gained  by  considering 
the  transition  from  Raman-Nath  diffraction  to  Bragg  diffrac¬ 
tion,  as  the  interaction  length  becomes  longer.  Raman-Nath 
diffraction  is  defined  as  resulting  from  modulation  that 
effectively  occurs  in  a  single  plane.  This  means  that  the 
acoustic  wavefronts  are  aligned  with  the  light  rays  through¬ 
out  the  interaction  region.  On  the  other  hand,  as  discussed 
in  Section  4. 1.3. 2,  Bragg  diffraction  depends  on  a  continuous 
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coherent  contribution  of  energy  transfer  from  one  order  to 
another  as  the  ray  travels  through  the  interaction  region. 
Implicit  is  the  traverse  of  at  least  one  acoustic  cycle  by 
the  optical  rays.  Thus  the  definition  of  the  Bragg  regime  is 
such  that  the  optical  ray  cannot  be  constrained  to  a  particu¬ 
lar  value  of  acoustic  phase,  and  hence  the  Raman-Nath  condi¬ 
tion  is  required  in  the  incoherent  case. 


Figure  4-10  shows  the  diffraction  amplitude  as  a  func¬ 
tion  of  modulation  index  4>m  for  various  values  of  the  in¬ 
teraction  length  expressed  in  normalized  form,  L/L0.  These 
curves  were  calculated  by  W.G.  Hoefer  of  the  General  Electric 
Company  Electronics  Laboratory.  [15]  The  reference  length  is 
defined  as  follows: 


Lo  = 


no  A 


(4-78) 


Consider  the  first  order  response.  For  normalized 
length  of  value  .316  and  less  the  response  corresponds  very 
closely  to  the  Bessel  function  ( <j> jp )  associated  with 
Raman-Nath  diffraction.  For  normalized  length  3.16  and  larg¬ 
er,  the  response  is  very  close  to  sin  (4> m/2 )  associated 
with  Bragg  diffraction.  In  the  curve  shown  for  L/L0  =  1.0, 
the  response  is  well-situated  in  between  these  limits.  This 
is  particularly  clear  in  the  region  of  the  first  null  in  the 


J,  response.  Thus  the  value  L/L0  =  1.0  can  be  taken  to  be 
tne  effective  limit  of  ray  registration  with  the  acoustic 


wavefront  for  the  entire  length  of  the  interaction  region. 
From  the  definition  of  the  reference  length  we  have: 


L  <  L0  = 


(4-79) 


Hence  for  this  case,  the  limiting  value  of  a/L  is: 
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Substituting  this  result  into  equation  (4-48),  with 
♦  m  a  */2,  as  before,  for  100%  modulation: 
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This  result  may  be  rearranged  into  a  form  which  is  simi¬ 
lar  to  equation  (4-64)  as  follows: 
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Using  the  values  of  the  parameters,  N  =  70,  6  =  1/6,  and 
p  E  =  0.1,  as  before: 
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Using  the  results  expressed  by  equations  (4-52),  (4-79)  and 
(4-82)  we  can  construct  the  following  table  (Table  VIII). 


TABLE  VIII.  ELECTRICAL  POWER  AND  GEOMETRY  FOR  100%  BRAGG 
MODULATION ,  1/6  FRACTION  BANDWIDTH,  BW  =  10  MHz,  T  =  7 
MICROSEC.,  ELECTRICAL  TO  ACOUSTIC  EFFICIENCY  =0.1,  FOR 

OPTIMIZED  GEOMETRY 


X  =  514  nm 


Electrical 

Power 

a 

L 

W 

Quartz 

50  watts 

2.03  mm 

8.9  mm 

41.6  m 

Dense 

Flint 

Glass 

5.9 

1 .06 

3.2 

21  .7 

LiNb03 

6.6 

2.24 

16.3 

46 

PbMo<\ 

2.1 

1.25 

5.5 

25.6 

Te02 

(Slow 

Shear ) 

0.59 

0.21 

0.15 

4.32 

It  is  seen  that  the  power  requirements  are  substantially 
more  for  the  polarized  readout  than  were  the  case  for  Bragg 
imaging,  due  to  the  more  restrictive  condition  on  the  maximum 
permissible  interaction  length. 

In  the  experimentation,  we  had  only  two  delay  lines  that 
were  anywhere  near  the  short  interaction  length  required  for 
sufficiently  small  optical  spreading.  These  were  the  quartz 
lines  previously  described,  with  an  interaction  length  of 
15  mm,  and  an  LiNb03  line  with  an  interaction  length  of  5  mm. 
We  now  evaluate  the  100%  modulation  acoustic  power  required 
for  these  two  lines.  With  specific  geometrical  dimensions, 
without  optimization,  we  have  to  go  back  to  equation  (4-48), 
using  the  assumed  electrical  efficiency  as  an  additional 
multiplier.  In  this  way  we  calculate  the  acoustic  power  for 

1)  Quartz  line:  a  =  2  mm  L  =  15  mm 


PE  = 


=  30  watts 


2)  LiNb03  line:  a  =  5  mm  L  =  5  mm 
Pg  =  47.5  watts 


From  these  results,  it  is  apparent  that  we  were  able  to 
achieve  only  a  small  optical  modulation  index.  The  30  watts 
for  quartz  is  slightly  deceptive,  since  its  interaction 
length  is  slightly  more  than  the  allowable  value  from  Table 
VIII.  The  power  entry  in  the  table  is  the  more  appropriate 
value,  i.e.  50  watts.  The  equipment  that  was  at  hand  for  the 
experiment  was  only  capable  of  delivering  3  watts.  Thus  from 
equation  (4-48)  we  can  calculate  the  following  values  of 
with  the  corresponding  contrast  ratios  calculated  from  equa¬ 
tion  (4-34),  as  being  induced  by  3  watts  for  the  two  lines: 

1 )  Quartz  line,  at  3  watts 

<|>m  =  0.5  radians,  contrast  ratio  =  0.48 

2)  LiNb03  line,  at  3  watts 

<|>m  =  0.4  radians,  contrast  ratio  =  0.39 

The  results  imply  a  somewhat  limited  intensity  modula¬ 
tion  of  the  light  passing  through  the  second  delay  line,  i.e. 
the  delay  line  associated  with  the  space  integrating  correla¬ 
tor,  in  the  incoherent  mode.  This,  together  with  other  limi¬ 
tations  which  will  be  discussed  shortly  has  resulted  in  dis¬ 
appointing  operation  of  this  part  of  the  system. 

4.2.3  Requirements  of  Ray  Registration 

As  has  been  previously  stated,  for  proper  efficiency  in 
the  incoherent  optical  mode,  with  polarized  readout,  it  is 
necessary  that  the  light  rays  entering  the  delay  line  stay 
collimated  or  registered  with  the  acoustic  wavefronts  through 
the  entire  interaction  thickness  of  the  line.  We  have  al¬ 
ready  investigated  the  limitation  this  imposes  on  the  inter¬ 
action  length  because  of  the  diffraction  spreading  of  the 
optical  rays.  In  addition  to  this,  the  orientation  of  all 
optical  rays  entering  the  delay  line  must  be  within  strict 
angular  limits  to  satisfy  this  registration  requirement. 

Thus  if  the  light  rays  cross  more  than  1/2  acoustic  spatial 
cycle  by  having  the  wrong  angular  orientation,  or  by  diffrac¬ 
tion  spreading,  the  modulation  efficiency  is  proportionally 
degraded.  The  situation  is  pictured  in  Figure  4-4. 

Thus,  the  incoherent  configuration  demands  a  very  large 
loss  of  light  in  the  system.  Of  the  approximate  full  hemi¬ 
sphere  of  angular  radiation  from  the  output  surface  of  the 
image  converter,  only  that  narrow  range  of  angle  can  be  used 
which  properly  aligns  with  the  acoustic  wavefronts.  This 
angle  is  approximately  given  by  the  ratio  of  1/2  acoustic 
wavelength  to  the  interaction  thickness.  In  the  case  of  the 
Andersen  line,  this  is  1/150  radians.  The  selection  of  this 
range  of  angle  is  accomplished  by  a  stop  in  the  second  opti¬ 
cal  correlator  frequency  plane.  This  narrow  selection  has  to 


occur  in  only  one  dimension,  and  the  full  speed  of  the  lenses 
can  be  used  in  the  orthogonal  direction,  which  in  our  experi¬ 
ment  is  1/4  radian.  Thus  the  system  is  equivalent  to  a  cir¬ 
cular  aperture  system  in  which  the  aperture  diameter  is  the 
geometrical  mean  of  these  two  numbers,  or  is  equivalent  to 
approximately  a  f/24  optical  system,  which  is  very  slow. 

With  the  maximum  intensity  level  available  from  the  image 
converter  (ITT  Model  4108)  a  theoretical  calculation  shows 
that  with  the  parameters  that  have  been  suggested,  a  positive 
signal  to  quantum  noise  ratio  should  be  available  at  the 
detector.  With  marginal  theoretical  performance,  and  addi¬ 
tional  degradation  encountered  in  the  practical  realization, 
satisfactory  performance  of  the  incoherent  section  has  not 
been  achieved.  Part  of  the  reason  for  the  poor  performance 
of  this  section  is  believed  to  be  curvature  of  the  acoustic 
wavefronts,  so  that  optical  registration  cannot  occur.  It 
must  be  remembered  that  at  the  60  MHz  carrier  selected,  the 
interaction  length  is  120  times  the  half  wavelength  of  the 
acoustic  propagation. 

4.3  LIGHT  BUDGET  ANALYSIS 

We  now  consider  the  light  budget  analysis  of  the  system, 
and  the  compatibility  of  the  light  levels  generated  by  each 
of  the  components  of  the  system  with  the  light  requirements 
of  the  next  stage. 

4.3.1  Time  Integrating  Correlator 

The  illumination  of  the  system  was  provided  by  the 
GaAlAs  laser  diode.  The  radiation  obtained  from  this  device 
was  used  without  any  spatial  filtering  with  a  pinhole.  Use 
of  such  a  pinhole  matched  in  size  to  the  diffraction-related 
divergence  of  the  beam  would  produce  a  uniform  beam  without 
significant  light  power  loss. 

The  effective  beam  spread  of  this  device  is  stated  by 
the  manufacturer  to  be  10°  *  25°.  This  corresponds  to  0.076 
steradians.  The  output  wavelength  is  830  nm.  The  maximum 
rated  output  power  of  the  diode  is  10  mW.  As  can  be  seen 
from  Figure  3-2  the  output  power  for  a  given  drive  current  is 
strongly  temperature  dependent.  Nevertheless,  the  output 
should  be  kept  to  10  mW,  regardless  of  the  temperature.  To 
exercise  a  comfortable  safety  margin,  the  device  was  biased 
in  actual  operation  to  produce  an  average  of  5  mW.  This 
corresponds  to  a  solid  angular  power  density  of 

5  mW  /  .076  steradians  =  .066  watts/steradian. 


Referring  to  Figure  2-1  we  see  that  the  light  collected 
by  the  first  collimator  in  the  system  is  projected  as  a  par¬ 
allel  beam  to  the  first  delay  line.  This  plane,  in  turn,  is 
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imaged  to  the  image  intensifier.  The  light  power  density  at 
the  delay  line  is  found  from  the  above  value  divided  by  the 
square  of  the  focal  length  of  the  first  collimator,  which  is 
approximately  15  cm.  Thus: 

IDL1  *  . 066/1 52  »  293  *  10"6  watts/cm2  (4-83) 

Since,  in  a  Bragg  modulator,  all  the  light  is  essential¬ 
ly  delivered  to  the  zero  and  first  order,  and  these  orders 
are  both  used  to  produce  an  image,  the  average  light  level  at 
the  image  intensifier  will  be  equal  to  the  above  value.  We 
now  calculate  the  response  of  an  intensifier  to  a  given  input 
light  level.  The  first  issue  is  the  conversion  efficiency  of 
the  light  to  emitted  photocathode  current,  which  is  released 
to  the  accelerating  field  in  the  device  by  the  photoemitting 
film. 


A  standard  material  which  is  used  for  this  purpose  is 
the  "S-20"  photocathode.  Since  we  are  using  infra-red  light 
at  830  nm  to  illuminate  the  intensifier,  the  device  was  pro¬ 
vided  with  an  "extended  red"  photocathode,  which  is  essen¬ 
tially  an  S-25  film.  The  spectral  response  of  these  two 
materials  are  shown  in  Figure  4-11.  These  curves  were  taken 


Figure  4-11.  Spectral  Response  Curves  for  S-20  and  S-25 
(extended  red)  Photocathodes.  Data  taken  from  Reference  [16] 
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from  data  published  by  ITT  [16].  In  this  figure,  the  sensi¬ 
tivity,  in  emitted  amperes  per  incident  optical  watt,  is 
plotted  vs.  the  optical  wavelength  in  nm.  The  locus  of  vari¬ 
ous  values  of  quantum  efficiency  in  units  of  electrons/photon 
are  also  plotted  on  the  same  graph.  It  is  seen  that  while 
the  quantum  efficiency  is  in  excess  of  10%  in  the  wavelength 
region  of  400  nm,  at  830  nm,  the  S-20  photocathode  is  nearly 
inactive,  and  the  S-25  photocathoda  exhibits  a  QE  of  approxi¬ 
mately  1%.  The  following  table  summarizes  the  QE  and  sensi¬ 
tivity  in  A/to  for  the  important  laser  wavelengths. 


TABLE  IX. 


QE 

S(A/W) 

S-25  at 

515  nm 

00 

o 

• 

.035 

633  nm 

.05 

.027 

830  nm 

.01 

.0067 

The  sensitivity  is  given  in  terms  of  the  QE  =  n  by  the 
following  formula,  which  can  be  regarded  as  the  definition  of 
QE: 

n  q  n  q  X 

S  =  -  =  -  amperes/watt  (4-84) 

h  u  he 

Thus,  for  example,  at  830  nm,  we  have: 

.01  x  1 .6  x  10“19  x  .830  x  10"6 

S  ■  -  =  .0067  amps/watt 

6.6  x  10-34  x  3  x  108 

agreeing  with  the  entry  in  Table  IX. 

Having  this  value,  we  may  now  find  the  maximum  allowable 
average  light  power  density  incident  on  the  photocathode  of 
the  intensifier,  in  terms  of  the  allowable  photocathode  cur¬ 
rent  density.  The  latter  is  the  best  quantity  to  properly 
define  the  onset  of  destructive  levels  of  operation  for  the 
device.  For  the  4108  intensifer,  ITT  recommends  operating 
the  photocathode  current  density  to  a  value  equal  or  less 
than  0.1  microamperes/cm2.  Since  the  sensitivity,  S,  as 
defined  before  relates  the  the  incident  light  power  density 
I image  to  fche  current  density  icath*  we  have 


13  H-ft  O 


information  presented  in  this  figure  is  basically  the  emis¬ 
sion  efficiency  of  the  phosphor  screen  of  the  intensifier 
expressed  in  optical  watts  radiated  per  watt  impinging  on  the 
screen  in  the  form  of  electron  flow.  This  is  shown  as  a 
function  of  optical  wavelength.  The  actual  number  plotted  on 
the  chart  is  watts/nm-watts.  This  data  was  taken  from  data 
published  by  ITT  [16]. 

Two  phosphors  are  shown  on  the  chart,  P-1  and  P-20.  The 
P-1  phosphor  is  the  material  used  on  our  larger  40  mm  device, 
and  is  the  choice  that  provides  a  good  integration  time  ( > 1 0 
milliseconds).  The  P-20  is  used  on  the  smaller  device 
(25  mm)  and  this  device  is  the  one  which  has  been  used  to 
report  the  present  results. 

As  seen  from  the  chart,  at  the  most  favorable  wave¬ 
lengths,  electron  energy  is  being  converted  to  light  at  bet¬ 
ter  than  50%  efficiency  (10“3  *  500  nm).  Over  the  active 
optical  band,  the  average  efficiency  can  be  estimated  to  be 
better  than  10-4  photons/nm-ev,  or  .05  photons/ev.  In  fact, 
the  overall  efficiency  published  by  ITT  is  .063  photons/ev. 
Remembering  that  1  coulomb-volt/sec.  *  1  watt,  we  can  calcu¬ 
late  the  output  optical  power  density  in  the  following  way: 


Output  Optical  coulombs 

=  *cath  ( 


Power  Density 


sec-cm 


1 


coulombs 

electron 


photons 


vbia«s  (volts)  X  Kn  — - — 

Dlas  P  electron-volt 


he  watts/sec 

\  photon 


(4-86) 


where  icath  =  photocathode  current  density 

(amps/cm2 ) 


q  =  electronic  charge  in  coulombs/electron 

Vbias  =  accelerating  potential  across  the 
cathode-phosphor  screen  gap  in  the 
intensifier. 


Kp  =  conversion  efficiency  of  the  phosphor 
in  photons/ev 

h  =  Planck’s  constant  =  6.6  x  10~34  joule 
sec. 


c 


velocity  of  light  =  3  x  io8  meters/sec 


Thus  for  the  P-20  photocathode: 


he 

po  =  ^-cath  vbias  Kp  x 

_7  6.6  x  10"34  x  3  x  io8 

-  10  7  x  6000  x  .063  x  - 

1 .6  x  10"19  x  .55  x  1 0-6 

=  85  x  io"6  watts/cm^ 


(4-87) 


Comparing  with  4-85,  the  optical  power  gain  is  seen  to 
be  5.7. 


This  result  may  be  checked  by  using  the  published  photo¬ 
metric  data  for  the  device.  The  rated  optical  gain  of  the 
device  at  the  favorable  input  wavelength  of  550  nm  is  50  ft- 
lamberts/ft-candle.  This  means  that  the  total  output  light, 
if  integrated  over  the  entire  solid  angle  of  a  hemisphere,  is 
50  lumens/ft2  for  an  illumination  of  1  lumens/ft2 .  The  rated 
photocathode  sensitivity  in  microamperes/lumen  is  225.  Thus 
assuming  a  maximum  photocathode  current  as  before  of  0.1  x 
1 0”6  amperes,  we  have: 


1  image 


10~7 

225  x  io-6 


444  x  10”®  lumens/cm2  incident 


•T' 

’.‘/■•.V/y 


>.r~  - 
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Therefore , 

*out  =  1 image  x  50 


.022  lumens/cm^ 


Over  the  effective  optical  wavelength  band  of  the  P-20 
phosphor,  the  conversion  between  lumens  and  watts  is  500. 
Thus  the  optical  output  power  density  can  be  calculated  as: 

5  1 

Iout  =  *022  lumens/cnr  x  — - -  watts/lumen 
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44  x  io“®  watts/cm^ 


(4-88) 


This  amount  of  power  corresponds  to  the  photometric 
level  of  20.6  ft-lamberts. 

This  result  is  reasonably  close  to  that  result  expressed 
by  equation  (4-87).  Conservatively,  we  will  use  44  micro¬ 
watts/cm2  as  the  allowable  output  level  of  the  device. 

4.3.2  Space  Integrating  Correlator,  Incoherent  Mode 

With  the  intensifier  optical  output  level  calculated  in 
the  previous  section,  we  now  evaluate  the  amount  of  light 
that  reaches  the  final  optical  detector,  and  the  resulting 
SNR,  based  on  quantum  noise,  or  the  NEP  (noise  equivalent 
power)  of  the  final  detector. 

The  radiation  from  the  phosphor  screen  can  be  modelled 
as  a  Lambertian  surface,  which  means  that  the  radiation  in 
any  direction  is  proportional  to  the  projected  surface  area 
causing  the  radiation.  Thus  the  optical  flux  is  proportional 
to  cos  0 ,  where  0  is  the  angle  from  normal  to  the  surface. 

The  subtense  of  a  sphere  in  solid  angle  is  4*  steradi- 
ans,  and  into  a  hemisphere  is  2n  steradians.  If  the  2*  angle 
is  integrated  with  a  weighting  function  cos  0 ,  the  result  is 
* .  Thus  the  radiation  normal  to  the  surface  is  as  strong  per 
steradian,  as  if  all  the  radiation  were  being  directed  uni¬ 
formly  into  a  solid  angle  of  *  steradians.  Hence  the  power 
luminance,  B  in  watts/steradian-cm-"2 ,  is  given  by 

rout  y 

B  =  Power  Luminance  =  -  watts/steradian-cm^ 

IT 


(4-89) 


where  IOU£  is  the  power  density  being  radiated  by  the 
surface,  i.e.  the  output  screen  of  the  intensifier. 

Due  to  the  requirements  of  waveform  registration  which 
has  been  explained  in  Section  4.2.3,  only  a  small  part  of 
this  radiation  can  be  accepted  in  angle.  The  desired  angle 
is  stopped  in  the  frequency  plane  of  the  space  integrating 
correlator,  and  this  amounts  to  an  effective  optical  speed  of 
f/22 . 


The  optical  f/  -  number  is  based  on  a  circular  aperture, 
and  is  equal  to  the  ratio  of  the  focal  length  to  the  diame¬ 
ter.  The  solid  collection  angle  is  thus  equal  to  the  area  of 
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such  an  aperture  divided  by  the  square  of  the  focal  length: 

Solid  Collection  Angle  (ir/4)D2  n 

“  f2  Q  =  ’  “  _______ 

In  Space-Int.  Correlator  °  p2  4(f/)2 

(4-90) 

where  f/  =  optical  f/  -  number  of  apertured  lens. 

The  power  density  in  the  image  reaching  the  second  delay 
line  DL2  is  then  the  product  of  the  solid  collection  angle 
and  the  luminance  B: 

*out 

*DL2  =  B  °S  - - T  (4"91) 

4 ( f/) 2 

We  are  now  interested  in  the  amount  of  light  that  passes 
through  the  delay  line  corresponding  to  a  spatial  degree  of 
freedom,  as  determined  by  the  bandwidth  of  the  acoustic  sig¬ 
nal  propagating  in  the  delay  line.  As  denoted  in  the  earlier 
discussion  on  delay  line  modulation,  the  spatial  element  is 
b,  which  if  the  bandwidth  is  10  MHz  and  the  carrier  frequency 
Ts  60  MHz,  is  6  acoustic  wavelengths.  Thus  the  area  of  the 
spatial  element  is  ab. 

Each  time  a  period  passes  equal  to  the  reciprocal  of  the 
signal  bandwidth,  another  independent  time  sample  has  been 
achieved.  Thus  the  single  degree  of  freedom  in  both  space 
and  time  is  associated  with  the  spatial  area  ab,  and  the  time 
increment  1/BW.  We  may  now  calculate  the  numiEer  of  elec¬ 
trons,  corresponding  tc  quantum  events  from  the  final  detec¬ 
tor  as  the  image  power  density  at  DL2  times  the  space  and 
time  elements,  divided  by  the  energy  per  photon,  or  quantum, 
times  the  quantum  efficiency  of  the  final  detector.  Thus  we 
have: 


No.  Circuit 

Electrons  TDL2  x  ab  *  nD  :out  ab  X  nD 

Degree  of 
Freedom 


h  u  (BW) 


4 ( f /  )2  h  c  (BW) 


The  notation  follows  that  given  in  equation  (4-86). 


For  the  dimensions  of  the  space  integrating  correlator 
being  described  we  have: 


44* 10~6* .2* . 06* (.55*  10~6)x  .1 
6  4x222*6.6*  10_3^*3x108x  10*  106 


=  8  electrons/degree  of  freedom 
IQut  watts/cm2 

a  =  .2  cm 
b  =  .06  cm 

The  power  SNR  based  on  quantum-limited  noise  is  simply 
the  number  of  events  per  degree  of  freedom.  Thus  the  result 
of  8  electrons  per  degree  of  freedom  corresponds  to  an  SNR  of 
9  dB.  This  result  assumes  a  detector  which  is  quantum  noise 
limited  such  as  a  photomultiplier,  and  the  assumed  QE  of  0.1 
is  typical  for  this  device. 

The  foregoing  analysis  is  based  on  a  detector  like  a 
photomultiplier  which  can  sense  groups  of  quanta,  according 
to  its  QE.  If  a  conventional  solid  state  detector  is  used 
the  available  optical  power  must  be  compared  to  the  NEP 
(noise  equivalent  power)  of  the  device.  From  the  preceding 
analysis  the  optical  power  into  the  final  detector  per  spa¬ 
tial  degree  of  freedom  is  clearly  given  by: 


Power/spatial  degree  of  freedom  =  PQ  =  IDL2  x  ab 


out 


ab 


4 ( f/  )2 


44  x  10~6  x  .2  x  ,Q6 
4  x  222 


2.7  x  10-1®  watts 


(4-93) 


The  rated  NEP  of  a  typical  avalanche  detector-amplifier 
combination,  such  as  the  RCA  Type  C30818  is  3.6* 10”14  watts/ 
/HzT 


i 


Hence,  for  a  bandwidth  of  107  Hz,  we  have: 


NEP 


3.6  x 


1.14  x  10_1°  watts 


Thus  for  PQ  =  2.7  x  10“10  watts,  the  SNR  is 
2 . 7 

SNR  =  — -  =  2.4  =*>  4  dB 

1.14 


(4-94) 


In  an  effort  to  pick  up  a  signal  in  the  incoherent  con¬ 
figuration,  a  spectrum  analyzer  was  used  on  the  detector 
output.  The  frequency  resolution  of  this  analyzer  was  ap¬ 
proximately  103  Hz.  Under  these  conditions  the  SNR  on  the 
spectrum  analyzer  display  should  have  been 

Po  2.7  x  10“ 10 

SNR  =  - - -  =  - - -  =  237  =>  24  dB 

NEP  /BW  3.6  x  10-14  / 103 


(4-95) 


Even  with  this  narrow  band,  a  signal  was  not  observable. 
This  result  helped  convince  us  that  the  curved  acoustic  wave- 
front  problem  described  earlier  is  the  most  likely  explana¬ 
tion  of  a  lack  of  signal. 

As  equation  (4-95)  indicates,  the  SNR  of  a  Johnson 
noise-limited  detector-amplifier  varies  as  1//BW,  not  1/BW  as 
in  a  quantum-limited  detector.  However,  for  a  high  enough  BW 
the  device  eventually  becomes  QN  limited,  when  the  quantum 
noise  exceeds  the  Johnson  noise.  This  happens  when: 

_ _  h  c  (BW)  ,  X  nD  <NEp)  2 

NEP  /BW  =  -  =$>  BW  =  f - ■) 

\  ri  p  he 


Assuming  ng  *  0.5  for  the  RCA  detector,  and  X  =  550 
nm,  we  have  a  cross-over  bandwidth  of  2.5  GHz. 

4.3.3  Coherent  Light  Modulator  Option 


As  mentioned  earlier,  it  is  believed  that  the  difficul¬ 
ties  with  the  coherent  optical  mode  in  the  space  integrating 
correlator  can  be  resolved  with  the  use  of  coherent  light. 
This  will  require  that  the  incoherent  light  developed  by  the 
image  intensifier  is  converted  to  coherent  light,  or  that  the 
coherent  light  delivered  to  the  intensifier  by  the  time-inte¬ 
grating  correlator  be  integrated  and  transformed  to  a  coher¬ 
ent  output. 
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In  this  section  we  investigate  the  use  of  a  coherent 
modulator  to  perform  the  functions  of  integration  and  the 
relaying  of  the  spatial  information  in  optically  coherent 
form  to  the  space-integating  correlator.  Two  devices  which 
are  candidates  for  this  purpose  are  the  Liquid  Crystal  Light 
Valve  (LCLV)  made  by  Hughes  Aircraft,  and  the  Itek  Pockel's 
Readout  Optical  Modulator  (PROM).  The  difficulty  with  using 
these  devices  directly,  without  an  image  intensifier  in  the 
system,  is  that  neither  the  LCLV  nor  the  PROM  is  sensitive  to 
the  infrared  wavelength  generated  by  the  input  laser 
(830  nm).  With  this  input  wavelength,  the  intensifier  is 
still  needed  to  convert  the  830  nm  pattern  to  light  in  the 
550  nm,  region,  which  can  then  energize  either  the  LCLV  or 
the  PROM. 

It  is  believed  that  the  LCLV  is  the  most  suitable  of  the 
two  devices  for  this  application,  being  a  simpler  device,  and 
having  a  moving  time  integration  window,  rather  than  the 
"store  and  dump"  mode  used  by  the  PROM.  Hence,  we  will  here 
analyze  the  compatibility  of  the  LCLV  with  the  rest  of  the 
proposed  system. 

Figure  4-13  shows  a  cross  section  of  the  LCLV.  The 
device  is  in  two  sections  which  are  optically  isolated  from 


Figure  4-13.  Cross-sectional  Schematic  of  the  Liquid  Crystal 
Light  Valve  (Taken  from  Reference  17) 


each  other.  The  writing  light  comes  in  from  the  right  and 
falls  on  a  CdS  photoconductive  layer.  This  layer  interacts 
with  a  bias  voltage  (a.c.)  to  produce  an  electric  field 
across  the  liquid  crystal  in  the  second  section,  which  is  an 
image  of  the  written  pattern. 

This  field  causes  a  change  in  the  polarization  of  the 
coherent  light  which  is  reflected  into  and  out  of  the  liquid 
crystal  layer.  This  effect  is  depicted  in  Figure  4-14.  It 
is  called  by  Hughes  a  "hybrid  field  effect  mode"  which  in¬ 
volves  both  the  effect  of  optical  activity  due  to  the  twisted 
nematic  effect  in  the  liquid  crystal  and  a  birefringent  ef¬ 
fect.  A  detailed  description  of  this  process  is  contained  in 
Reference  17. 


Figure  4-14.  Operation  of  the  Hybrid  Field  Effect  Mode  Light 
Valve:  (a)  the  off-state;  (b)  the  on-state 

(Taken  from  Reference  17) 


Figure  4-15  shows  the  spatial  frequency  response  of  the 
LCLV  (MTF  curve)  as  measured  by  the  first  order  diffraction 
in  a  coherent  readout.  A  grating  of  variable  spatial  fre¬ 
quency  was  written  into  the  device,  and  the  value  of  this 
spatial  frequency  is  the  abscissa  of  this  graph.  The  inten¬ 
sity  of  the  diffracted  order  is  down  by  a  factor  of  1/2  (3 
dB)  at  approximately  18  cycles  (line  pairs)  per  mm.  Thus  the 
performance  of  the  device  at  10  cycles/mm,  which  is  the  grat¬ 
ing  frequency  generated  by  the  image  of  the  delay  line, 
should  be  excellent. 


Figure  4-15.  Relative  Intensity  of  the  First  Order 
Diffraction  Spot  as  a  Function  of  Spatial  Frequency. 
(Intensity  is  normalized  to  one  at  zero  spatial  frequency. 
One  square  inch  of  LCLV  was  used  in  the  measurement.  Second 
order  diffraction  spot  intensity  was  <1%  of  first  order  at 
all  frequencies.  Readout  wavelength  =  514.5  nm. 
(Taken  from  Reference  17) 


Finally,  the  sensitometry  of  the  LCLV  is  shown  in  Figure 
4-16.  Both  scales  in  this  presentation  are  log-intensity  in 
the  manner  of  a  sensitometry  curve  for  photographic  film. 

The  curve  exhibits  a  threshold  knee  somewhat  below  1  micro¬ 
watt/cm2  and  a  saturation  shoulder  in  the  region  of  100  mi¬ 
crowatts/cm2  .  For  the  best  and  most  linear  operation,  the 
normal  light  level  without  signal  modulation  should  be  in  the 
middle  of  the  linear  range  of  this  curve,  which  is  approxi¬ 
mately  10  microwatts/cm2 .  Thus  the  44  microwatts/cnr  which 
was  found  to  be  available  from  the  intensifier  (and  which  may 
be  substantially  higher)  in  the  previous  section,  should  be 
quite  adequate  to  drive  the  LCLV.  This  should  be  true  wheth¬ 
er  the  direct  coupling  through  the  fiber  optic  faceplates  of 
the  two  devices  is  used,  or  the  output  screen  of  the  intensi¬ 
fier  is  optically  imaged  with  a  reasonably  fast  lens  to  the 
input  screen  of  the  LCLV.  To  image  with  a  loss  of  1/4  would 
require  an  f/1  lens.  This  is  fast,  and  with  high  resolution 
would  be  a  relatively  expensive  lens,  but  it  is  quite  practi¬ 
cal  for  incoherent  imaging. 


Figure  4-16.  Sensitometry  Curve  for  LCLV  over  a  3  mm 
Diameter  Region  of  Aperture 
(Taken  from  Reference  17,  by  Casasent) 


5  RECOMMENDED  AOAP  DESIGN 


5.1  COHERENT  AND  NON-COHERENT  OPTIONS 

In  the  previous  sections  we  have  considered  both  opti¬ 
cally  incoherent  and  coherent  operation  of  the  space-inte¬ 
grating  correlator  of  the  AOAP. 

The  incoherent  mode  is  a  rather  natural  conceptual  pos¬ 
sibility,  since  the  input  optical  pattern  to  this  section  of 
the  AOAP  is  developed  on  the  phosphor  screen  of  the  image 
intensifier,  which  radiates  totally  incoherent  light.  The 
original  hope  was  that  an  optically  incoherent  correlator 
could  accept  this  input  directly  and  thus  transmit  a  carrier 
borne  correlation  function  to  the  output  detector  as  an  in¬ 
tensity  modulation. 

It  was  found  that  with  the  maximum  intensity  that  can 
safely  be  derived  from  the  output  screen  of  the  intensifier, 
there  is  a  fundamental  light  loss  which  results  in  an  output 
SNR  of  9  dB,  for  a  bandwidth  of  10  MHz,  and  a  photomultiplier 
detector  with  a  QE  of  10%.  Alternatively  an  SNR  of  4  dB  is 
achieved  with  a  sensitive  avalanche  silicon  detector. 

These  SNR's  come  about  because  of  the  narrow  restriction 
on  light  rays  that  can  be  used  in  the  acoustic  delay  line 
modulator,  and  the  requirement  of  ray  registration  with  the 
acoustic  wavefronts.  The  indicated  SNR's  were  not  achieved 
in  the  experiment.  Narrow  band  tests  were  also  performed  on 
the  system  to  enhance  the  SNR  in  the  incoherent  mode,  also 
without  success.  Even  if  the  incoherent  mode  can  be  made  to 
work  according  to  the  theoretical  SNR,  it  appears  the  perfor¬ 
mance  will  be  marginal. 

On  the  other  hand,  it  has  been  found  to  be  easy  to 
achieve  good  performance  with  a  coherent  optical  system.  To 
show  this,  an  optical  grating  of  the  same  spatial  frequency 
as  the  correlation  pattern  achieved  on  the  image  converter 
was  placed  in  the  position  of  the  image  converter  and  illumi¬ 
nated  with  an  expanded  HeNe  laser  beam.  Now  little  or  no 
light  has  to  be  sacrificed  to  match  the  required  angular 
registration  with  the  acoustic  wavefronts.  A  Bragg  optical 
delay  line  can  be  used,  as  was  done  in  the  first  optical 
section,  and  with  the  generation  of  only  two  optical  orders, 
the  full  potential  of  the  Bragg  thickness  can  be  used  to 
generate  an  amplitude-modulated  signal.  A  slight  curvature 
of  the  acoustic  wavefronts  does  not  degrade  the  amount  of 
diffraction,  but  merely  causes  a  slight  curvature  of  the 
diffracted  optical  wavefronts.  The  concept  is  essentially 
that  of  a  volume  hologram.  The  signal  in  both  the  incoherent 


and  coherent  cases  is  generated  by  the  action  of  a  moving 
grating  (the  acoustic  wave)  passing  by  the  image  of  the  cor¬ 
relation  pattern  (a  fixed  grating).  In  the  coherent  case  the 
signal  carrier  can  be  alternatively  thought  of  as  being  gen¬ 
erated  by  a  heterodyne  action.  This  is  the  interference 
between  an  optical  order  not  diffracted  by  the  delay  line 
becoming  coincident  with  an  order  which  is  diffracted  by  both 
the  input  grating  and  the  delay  line.  The  light  which  is 
diffracted  by  the  delay  line  has  a  doppler  shift  equal  to  the 
acoustic  frequency,  and  the  acoustic  carrier  recovered  in  the 
detector  is  the  beat  frequency  between  the  two  orders.  The 
device  which  we  have  selected  to  provide  spatial  modulation 
of  a  coherent  optical  beam,  thus  serving  the  role  of  the 
optical  grating  in  the  experiment  which  has  been  described, 
is  the  liquid  crystal  light  valve. 

An  important  restriction  on  the  use  of  an  LCLV  in  the 
AOAP  system  is  the  spectral  range  of  the  photoconductor  used 
as  an  input  receptor  for  the  device.  In  the  case  of  the 
Hughes  LCLV  this  is  a  film  of  cadmium  sulfide,  which  does  not 
respond  well  at  optical  wavelengths  much  longer  than  550  nm. 
Thus  if  an  infrared  illumination  is  used,  as  has  been  done 
with  the  present  configuration,  in  which  a  GaAlAs  laser  diode 
is  the  light  source  operating  at  830  nm,  there  is  no  possi¬ 
bility  that  the  LCLV  can  respond  to  this. 

However,  we  have  already  noted  that  the  required  level 
of  input  light  for  this  device  is  well  matched  to  the  expect¬ 
ed  output  level  of  the  Image  intensifier  which  has  already 
been  used  in  the  experiments  performed.  An  output  level  of 
at  least  44  microwatts/cm2  is  obtainable  from  the  image  in¬ 
tensifier,  and  an  input  level  of  10  microwatts/cm2  puts  one 
in  the  middle  of  the  linear  range  of  operation  of  an  LCLV. 

Thus  since  the  input  level  of  the  LCLV  and  the  output 
level  of  the  image  intensifier  are  compatible,  the  retention 
of  the  intensifier  in  the  optical  design  offers  a  solution  to 
the  wavelength  compatibility  problem.  The  intensifier,  as 
explained  in  Section  4.3.1,  uses  an  S-25  photoemitter,  which 
is  an  "extended  red"  material,  and,  as  indicated  in  Figure  4- 
11,  responds  quite  satisfactorily  to  input  radiation  at 
830  nm.  Also  as  previously  noted,  an  abundance  of  optical 
power  at  830  nm  is  available  from  the  time  integrating  corre¬ 
lator,  being  a  coherent  system,  with  a  high  degree  of  light 
utilization  efficiency  from  input  to  output. 

It  is  estimated  that  the  output  of  the  image  intensifier 
can  be  relayed  to  the  input  plane  of  the  LCLV  by  either  a 
fast  high  quality  1:1  lens,  or  by  intimate  contact  between 
the  fiber  optical  faceplates  provided  on  both  devices.  In 


either  case,  the  coupling  should  be  possible  with  a  loss  of 
not  more  than  4:1,  a  result  which  is  compatible  with  the 
optical  power  levels  cited  above. 

It  has  been  shown  that  both  the  spatial  frequency  re¬ 
sponse  of  the  image  intensifier  and  the  LCLV  are  well  within 
the  carrier  spatial  frequency  developed  by  the  AO  delay  lines 
of  20  cycles/mm. 

Finally,  the  spectral  response  of  the  LCLV  and  the  out¬ 
put  optical  spectrum  of  the  image  intensifier  phosphor  are 
well  matched.  In  Figure  5-1  we  show  a  spectral  response 
curve  of  the  cadmium  sulfide  as  measured  by  its  effect  on  the 
LCLV.  It  is  seen  that  the  primary  response  band  is  between 
500  and  600  nm.  Referring  to  Figure  4-12  it  can  be  observed 
that  the  emission  of  both  the  P-1  and  the  P-20  phosphors  are 
well  matched  to  the  LCLV  photoconductor.  In  particular,  the 
P-1  phosphor  straddles  almost  precisely  the  band  of  the  LCLV. 


Figure  5-1 .  Light  Valve  Photo  Response  Characteristics 
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In  view  of  these  considerations,  our  conclusion  is  that 
a  coherent  optical  mode  for  the  space- integrating  correlator 
will  produce  far  superior  performance  of  the  AOAP  system,  and 
this  is  the  configuration  that  is  planned  for  the  next  exper¬ 
iment  . 

5.2  PROPOSED  REDESIGN  OF  THE  AOAP  OPTICAL  CHANNEL 

In  the  next  phase  of  this  effort,  we  will  implement  the 
system  shown  in  Figure  5-2.  This  design  is  very  similar  to 
the  experiment  which  has  been  accomplished,  with  the  impor¬ 
tant  addition  of  a  liquid  crystal  modulator  to  convert  the 
space  integrating  correlator  to  an  optically  coherent  mode. 

A  spatial  filter  has  been  added  to  the  GaAlAs  laser  to 
produce  a  "clean"  spatial  illumination  for  the  first  delay 
line.  Imaging  from  the  first  delay  line  to  the  image  inten- 
sifier  will  be  accomplished  by  a  telecentric  system  as  be¬ 
fore.  Wollensak  f/4.5,  241  mm  lenses  obtained  for  this  pur¬ 
pose  will  be  used.  Again,  no  spatial  optical  stop  will  be 
required  in  the  frequency  plane,  since  the  delay  line  oper¬ 
ates  in  the  Bragg  mode,  and  the  resulting  two  diffraction 
orders  provide  the  required  spatial  carrier  based  on  light 
amplitude  variation.  In  the  first  phase  of  the  new  experi¬ 
ment,  an  Isomet  glass  Bragg  cell  will  be  used  as  before  for 
an  initial  demonstration,  since  a  high  degree  of  modulation 
can  be  obtained  with  these  devices.  They  do  not  provide  the 
7  microseconds  time  window  that  is  the  goal  of  this  effort, 
but  that  will  be  achieved  later  with  the  Andersen  lines  pre¬ 
viously  obtained,  or  with  new  lines,  possibly  constructed  of 
PbMo04 .  In  any  case,  a  spatial  carrier  frequency  of  10  cy¬ 
cles/mm  will  be  used.  For  the  glass  or  the  PbMo04  lines  this 
required  a  temporal  carrier  frequency  for  the  acoustic  signal 
of  approximately  35  MHz.  The  delay  line  will  be  mounted  on  a 
transverse  slide  with  a  differential  micrometer  control  for 
careful  phase  adjustment. 

For  the  initial  experiment  the  25mm  diameter  image  con¬ 
verter  will  be  used,  which  has  the  P-20  phosphor,  and  a 
fiber-optic  faceplate  for  coupling  from  the  phosphor  screen 
to  the  next  device. 

Following  the  image  converter  will  be  the  LCLV,  to  be 
obtained  from  Hughes  Aircraft.  Although  the  P-20  phosphor  is 
not  quite  as  good  a  spectral  match  to  the  LCLV  as  the  P-1,  it 
is  a  reasonably  good  match.  The  P-1  can  be  used  in  later 
experiments.  The  P-20  does  not  provide  any  significant  time 
integration  (on  the  order  of  microseconds),  but  the  response 
characteristics  of  the  LCLV  will  provide  a  very  satisfactory 
integration  time. 
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The  readout  liquid  crystal  plane  is  located  on  the  other 
side  of  the  LCLV  device  from  the  input  photoconductor.  A 
collimated  laser  beam  (expected  to  be  a  HeNe  laser  providing 
633  nm)  will  be  introduced  into  this  layer  with  a  beam  split¬ 
ter  oriented  at  45*  to  the  optical  axis,  as  shown  in  Figure 
5-2.  This  will  inject  the  beam  parallel  to  the  optical  axis 
into  the  liquid  crystal  (i.e.  perpendicular  to  the  liquid 
film)  and  the  light  will  be  reflected  and  will  pass  out  of 
the  liquid  film.  The  passage  through  the  liquid  film,  and  a 
polarizer,  will  provide  the  desired  amplitude  modulation  in 
accordance  with  the  pattern  relayed  from  the  image  converter. 
The  beam  carrying  this  pattern  will  be  travelling  parallel  to 
the  optical  axis  once  again  and  will  pass  through  the  second 
delay  line  DL2,  which  is  expected  to  be  identical  to  the 
first  delay  line.  Bragg  modulation  will  occur  in  this  delay 
line  as  in  the  first,  and  the  product  beam  will  then  pass  to 
the  output  detector.  Imaging  from  the  LCLV  to  DL2  will  occur 
by  means  of  a  telecentric  optical  system  identical  to  the 
first  section,  using  two  more  Wollensak  lenses.  Again,  no 
stop  will  be  required  in  the  frequency  plane,  because  of  the 
Bragg  operation  of  DL2. 

Of  course,  in  closed  loop  operation,  the  output  of  the 
final  detector  will  be  electronically  subtracted  from  the 
main  channel  signal  which  was  used  to  modulate  the  laser 
source,  to  complete  the  adaptive  loop.  As  pointed  out  in 
Section  4.1.4  all  of  the  operations  described  are  theoreti¬ 
cally  linear.  With  the  introduction  of  the  two  diffracted 
orders  from  the  Bragg  line  to  the  image  intensifier,  a  linear 
carrier  frequency  is  realized  in  the  intensity  variation. 

The  image  intensifier  is  sensitive  to  input  intensity  and 
converts  this  to  an  emitted  intensity.  This,  in  turn  is 
preserved  as  an  intensity  variation  into  the  photoconductor 
of  the  LCLV.  Here  the  linear  intensity  variation  is  scaled 
into  a  coherent  light  amplitude  modulation  of  the  laser  beam 
passing  through  the  liquid  crystal.  This  is  multiplied  in 
amplitude  by  the  diffraction  effect  of  DL2,  and  the  final 
product  is  summed  and  delivered  to  the  final  detector.  This 
detector  is,  like  the  intensifier,  sensitive  to  input  inten¬ 
sify,  and  the  carrier  realized  by  the  interference  of  the  two 
ditfracted  orders  becomes  a  linear  function  in  the  detector 
current  output.  It  is  obvious  that  the  loop  gain  can  be 
controlled  over  wide  limits  by  electronic  control  in  the 
feedback  path. 


6  AOAP  SYSTEM  TEST  STUDY 


6.1  TRANSIENT  RESPONSE  OF  AN  ADAPTIVE  LINEAR  PREDICTION 

FILTER  TO  A  SINUSOID  IN  NOISE 

6.1.1  Introduction 

Linear  prediction  filters  have  many  applications  in 
diverse  fields,  such  as  spectral  analysis,  system  modeling, 
speech  encoding,  and  noise  cancellation  [18] -[21].  One  ap¬ 
plication  of  interest  is  the  cancellation  of  a  sinusoidal 
interference  in  white  noise.  A  comprehensive  steady-state 
analysis  of  this  problem  appears  in  [2  „  Of  interest  here, 

however,  is  the  transient  response  to  sinusoid  in  noise. 

In  particular,  the  behavior  of  the  tre  -lent  response  will  be 
examined  as  the  number  of  taps  N  appn  hes  infinity,  but  the 
total  delay  remains  constant.  The  mo  /  ion  for  this  inves¬ 
tigation  stems  from  a  desire  to  charac  ize  the  transient 
performance  of  the  AOAP  which  realizes  a  continuum  of  adap¬ 
tive  weights  [4]  —  [6 ) .  The  transient  response  in  the  noise- 
free  case  was  described  on  a  very  simplified  basis  in  Section 
2-2,  and  the  conclusions  here  will  reinforce  the  essentials 
of  that  discussion. 

A  general  discussion  of  the  ALP  filter  and  its  analyti¬ 
cal  properties  will  be  presented  first,  followed  by  a  tran¬ 
sient  analysis  for  the  case  of  a  sinusoid  in  white  noise.  A 
particular  example  will  then  be  considered  and  the  asymptotic 
behavior  examined  as  the  number  of  taps  increases  to  infini¬ 
ty.  Experimental  results  are  also  obtained  using  a  real  time 
ALP  filter  which  was  previously  designed  and  constructed 
[18]. 


6.1.2  General 


Figure  6-1  shows  an  adaptive  transversal  filter  struc¬ 
ture  whereby  an  output 

N 

y(t)  =  x ( t )  -  l  W  x ( t-T  )  (6-1  ) 

n=  1  n  n 

is  formed  by  subtracting  from  the  input  x(t),  a  sum  of  N 
delayed  input  signals  x(t-Tn),  n  =  1,2,...,N,  weighted  by 
the  adaptive  weights  WlfW2,...WN.  The  weighted  sum  is 
denoted  by  x(t)  and  is  intended  to  form  a  least-mean-square 
estimate  of  the  present  value  of  x(t)  using  a  linear  sum  of 
past  values.  Hence,  this  configuration  is  known  as  an  adap¬ 
tive  linear  predictor  (ALP)  filter  [18].  The  adaptive  weight 
algorithm  by  which  this  is  accomplished  is  known  as  the 
"least-mean-square"  (LMS )  algorithm  [19],  which,  in  discrete 
time,  can  be  written  [18] 


Figure  6-1.  Block  Diagram  of  Adaptive  Linear  Prediction 
Filter  Using  Correlation  Cancellation  Loops 


Wn(t+1)  =  Wn(t)  +  u  y(t)x(t-Tn),  n= 1 , 2 , . . . , N .  (6-2) 

As  is  well-known  in  linear  estimation  theory,  LMS  estimation 
error  requires  that  the  error,  or  output  y(t),  in  this  case, 
be  uncorrelated  with  each  auxiliary  signal,  and  is  known  as 
the  "orthogonality  condition"  [23].  This  is  exactly  the  con¬ 
dition  that  (6-2)  converges  to  in  mean,  and  hence  is  also 
known  as  a  "multiple  correlation  cancellation  loop"  [18]. 
Substituting  (6-1)  into  (6-2)  and  solving  for  this  condition 
then  determines  the  mean  steady-state  weight  vector 

E  W»  =  R-1  r  (6-3) 

where  R  is  an  N  *  N  covariance  matrix  with  elements 


=  E  {  x  ( t-T £ ) x  ( t-T  j  ) } 


(6-4) 


and  r  is  an  N  *  1  column  vector  with  components 


ri  “  E  {  x ( t ) x ( t-T £ ) }  .  (6-5) 

The  mean  transient  solution  of  (6-2)  can  be  written  [18] 

E  W(t)  =  (I-yR)  fc  W(0)  +  R-1[  I_-(I_-yR)  fc]  r  ,  (6-6) 

which  can  be  shown  to  converge  to  (6-3)  provided 

0  <  “  <  2/xn,a*'  <6‘7> 

where  Xmax  is  the  largest  eigenvalue  of  the  covariance 
matrix  R.  The  eigenvalues  of  R  also  characterize  the  tran¬ 
sient  solution  (6-6),  each  one  being  associated  with  an 

-yX  .  t 

exponential  convergence  mode  of  the  form  e  1  [19]. 

Thus,  the  convergence  time  constant  of  the  ifch  mode  is 
given  by 


where  Ts  is  the  sampling  interval. 

In  many  cases  of  practical  interest,  the  so-called 
"clipped-LMS"  algorithm  is  employed  whereby  (6-2)  becomes* 

W^(t+1)  =  W^(t)  +  y  y ( t )  sgn  [x(t-Tn)].  (6-9) 

In  this  variation,  the  correlation  multiplication  only  in¬ 
volves  the  sign  of  the  auxiliary  tap  signal,  and  results  in  a 
simplified  hardware  implementation.  It  can  be  shown  that  for 
sinusoidal  signals  and/or  Gaussian  processes,  the  weights 
will  converge  to  exactly  the  same  steady  state  values,  pro¬ 
vided  that 

y  <  - - -  (6-10) 

J  2*  *max 

IT  O  „ 


*  An  apostrophe  (’)  will  be  used  throughout  to  indicate 
quantities  associated  with  the  clipped-LMS  algorithm. 


However,  the  modal  time  constants  in  this  case  are  modified 
to  [24] 


/!> 


(6-11  ) 


6.1.3  Transient  Response  for  Sinusoid  in  Noise 


An  analysis  of  the  transient  response  of  an  ALP  filter 
to  a  sinusoid  in  white,  or  uncorrelated,  noise  appears  in 
[25]  and  will  be  summarized  here.  A  single  sinusoid  of  fre¬ 
quency  fQ  =  uQ/2n  in  white  noise  can  be  expressed  as 


x(t)  *  a  cos  uQt  +  n(t). 


(6-12) 


where  n(t)  is  white  noise  with  variance  a2.  Substituting 
this  expression  into  (6-4)  and  (6-5)  then  gives  the  elements 
of  the  covariance  matrix  as 


‘ij  =  Y  COs[a>o(Ti-Tj)]  + 


(6-13) 


rA  =  y  cos  (0)oTi) , 


(6-14) 


A  o 

where  y  =  A  /2  is  the  power  of  the  sinusoid.  If  the  taps 


are  uniformly  spaced  as 


Tn  =  T0  +  nT  r 


(6-15) 


the  eigenvalues  of  (6-8)  can  be  simply  expressed  as  [24 


x1,2  -  *5  <u8>  +  °2 


( 6-16a) 


Xi  =  o  ,  2<  i  <  N 


( 6-1 6b ) 


where 


sin  Nut 
o 


N  sin  uQt 


(6-17) 
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It  is  also  shown  in  [24]  that  the  sinusoidal  component  of  the 
output  corresponds  to  the  convergence  modes  associated  with 
the  two  largest  eigenvalues  (6-1 6a).  If  6  *  0,  then  the 
composite  convergence  will  be  dominated  by  the  smaller  of 
these  two  eigenvalues.  Therefore,  combining  (6-8)  and  (6- 
16a)  gives  the  dominant  time  constant  as 


T max  "  „  [-fL(i-|B|)  +  a2]  ‘  (  " 

For  the  clipped-LMS  algorithm,  a  similar  expression  is  ob¬ 
tained  from  (6-11)  as 


t  ' 

max 


+  a 


2 


( 1  -  1 6  |  ) 


a2] 


(6-19) 


Special  cases  of  the  above  occur  for  a  noiseless  sinu¬ 
soid,  for  which 


T  max 


T 

s 


pN  ^  ( 1- 1  8 


) 


and 


T  ' 

max 


‘N  J  W  (1-|6  |) 


a 

n 


0  . 


(6-20) 


(6-21  ) 


Examination  of  these  last  two  expressions  incidentally 
shows  how  the  clipped-LMS  algorithm  reduces  the  variation  of 
time  constant  with  Y . 


It  does  not  appear  to  be  analytically  tractable,  at 
least  within  the  scope  of  the  present  effort,  to  extend  this 
analysis  to  the  case  of  a  sinusoid  in  correlated  noise. 
Therefore,  examination  of  this  effect  will  be  deferred  to  the 
section  on  experimental  results. 

6.1.4  An  Example 

An  example  will  now  be  given  to  illustrate  the  foregoing 
theory.  Suppose  that  the  tap  delays  are  chosen  such  that 
foTi=1/8  and  f0TN=3/8,  i.e.,  the  taps  span  90  degrees 
of  phase,  as  illustrated  in  Figure  6-2.  Then,  the  spacing 
between  taps  T  is  determined  by  the  expression 


Figure  6-2. 


Phasor  Diagram  of  Main  and  Auxiliary  Tapped 
Delay  Line  Signals 


f  T 


1 


4 (N— 1  ) 

Substitution  into  (6-17)  then  gives 


(6-22) 


0  = 


sin  [ - J 

2(N-1 y 

N  sin  ( 


2  (N- 1  ) 


) 


r 


2 

U 


(6-23) 


Equation  (6- 16a)  shows  that  in  the  limit  as  N  ♦  • ,  the 
eigenvalues  increase  to  infinite  values,  which  will  eventual¬ 
ly  result  in  instability  according  to  (6-8)  or  (6-11). 
Therefore,  in  reality  the  gain  must  be  reduced  as  N  increases 
so  as  to  keep  the  value  of  yN  a  constant.  This  is  also  a 
natural  physical  constraint  for  an  implementation  involving  a 
continuum  of  weights  such  as  AOAP.  It  can  be  seen  from  ( 6- 
20)  and  (6-21)  that  if  yN  is  a  constant,  then  the  time  con¬ 
stant  for  the  noiseless  sinusoid  varies  with  N  only  as  a 


function  of  0  which  approaches  the  asymptotic  value  of  2/ir 
given  in  (6-23). 


Consider,  for  example,  the  clipped-LMS  algorithm  operat 
ing  with  a  noiseless  (o=0)  sinusoid  of  rms  value  0.1  corre¬ 
sponding  to  y  =  0.01,  a  loop  gain  y N  *  1/160  and  sampling 
frequency  of  6500  Hz.  Substitution  of  these  values  into  (6- 
23)  and  (6-21)  then  gives  the  values  in  Table  X  for  various 
values  of  N. 


TABLE  X.  CALCULATED  VALUES  FOR  y  «  0.01,  an  =  0, 
Ny  *  1/160,  Tg  *  1/6500 


N 

6 

*min 

T  max 

2 

0 

0.0100 

0.6170 

4 

.4330 

0.0113 

1  .0882 

8 

.5477 

0.0181 

1.3641 

16 

.5946 

0.0324 

1.5222 

00 

.6366 

00 

1  .6980 

As  expected,  the  time  constant  approaches  a  limiting  value  as 
N  ♦  6  . 

6.1.5  Experimental  Results 

The  ALP  filter  hardware  described  in  [18]  was  restored 
to  operation  and  was  used  to  obtain  the  experimental  results 
reported  here.  The  purpose  of  these  experiments  was  to  veri¬ 
fy  the  analytical  results  obtained  in  the  last  section  for 
uncorrelated  noise  and  to  determine  the  effects  of  correlated 
noise  on  this  performance. 


Figure  6-3  shows  the  experimental  transient  response 
obtained  using  the  same  parameters  as  the  example  in  the  last 
section  (f0=446.43  Hz,  T,=0.28  ms,  Tn=0.84  ms,  y=0.01, 
yN»1/160,  Ts=1/6500).  Tne  upper  trace  shows  the  adaptive 
weight  trajectories  for  N=8  and  N=2,  respectively,  referred 
to  a  10V  reference  (2  cm/unit).  These  curves  can  be  seen  to 
approach  the  theoretical  steady  state  values  given  in  [22]. 
For  N*2,  these  values  are  simply  ±  /  2/2  2  ±  0.707,  as  can  be 
appreciated  from  the  diagram  in  Figure  6-2. 


Figure  6-3.  Transient  Response  for  yN=1/160,  Noiseless  Case 


The  lower  trace  in  Figure  6-3  shows  the  transient  re¬ 
sponse  of  the  output,  which  starts  initially  at  the  input 
level  of  0.1V  rms  s  282  mV  peak-to-peak* ,  and  decays  down  to 
a  residual  noise  level  caused  by  digital  processing  noise 
which  was  inadvertently  picked  up  (for  a  noiseless  sinusoid, 
the  output  voltage  should  ideally  converge  to  zero). 

Both  of  the  exponential  curves  in  Figure  6-3  exhibit  a 
time  constant  which  can  be  seen  to  be  in  rough  agreement  with 
the  theoretical  values  in  Table  I  for  N=8  and  N=2. 

For  the  remaining  experiments,  it  was  more  convenient  to 
maintain  the  same  loop  gain  which  was  held  constant  at 
ii  =  1/32t>.  Figure  6-4  shows  a  replay  of  Figure  6-3  with  con¬ 
stant  loop  gain,  and  will  be  used  to  compare  with  plots  to 
follow. 

Figure  6-5  shows  the  convergence  for  the  signal  plus 
noise  case.  The  noise  level  was  31.6  mV  rms  (-30  dBV)  which 
corresponds  to  a  10  dB  SNR.  The  noise  bandwidth  for  this 
equipment  was  5  kHz,  which  resulted  in  nearly  uncorrelated 


*The  calibration  shown  for  the  lower  traces  in  Figures  6-3 
through  6-6  is  in  error;  the  actual  value  should  be  100 
mV/ cm. 
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6.2  STEADY-STATE  PERFORMANCE  OF  A  REAL  ADAPTIVE  PROCESSOR 
IN  A  MULTI-PATH  ENVIRONMENT 

6.2.1  Introduction 


The  multiple  "correlation  cancellation  loop"  (CCL)  [1],- 
[2]  or  "least-mean-square"  ( LMS )  algorithm  [3]  [26]  [27] 

technique  is  widely  employed  in  adaptive  arrays  for  radar, 
sonar,  and  communication  systems  as  well  as  in  many  other 
signal  processing  applications,  such  as  adaptive  data  com¬ 
pression  [28] ,  adaptive  echo  cancellation  [29] ,  adaptive 
equalization  [30],  adaptive  noise  cancellation  [19],  adaptive 
linear  prediction  [18]  and  adaptive  multipath  cancellation 
[31]. 


Typical  radar  systems  employ  multiple  CCL's  as  sidelobe 
cancellers  to  null  the  effects  of  interfering  signals.  In 
some  cases,  the  interfering  signals  are  received  over  multi¬ 
ple  transmission  paths  due  to  reflections  from  terrain  fea¬ 
tures  and/or  man-made  objects  such  as  towers  or  buildings. 
Sidelobe  cancellation  in  the  multipath  environment  is  then 
complicated  by  the  increased  degrees  of  freedom  required  to 
cancel  each  path.  However,  in  this  case,  the  multiple  re¬ 
ceived  signals  are  correlated  in  time  and  can  therefore  be 
removed  by  incorporating  tapped  delay  lines  in  each  auxiliary 
channel  and  providing  an  adaptive  cancellation  loop  for  each 
tap. 


The  question  naturally  arises  as  to  what  length  of  delay 
line  and  how  many  taps  are  required  for  effective  cancella¬ 
tion  of  multipath  narrowband  or  wideband  interference.  Sev¬ 
eral  ad-hoc  studies  have  been  performed  to  partially  answer 
this  question  for  specific  applications.  However,  no  one  to 
the  authors'  knowledge,  has  undertaken  a  definitive  study  of 
the  problem  or  has  derived  any  general  guidelines  to  aid  in 
this  question.  It  is  the  purpose  of  this  report  to  develop  a 
methodology  for  dealing  with  this  important  problem  and  to 
provide  a  simple  model  for  understanding  the  behavior  of 
cancellation  performance  and  how  it  relates  to  the  various 
parameters  involved. 

It  is  emphasized  that  this  study  concerns  the  behavior 
of  a  CCL  processor  for  real  signals  as  opposed  to  complex 
baseband  signals  derived  from  in-phase  and  quadrature  demodu- 
lation  of  a  narrowband  signal,  as  is  usually  encountered. 

The  reason  for  this  departure  relates  to  the  optical  configu¬ 
ration  of  the  AOAP,  in  which  the  optical  analog  signals  are 
inherently  real  quantities.  The  "continuum"  of  signals  or 
weights  are  contained  on  an  optical  traveling  spatial  carrier 
frequency,  reflecting  the  character  of  the  narrow  band 


signals  being  processed.  In  this  way,  the  real-signal  opti¬ 
cal  channel  becomes  equivalent  to  a  complex  signal  processor, 
handling  the  equivalent  of  in-phase  and  quadrature  compo¬ 
nents. 


As  will  be  seen,  the  chief  distinction  in  this  case  is 
that  more  delay  line  taps  must  be  employed  for  adequate  quad¬ 
rature  sampling.  This  is,  however,  no  hardship  to  bear  for 
the  AOAP  which  realizes  an  almost  infinite  effective  number 
of  taps,  limited  only  by  optical  resolution. 

The  organization  of  this  section  is  as  follows:  A  brief 
description  of  the  multiple  CCL  technique  is  given  in  Sec¬ 
tion  6.2.2.  In  Section  6.2.3,  adaptive  cancellation  is 
studied  for  wideband  and  narrowband  jammers  in  a  single  mul¬ 
tipath  environment.  An  impulse  response  model  is  also  devel¬ 
oped  in  order  to  qualitatively  explain  the  performance. 
Section  6.2.4  presents  a  few  representative  computer  simula¬ 
tion  results  and  Section  6.2.5  contains  conclusions  and  gen¬ 
eralizations  of  the  multiple  CCL  simulation  study.  Several 
appendices  are  included  for  various  special  analysis  topics, 
additional  simulation  plots,  and  computer  program  listings. 

6.2.2  General  Theory 

A  diagram  of  a  multiple  real  CCL  processor  to  be  used 
for  interference  cancellation  in  a  radar  or  sonar  system  is 
shown  in  Figure  6-7.  For  simplicity,  only  a  single  auxiliary 
channel  with  multiple  taps  is  considered  here?  the  techniques 
and  results  derived  can  be  readily  extended  to  the  multiple 
auxiliary  case. 

The  output  signal  r(t)  is  formed  by  subtracting  from  the 
main  input  m(t)  a  weighted  sum  of  N  delayed  versions  of  the 
auxiliary  signal  a(t),  viz., 

N-1 

r ( t )  =  m(t)  -l  wA  a ( t  -  iT/N),  (6-24) 

i=0 

where  T  is  the  length  of  the  delay  line.  Cancellation  of  the 
main  signal  is  achieved  by  adjusting  each  adaptive  weight  so 
as  to  reduce  the  correlation  between  the  corresponding  de¬ 
layed  auxiliary  signal  and  the  residue.  This  procedure  can 
be  written  as 


MAIN  INPUT 


AUX  INPUT 


a  ('$)  •••  a(t(VT) 


OUTPUT 

RESIDUE 


-.-v- - 1 


WN-1 


Figure  6-7.  Adaptive  Array  Processor  with  a  Single 
Auxiliary  Tapped  Delay  Line 


where  p  is  the  loop  gain  constant.  Adaptation  then  proceeds 
until  the  .residue  is  uncorrelated  with  each  auxiliary  signal. 


i.e. 


E  { r ( t )  a( t-iT/N)}  =  0, 


i  =  0,1,... , N-1 .  (6-26 ) 


This  process  is  equivalent  to  applying  the  orthogonality 
principle  [23]  which  minimizes  the  rms  residue  power.  Sub¬ 
stitution  of  (6-24)  into  (6-26)  then  determines  the  optimal 
steady  state  adaptive  weights  (denoted  by  the  subscript  « )  as 
the  solution  to 
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(6-27) 


l  rAj  w„ j  =  i  =  0,1,... ,N-1 

j=0 


where 


r  ,  IT  .  ,  JT  .  ,  A  .  .  T  , 

r4j  -  E  [a(«;-  —  )  a(t-  —  )]  ■=  Ra[  U-3)  -  ] 

iT  A  iT 

ri  =  E  [m(t)  a(  t-  — )  ]  =  Rma  (-  —  )  , 


( 6-28a ) 


( 6-28b ) 


and  Ra(x  )  and  R^fT  )  are  the  autocorrelation  and 
cross-correlation  functions  of  the  processes. 

In  matrix  notation  (6-27)  becomes 


R  w„  =  r 


(6-29) 


where  R  is  the  N  *  N  auxiliary  signal  covariance  matrix  with 
elements  r^j  and  r  is  N  *  1  column  matrix  with  elements 
ri/  appropriately  defined  in  (6-28).  Since  uniform  tap 
spacing  was  assumed,  the  matrix  R  is  Toeplitz,  i.e.,  its 
element  values  r^j  are  a  function  of  i-j  only.  Also, 
since  R  is  a  covariance  matrix,  it  is  positive  definite  and 
hence  nonsingular.  Therefore,  a  unique  solution  exists  for 
the  weight  vector  which  can  be  obtained  as 

Woe  =  R-1  r.  (6-30) 

The  most  common  computational  approach  in  solving  (6-29) 
is  the  Levinson-Robinson  technique  [32]  which  requires 
only  O (N2 )  computations,  and  is  the  one  employed  here.  How¬ 
ever,  for  large  N,  R  becomes  so  ill-conditioned  that  computa¬ 
tions  even  with  double  precision  computer  arithmetic  can  at 
times  lead  to  inaccurate  results.  Therefore,  care  must  be 
exercised  in  applying  these  calculations. 

6.2.3  Adaptive  Cancellation  in  a  Multipath  Environment 

6 .2 .3 . 1  General 

A  major  goal  of  this  report  was  to  define  a  multipath 
interference  cancellation  test  simulation  that  encompasses 
all  the  important  parameters  of  the  problem,  yet  is  simple 
enough  so  that  the  effect  of  various  parameters  on  the  multi¬ 
ple  CCL  processor  cancellation  performance  can  be  effectively 
studied. 


't. 


WS'i 


A  simplified  system  diagram  in  the  presence  of  a  wide¬ 
band  or  narrowband  jammer  in  a  single  multipath  environment 
is  shown  in  Figure  6-8.  Each  weight  w^  is  calculated  by  a 
correlation  cancellation  loop  as  shown  in  Figure  6-7,  such 
that  the  output  signal  r(t)  is  uncorrelated  with  the  delayed 
auxiliary  signals  a(t-iT/N),  i=0 , 1 , . . . ,N-1 .  Without  loss  of 
generality,  the  noise  entering  the  main  channel  can  be  ne¬ 
glected  since  it  will  be  uncorrelated  with  the  auxiliary 
signals  and  therefore  will  propagate,  unaffected,  to  the 
output  of  the  processor.  The  processor  receiving  front  end 
contains  a  bandpass  filter  of  bandwidth  B  and  center  frequen¬ 
cy  fp  such  that  the  noise  n(t)  entering  the  auxiliary  chan¬ 
nel  is  colored  with  power  spectrum  Sn(f)  shown  in  Figure  6- 
8  and  given  as 


sn(f) 


1 

B 


rect  ( 


f-f 

B 


( 6-31  a ) 


where  rect(x)  is  defined  as 


A  1 
rect(x)  =  { 

0 


I  x | <  1/2 
I  x | >  1/2 


Taking  the  inverse  Fourier  transform  of  the  power  spectrum 
then  gives  the  noise  autocorrelation  function 


Rn(T  )  =  sine  Bt  cos  o>0t  ,  ( 6-31b) 

where  oj0  =  2u  fQ,sinc(x)  is  defined  as 
A  sin  it  x 

sinc(x)  =  -  , 

IT  X 

and  for  convenience  the  noise  power  is  normalized  to  unity, 
i .  e . ,  Rn ( 0 )  =  1 . 

The  received  signal  in  the  main  channel  m(t)  is  given  by 


m(t)  =  s ( t )  +  <*m  s(t-TQ)  (6-32) 

where  am  is  the  multipath  component  strength  (normally 
aq,<1  )  and  t0  is  the  excess  time  delay  relative  to  the 
direct  path.  The  corresponding  received  signal  in  the  auxil¬ 
iary  channel  is  given  by 

a ( t)  =  s ( t )  +  o a  s(t-T0— t,)  +  n(t),  (6-33) 
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Wideband  Jammer 


Narrowband  Jammer 


Figure  6-8.  Simplified  Adaptive  Array  Processor  in 
Multipath  Environment.  (a)  Block  Diagram,  (b)  Noise  and 
Jammer  Spectra,  and  (c)  Autocorrelation  Functions 
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where  Pm  »  R^O)  is  the  power  in  the  main  channel/  and 
Ra  and  Rma  are  defined  in  (6-36)  and  (6-37).  In  vec¬ 
tor  form,  (6-38)  can  be  expressed  as 

Pr  -  Pm  +  wT  (R  w  -  2  r ) , 

where  use  of  (6-28)  was  employed.  The  residue  power  is  mini¬ 
mized  for  w  *  Ho  as  defined  in  (6-29)  giving 


Pr  .  -  Pm  -  3*  r  . 

rmin  m  — 

6. 2. 3. 2  Wideband  Jammer 


(6-39) 


For  a  wideband  jammer,  the  interfering  signal  s(t)  can 
be  represented  as  a  random  process  with  a  flat  bandpass  power 
spectrum 


S_(f)  =  -  rect  ( 

s  B 


( 6-40a ) 


as  shown  in  Figure  6-8.  Its  corresponding  autocorrelation 
function  Rs(t),  which  is  the  inverse  Fourier  transform 
of  Ss(f),  is  given  by 


rs(t)  =  y  sine  Bt  cos  uiQr  . 


6. 2. 3. 3  Narrowband  Jammer 


( 6-40b ) 


Under  the  narrowband  jammer  assumption,  the  interfering 
signal  s(t)  is  represented  as  two  sinusoids  of  corresponding 
signal-to-noise  power  ratios  y1  ,  y2 ,  that  is. 


s  ( t )  =  / 2y  .  cos  a) .  t  + 


2y  j  cos  w9t  , 


(6-41  ) 


where  fl  =  u^/2*,  f,  =  u2/2t<  are  frequencies  contained  in  the 
bandpass  filter  of  bandwidth  B.  The  corresponding  one-sided 
power  spectrum  is  given  as 


Ss(f)  =  t  ( f-f  t  )  +  y  2s  ( f-f  2 ) 


( 6-42a ) 


Taking  the  inverse  Fourier  transform  of  the  two-sided  power 
spectrum,  the  autocorrelation  function  Rs(t )  is  given  as 


Rg(T)  =71  COS  (i>  i  x  +  Y2  COS  U>2T 


( 6-42b) 
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The  total  output  power  can  be  obtained  as  in  ( 6—38 ) — ( 6- 
39).  However,  in  this  case,  an  alternate  expression  can  be 
obtained  by  substituting  (6-32),  (6-33),  (6-41)  in  (6-38a) 
and  taking  Fourier  transforms,  giving 


Py.  =  Py.  +  P_  +  P  y. 

r  r1  r2  rc 


(6-43) 


where 


-j0,kTo  N-1  -j«k^  ■3“k(l+To+Tl)112 

Prk^kl1+anie  -  l  wite  +  aae  ]|  » 

K  i=0 


k=  1  ,2 


(6-44) 


are  sinusoidal  components  of  the  output  power  and 
N-1  N-1 

Pr  -  l  l  w£  w,  Rn  [ ( i- j )  T/N]  (6-45) 

°  i=0  j=0 

is  the  noise  output  power.  The  function  Rn  (t )  in  (6-45) 
is  evaluated  using  (6-31). 

6. 2. 3. 4  Impulse  Response  Model 

An  equivalent  representation  of  (6-32),  (6-33)  is  given 
as 


m(t)  =  s ( t )  *  hm(t) 

a(t)  *  s ( t )  *  h  ( t)  +  n(t)  , 


where 


hm(t)  -  6  (t)  +  aai(t-T0) 
ha(t)  -  «  (t)  +  a a5 ( t_T o“T  1 ) 


( 6-46a ) 
( 6-46b ) 

( 6-47a ) 
(6-47b) 


represent  the  corresponding  main  and  auxiliary  channel  im¬ 
pulse  responses.  The  corresponding  channel  frequency  re¬ 
sponses  Hro(f),  Ha(f)  are  the  Fourier  transforms  of 
hm(t),  ha(t)  in  (6-47),  and  are  given  as 


Vf)  -  1  +  <>me 


(6-48a) 


Ha(f)  =  1  +  aae 


-j“(T0+  T!  ) 


( 6-48b) 


where  w  =  2*  f . 


An  ideal  processor  with  knowledge  of  the  channel  fre¬ 
quency  responses  Hm(f)  and  Ha(f)  would  process  the  auxil¬ 
iary  signal  a(t)  by  first  applying  the  inverse  filter 
Ha-'  in  order  to  remove  the  effect  of  multipath  in  the 
auxiliary,  and  then  apply  the  filter  H^,  to  introduce  the 
multipath  components  associated  with  the  main  signal  m(t). 
Thus,  the  composite  filter  that  would  be  applied  to  the  aux¬ 
iliary  signal  can  be  written 


H(f)  =  H  "'(f)  H_( f ) , 


( 6-49a ) 


or  in  the  time  domain 


h(t)  =  ha_1(t)  *  hm(t). 


( 6-49b ) 


Finally,  subtraction  of  the  resultant  signal  from  m(t)  would 
give  the  output  signal  r(t)  as  shown  in  Figure  6-9.  Also 
shown  are  the  corresponding  filter  impulse  responses  hm(t), 
ha(t),  ha-1  (t)  and  h(t),  where 


*>a  (fc)  =  F_1[Ha  7(f)]  =  I  ( -a a ) n  6[  t-n(T0+Tl  )]  . 

n=Q 


(6-50) 


For  t1=o,  substitution  of  (6-50)  and  (6-47a)  into  (6-49b) 
gives 


h  ( t)  =6(t)  +  (om-aa)  l  { -a  a)  n~ 1  5(t-nx0)  .  (6-51) 

n=1 


This  response  will  be  appropriate  as  an  approximation  to  the 
ideal  filter  for  an  actual  system  if  Bt1 <<  1,  a  condition 
that  is  almost  always  realized  in  practice. 


The  actual  adaptive  CCL  processor  on  the  other  hand  has 
no  knowledge  of  the  channel  filters  Hu,  and  Ha  but  adapts 
the  weights  w^  such  that  the  estimate 


m(t)  =  l  wi  a(t - ) 
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Figure  6-9.  Ideal  Multipath  Cancellation  System 
(a)  Block  Diagram,  (b)  Filter  Impulse  Responses 


is  optimum  in  the  sense  of  minimizing  the  mean-squared  error. 
In  the  absence  of  multipath,  the  steady  state  performance  of 
the  processor  would  be  the  limiting  behavior  of  (6-39),  being 
a  function  of  the  jammer  signal  s(t),  and  the  number  of  adap¬ 
tive  taps  N.  In  the  presence  of  multipath,  the  processor  in 
essence  approximates  the  impulse  response  of  the  infinite 
impulse  response  (HR),  or  recursive,  filter  of  (6-51)  with 
the  finite  impulse  response  (FIR),  or  nonrecursive,  N  tap 
filter 


^  r  ,  IT, 

h ( t )  =  l  6(  t - )  .  (6-52) 

i=0  N 

/\ 

Thus,  the  problem  can  be  visualized  as  matching  h(t)  to  h(t) 
such  that  the  output  power  is  minimized. 

Comparison  of  (6-51)  with  (6-52)  shows  that  the  best 
filter  match  that  can  be  achieved  is  when  h(t)  =  h(t),  0<t<T. 
This  ideal  will  be  approximately  realized  if  the  number  of 
taps  N  is  large  enough  for  good  interpolation  of  the  multi- 
path  signal  components  located  at  multiples  of  tg*.  Thus, 
under  these  conditions,  and  neglecting  the  additive  noise, 
the  residue  will  be  approximated  by  the  expression 

r ( t )  =  m(t)  -  h(t)*a(t) 

-  [h(t)  -  h(t)]  *  a(t ) 

oo 

=  ( ani'a  a)  I  a(fc-nTo)  '  (6-53) 

n=M+l 

where  M  is  a  function  of  tq  defined  by 

Mtq<  T  <  (M+1  )tq  ,  ( 6-54a ) 

i. e. , 

M  =  [t/tJ  (6-54b) 

where  [  x]  is  defined  as  the  greatest  integer  less  than  or 
equal  to  x. 


*  An  analysis  of  interpolation  of  a  bandlimited  random 
process  appears  in  Appendix  B. 


The  residue  power  can  then  be  expressed  using  (6-53)  as 


Pr=  e{  r2  (t )} 

2  tam~aa)2  l  l  (“aa)n+m  RaKn_m)To]  (6-55a) 

n,m=M 


=  a 


(a  -a_) ' 
2M  a _ ® 


[P,+  2 


1  -a 


l  (-1 ) Pcosh 
P=1 


(p  A  na  a )  RQ ( pt  Q )  ]  , 
( 6-55b ) 


where  Ra  is  defined  in  (6-36). 


Examination  of  (6-55)  shows  that  this  expression  can  be 
usefully  applied  to  determine  most  of  the  behavioral  charac¬ 
teristics  of  adaptive  cancellation  in  multipath.  First  of 
all,  it  can  be  seen  that  in  the  unlikely  event  that 
am  =  aa»  residue  power  reduces  to  zero  as  would  be 

expected  since  both  signals  are  exactly  matched.  Further¬ 
more,  if  a a  =  0,  Pr  also  reduces  to  zero  as  would  be 
expected  since  in  this  case  there  is  no  inverse  auxiliary 
multipath  filter  to  approximate. 


A  special  case  of  the  above  for  the  wideband  jammer 
occurs  for  Btq  >>  1  and  =0,  for  which 


Ra[  (  n~m]  T0]  •  [(l+«a2)«mn+aa6m,n+1+aa6m,n-l]  ps' 

BTo^I, 


(6-56) 


and  therefore,  (6-55)  reduces  to 

pr  "  “a2”  <anr°a>2  ps  '  Bto»1  •  (6-57) 

In  this  case,  the  residue  power  will  increase  from  one  pla¬ 
teau  to  the  next  as  tq  is  increased,  where  each  successive 
level  change  occurs  with  an  integer  value  change  of 
M*[t/tq].  A  non-mathematical  way  of  viewing  this  behavior 
is  to  realize  that  the  approximating  FIR  filter  $  is  equiva¬ 
lent  to  the  HR  filter  h  depicted  in  Figure  6-9  over  a  window 
of  fixed  length  T.  Thus,  as  tq  is  increased,  the  number  of 
impulses  contained  in  the  window  will  fall  off  one-by-one. 
These  transitions  will  occur  for  t0/T  =  1,  1/2,  1/3, 

1/4,...,  i.e.,  where  M  changes  value. 


For  any  given  application,  (6-57)  can  be  utilized  to 
specify  the  length  of  delay  line  required  for  effective  can¬ 
cellation  of  multipath  interference  of  given  maximum  delay 
and  strength.  Note  that  it  is  generally  not  sufficient  to 
employ  a  delay  line  whose  length  merely  equals  the  longest 
multipath  delay.  As  an  example,  consider  a  multipath  of 
strength  aa  =  0.5  (-6  dB)  in  the  auxiliary  antenna  and  (as 
a  worst  case)  am  =  o  in  the  main  antenna.  For  these  param¬ 
eters,  the  cancellation  ratio  is  then 


^r  _  -2  [  T/t  Q+ 1  ] 

Pc 


(6-58) 


i.e.,  6  dB  for  each  multiple  of  T  to  tq.  Thus,  in  the  ab¬ 
sence  of  noise,  30  dB  of  cancellation  would  require 
T  >  5  t  i.e.,  the  required  delay  is  five  times  the  lon¬ 
gest  multipath  delay. 


6.2.4  Computer  Simulation  Results 
6. 2. 4.1  Wideband  Jammer 


A  computer  program  was  written  to  evaluate  the  perfor¬ 
mance  of  a  single  auxiliary  channel,  tapped  delay  CCL  proces 
sor  in  regard  to  its  ability  to  cancel  a  wideband  jammer 
signal  in  a  single  multipath  environment.  The  program  was 
written  in  FORTRAN  IV  code  and  was  compiled  and  executed  on 
VAX— II  computer.  All  plots  were  generated  on  a  PLXY  plotter 
Equations  (6-28),  ( 6-35 )-( 6-38 )  were  written  in  FORTRAN  code 
and  the  values  of  the  weight  vectors  in  (6-30)  were  obtained 
not  by  matrix  inversion  but  by  the  Levinson-Robinson  algo¬ 
rithm  which  requires  only  0  (N2 )  computations  [32] . 

Referring  to  the  simplified  system  diagram  of  Figure  6- 
8,  the  simulation  parameters  are  the  system  bandwidth  B, 
center  frequency  f^,  jammer  power  y  (also  power  to  noise 
ratio  since  the  noise  power  is  normalized),  multipath  compo¬ 
nents  strengths  am  and  aa,  excess  multipath  delay  in  the 
auxiliary  channel  ,  total  auxiliary  delay  T,  and  number  of 
adaptive  taps  N.  A  set  of  nominal  dimensionless  parameters 
values  was  selected  as  BT  =  2,  f0T  =  10,  t1/T  =  10~3 , 

Y  =  1000  (30  dB),  corresponding  to  B  =  20  MHz, 
fQ  =  100  MHz,  T  =  .1  Msec,  t 1  =  10“4  Msec*  and 
°a'  am  =  .5,  1.  Plots  of  Pr(dB)  vs. 

t 0/T  =  0,  .01,  .02,. ..,1.8  for  N  =  4,8,16,32,64  were  gen¬ 
erated. 


♦This  corresponds  to  an  angle  0=2°  in  Figure  6-8  for  an 
antenna  separation  of  1m. 


Initial  runs  using  single  precision  arithmetic  resulted 
in  poor  accuracy  for  N  >  8.  Therefore  all  calculations  were 
done  in  double  precision  arithmetic.  An  additional  problem 
was  experienced  for  large  values  of  N,  whereby  the  covariance 
matrix  R  of  (6-28a)  used  in  (6-29)  becomes  highly  ill- 
conditioned,  due  to  a  high  ratio  of  largest  to  smallest 
eigenvalues  of  R.  This  ratio  of  largest  to  smallest  eigen¬ 
values  is  sometimes  referred  to  as  the  condition  number  of  a 
matrix  and  can  be  related  to  the  spectral  variation  of  the 
sampled  power  spectrum  [33]  [20], 


Sa'(f)  =  l  Sa( f-nN/T )  .  (6-59) 

n=-» 

For  a  bandpass  spectrum,  there  is  no  overlap  in  (6-59) 
for  N  >  BT,  and  hence  the  eigenvalue  spread  is  unbounded  and 
the  matrix  is  inherently  ill-conditioned.  The  ratio  of  spec¬ 
tral  variation  will  also  be  large  for  oa  -  1  due  to  multi- 
path  cancellation  of  some  frequency  components. 

To  counteract  the  ill-conditioning  of  R  which  results  in 
inaccuracies  in  Pr,  a  small  value  AR  was  added  to  the  diag¬ 
onal  elements  of  R.  An  optimum  value  of 

AR  =  1 0“9 

was  determined  empirically. 

A  plot  of  Pr  (dB)  versus  tq/T  for  am=0,  aa=*.5 
and  N=4 ,8, 16,32 ,64  is  shown  in  Figure  6-10.  From  sampling 
theory,  it  is  known*  that  for  an  infinite  number  of  samples, 
optimum  performance  in  Pr  will  be  obtained  provided  the 
sampling  rate  (Rs  =  N/T)  is  greater  than  the  Nyquist  rate, 
which  is  defined  as  twice  the  largest  significant  frequency 
in  the  signal  spectra  (Rjj  =  2fQ  +  B).  The  deterioration 
in  performance  for  a  finite  number  of  samples  that  is  pro¬ 
gressively  getting  smaller  is  addressed  in  detail  in  Appendix 
B.  From  Figure  6-10  it  is  seen  that  for  t0/T<.3,  near¬ 
optimum  performance  is  attained  for  as  few  as  N=16  taps, 
corresponding  to  a  relative  sampling  ratio 

Rs  N  16 

-  =  -  =  -  =  0.73  . 

Rn  2  f0T  +  BT  2(10)  +  2 

Also,  the  performance  degradation  due  to  increasing  multipath 
delay  length  tq  is  evident  as  predicted  by  (6-55). 

*  The  relationship  of  sampling  theory  to  the  multipath 
cancellation  problems  was  suggested  by  Dr.  Frank  Dickey  of 
the  General  Electric  Co.,  Syracuse,  N.  Y. 


Figure  6-10.  Output  Residue  Power  vs.  Normalized  Multipath 
Delay  for  a  Wideband  Jammer  and  T  =  0.1  usee. 


The  fine  structure  in  the  response  is  due  to  spectral 
interference  of  the  multipath  components.  As  shown  in  Figure 
6-9,  the  auxiliary  channel  response  Ha(f)  contains  peaks  at 
integer  multiples  of  1/t0.  Therefore,  the  auxiliary  band¬ 
pass  signal  spectrum  will  be  enhanced  by  the  multipath  if 
t0/T  is  an  integer  multiple  of  1/f0T  =  0.1  and  hence 
improve  performance  at  these  points.  Thus,  the  valleys  in 
Figure  6-10  for  large  N  which  occur  for  tq/T  =  0.1, 0.2,... 
are  a  result  of  auxiliary  signal  enhancement.  On  the  other 
hand,  for  points  in  between  tq/T  =  .05,  .015,...,  auxiliary 
signal  is  reduced,  thereby  increasing  the  residue  power,  and 
is  also  evident  in  Figure  6-10. 

Another  plot  is  shown  in  Figure  6-11  for  all  the  same 
parameters  as  in  Figure  6-10  except  that  the  length  of  the 

delay  line  is  increased  by  a  factor  of  8  to  T  =  0.8  us.  It 

is  expected  that  the  Pr  response  for  N=64,128  will  be  at 
least  as  good  as  the  one  for  N=8,16  of  Figure  6-10  since  the 

sampling  rate  is  the  same,  although  more  total  taps  have  now 

been  added.  In  fact  it  is  noted  that  the  Pr  response  in 
Figure  6-11  has  now  achieved  its  limiting  behavior  for  N=64 
or  a  relative  sampling  ratio 

Rs  64 

—  =  -  =  0.36  , 

Rn  2(80)  +  16 
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Figure  6-11.  Output  Residue  Power  vs.  Normalized  Multipath 
Delay  for  a  Wideband  Jammer  and  T  -  0.8  vsec 


as  opposed  to  the  Figure  6-10  Pr  response  achieving  its 
limiting  behavior  at  Rs  =  .73  R*j  (N=16).  The  limiting 
performance  of  the  Pr  response  of  Figures  6-10  -  6-11  can 
be  explained  in  terms  of  the  taps  occurring  in  quadrature 
pairs.* 

The  first  four  cycles  of  the  fine  structure  in  Figure  6- 
11  are  similar  to  Figure  6-10  over  an  absolute  time  delay 
interval  of  0  <  To<  0«05»  as  would  be  expected.  For  larg¬ 
er  delay,  the  plateauing  effect  predicted  by  (6-57)  is  clear¬ 
ly  evident.  The  entire  response  in  fact  is  very  close  to 
that  predicted  by  the  impulse  response  model,  whereby  the 
residue  jumps  in  6  dB  steps  for  t0/T  =  1 , 1/2 ,1/3,1/4 , . . . 
and  eventually  reaches  a  level  of  1000/4  (24  dB).  Actually 
each  plateau  holds  on  a  little  longer  than  expected  because 
of  the  ability  to  predict  the  signal  value  somewhat  beyond 
the  length  of  the  delay  line. 


*  The  effect  of  quadrature  tap  spacing  on  the  processor 
performance  is  analyzed  in  Appendix  C. 
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6. 2. 4. 2  Narrowband  Jammer 

The  effect  of  multiple  narrowband  jammers  distributed 
over  a  band  of  frequencies  will  now  be  examined.  The  sim¬ 
plest  situation  to  consider  is  two  sinusoidal  jammers  located 
at  the  band  edges.  Referring  to  Figure  6-8,  the  additional 
simulation  parameters  required  to  specify  the  two  tone  jammer 
are  the  frequencies  fx  ,  f2  and  the  corresponding  signal 
strengths  (signal  to  noise  power  ratios)  Yi  »y2 •  The  nominal 
dimensionless  parameters  that  were  selected  are  f 2 T  =  9, 
f2T  =  11,  =  1000  (30  dB)  and  y,  =  100  (20  dB),  correspond¬ 
ing  to  fj  =  90  MHz,  f2  =  110  MHz  for  T  =  0.1  usee. 

A  plot  of  Pr  versus  tq/T  for  am  =  0,  =  .5, 

T=0.1  usee  and  N=4 , 8 , 1 6 , 32 , 64  is  shown  in  Figure  6-12.  The 
ripples  in  these  plots  can  be  explained  by  interference  of 
the  direct  and  multipath  sinusoids  at  the  main  and  auxiliary 
antennas:  variations  of  signal  level  at  the  main  input  will 

of  course  directly  translate  to  the  output;  on  the  other 
hand,  variations  at  the  auxiliary  input  will  affect  the  SNR, 
thereby  enhancing  or  degrading  the  cancellation  performance. 
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Figure  6-12.  Output  Residue  Power  vs.  Normalized  Multipath 
Delay  for  a  Narrowband  Jammer  and  T  =  0.1  usee. 


Another  point  of  interest  concerning  Figure  6-12  is  that 
the  plot  is  periodic,  repeating  at  tq/T  =  1.0  since  each 
sinusoid  has  then  gone  through  an  integer  number  of  cycles, 
viz.,  f^Q  =  9  and  f2x0  =  11. 

As  can  be  seen,  contrary  to  the  wideband  case,  the  per¬ 
formance  with  the  narrowband  jammer  does  not  degrade  as  t0 
is  increased.  This  behavior  is  a  consequence  of  the  periodic 
nature  of  the  signal  autocorrelation  function.  Therefore,  in 
this  case,  it  is  not  necessary  to  have  a  delay  line  length 
which  exceeds  the  longest  multipath  delay  as  was  the  case  for 
wideband  interference;  the  delay  only  has  to  encompass  one 
period  of  the  autocorrelation  function  envelope.  For  the 
case  of  two  sinusoids,  this  period  for  the  particular  param¬ 
eters  chosen  is  given  simply  by  T=2/(f2-f1  ),  which  is  equal 
to  0.1  y  sec  for  the  previous  example. 

Further  examination  of  Figure  6-12  shows  that  near¬ 
optimum  performance  is  attained  for  a  normalized  sampling 
rate 


2fJT  +  BT  2(10  )  +  2 


=  1.46 


corresponding  to  N  =  32,  as  compared  to  the  value 
RS/RN  *  *73  required  for  wideband  jammer  cancellation 
shown  in  Figure  6-10.  However,  in  this  case,  the  limiting 
residue  power  is  much  lower  since  the  sinusoidal  signals  are 
correlated  over  the  entire  delay  line,  thereby  distributing 
the  noise  power  over  a  larger  number  of  adaptive  weights.  In 
fact,  the  average  near-optimum  output  power  for  the  narrow- 
band  jammer  is  approximately  15  dB  lower  than  the  correspond¬ 
ing  wideband  power  for  approximately  equal  jammer  signal 
strength  (Yt  +  y2  -  30.4  dB,  whereas  y  =  30  dB).  All  of  this 
seems  to  suggest  that  wideband  jamming  is  the  most  difficult 
problem  to  be  solved,  presumably  due  to  the  higher  number  of 
degrees  of  freedom  involved. 

6.2.5  Conclusions 

The  steady-state  cancellation  performance  of  an  adaptive 
array  processor  using  a  single,  real,  auxiliary  tapped  delay 
line  with  a  multiple  CCL  has  been  analyzed  for  a  single  mul¬ 
tipath  interference  environment.  Both  wideband  random  noise 
jammers  and  multiple  narrowband  jammers  were  considered.  The 
techniques  and  results  developed  can  be  easily  extended  to 
multiple  auxiliary  channels  and  multiple  multipath  compo¬ 
nents  . 


For  all  types  of  interferences  good  cancellation  is 
generally  achieved  provided  that  the  auxiliary  signal  delay 
line  taps  are  spaced  closer  than  dictated  by  the  Nyquist 
sampling  rate.  For  the  wideband  jammer,  near-optimum  perfor¬ 
mance  was  attained  in  some  cases  at  half  the  Nyquist  sampling 
rate  or  less  due  to  the  taps  occurring  in  quadrature  pairs. 

For  small  multipath  delay,  Btq  <<  1,  the  CCL  processor 
responses  were  explained  in  terms  of  auxiliary  or  main  signal 
enhancement  or  reduction  due  to  multipath.  An  impulse  re¬ 
sponse  model  was  developed  to  qualitatively  explain  the  sys¬ 
tem  behavior  for  large  multipath  delay,  Bx0  >>  1.  This 
model  relates  the  performance  to  a  processor  which  employs  a 
finite  impulse  response  (FIR)  filter  in  the  auxiliary  channel 
to  match  an  ideal  infinite  impulse  response  (HR)  filter 
that  would  be  required  for  exact  cancellation. 

For  the  wideband  jammer,  good  cancellation  is  obtained, 
provided  that  the  length  of  the  delay  line  T  is  considerably 
longer  than  the  longer  multipath  delay  Xg,  depending  on 
multipath  components  strengths  of  the  main  and  auxiliary 
channels,  denoted  by  am  and  aa  respectively.  For  large 
multipath  delays,  Bt Q  >>  1,  the  limiting  output  residual 
power  can  be  simply  expressed  using  this  model  as 

pr  =  “a2”  (am"aa)2  ps  ' 

where  Ps  is  the  jammer  signal  power  and  M  =  [ T/x G]  is  the 
number  of  integer  multiples  of  xQ  contained  in  the  total 
delay  T. 

For  narrowband  jamming,  the  cancellation  performance  was 
shown  to  be  unaffected  by  x 0  provided  that  the  delay  line 
exceeds  the  period  of  the  autocorrelation  function. 

In  all  cases,  the  effect  of  the  auxiliary  channel  excess 
multipath  delay  length  t  ^  on  the  processor  performance  is 
minimal  since  it  merely  shortens  the  effective  delay  length 
to  T  -  Xj  =  T  for  x  1  <<  T,  and  is  negligible  for  any  prac¬ 
tical  application. 

The  main  results  of  this  report  concern  the  performance 
of  real  multiple  CCL,  which  relates  to  a  specific  hardware 
application.  In  general,  adaptive  cancellation  systems 
process  complex  signals  derived  from  in-phase  and  quadrature 
demodulation  of  a  narrowband  signal.  The  present  analysis 
approach  could  be  modified  to  accommodate  this  situation  by 
redefining  all  signals  as  complex  quantities  and  by  employing 
complex  arithmetic  in  solving  for  the  steady-state  residue 
power. 


APPENDIX  A 


ADDITIONAL  SIMULATION  RESULTS 


Additional  simulation  results  are  contained  in  this  Appendix 
which  further  demonstrates  the  behavior  of  the  processor  cancella 
tion  performance  in  a  multipath  environment,  as  a  function  of  the 
various  system  parameters.  Plots  of  processor  output  power  Pr 
(dB)  versus  normalized  multipath  delay  t0/T  for  a  wideband  jammer 
are  shown  in  Figures  A1-A16,  and  the  corresponding  values  of  the 
system  parameters  are  summarized  in  Table  AI . 

TABLE  AI.  WIDEBAND  JAMMER  SYSTEM  PARAMETERS 


Fig.#  f„  (MHz)  B (MHz ) 


T(y  sec ) 1 t  j /T 


IO-3  10“9  8,16, 

.  32,64 


+ 

0 

.  1 
.5 

1 

10-3 

10“3 

IQ"3 


64,128 


The  plots  in  Figures  A1-A9  demonstrate  the  effect  of  the  mul¬ 
tipath  component  strengths  for  am  =  aa  =  0,  .5,  1  on  the  Pr 
response.  As  discussed  in  the  main  body  of  this  report,  the  multi 
pie  CCL  processor  in  the  presence  of  multipath  approximates  the 
impulse  response  of  an  infinite  impulse  response  (HR)  deconvolu¬ 
tion  filter  with  one  of  a  finite  impulse  response  (FIR)  filter.  A 
predicted  by  (6-55)  and  demonstrated  by  Figures  A2-A3,  this 
approximation  is  hopeless  for  =  1,  since  the  deconvolution 
filter  has  a  constant  infinite  Impulse  response.  In  the  unlikely 
event  that  om  =  a-,  the  residue  power  contains  no  signal 
component  as  predicted  by  (6-55)  and,  thus,  reduces  to  the 
auxiliary  noise  level  as  seen  in  Figures  AI ,  A5  and  A9.  The 
similarity  of  the  Pr  responses  in  Figures  A4  and  A6  can  be  easily 


explained  by  (6-55)  since  otm-aa  is  the  same  for  both 
cases.  In  the  absence  of  the  auxiliary  multipath  component 
(aa«0),  the  processor  task  is  much  easier,  since  it  does 
not  now  have  to  synthesize  the  impulse  response  of  an  HR 
filter,  but  simply  introduces  the  main  multipath  component  in 
the  auxiliary.  Thus,  good  cancellation  performance  is  pos¬ 
sible  in  this  case  as  evidenced  in  Figures  A7-A8.  The  beha¬ 
vior  of  the  residue  power  in  Figures  A7-A8  for  x q/t  >  8  can 
be  explained  in  terms  of  good  extrapolation  being  possible 
provided  the  number  of  taps  N  is  large  enough*.  The  nominal 
parameter  values  am=o,  ot  ^,=  .5  were  selected,  and  the  re¬ 
maining  plots  are  variations  of  this  case. 


The  effect  of 
response  is  shown 
result  is  essentia 
tive  length  T-x,  . 
Figure  A6  shifted 
Figure  A10  plot  is 
plot.  Therefore, 
effect  of  x j  is  mi 
could  be  set  to  0. 


the  delay  across  the  array  x1  on  the  Pr 
in  Figures  A10-A13.  As  can  be  seen,  the 
lly  to  shorten  the  delay  line  to  an  effec- 
Thus,  the  last  three  plots  are  similar  to 
to  the  left  by  the  amount  Xj/T  and  the 
essentially  the  same  as  the  Figure  A6 
for  the  usual  system  condition  x,<<T,  the 
nimal  and  without  loss  of  generality  Xj 


The  response  of  Figure  A14  demonstrates  the  ill-condi¬ 
tioned  nature  of  the  covariance  matrix  R  of  (6-29)  for  large 
N  and  the  resulting  numerical  inaccuracTes  when  no  small  term 
( A  R=  0 )  is  added  to  its  diagonal  elements  to  correct  the  situ¬ 
ation.  Superposition  of  Figure  A14  on  Figure  A6  shows  that 
adding  ar=10~^  to  the  diagonal  of  R  has  introduced  an  error 
less  than  1  dB  for  large  xQ/T  and  significantly  less  than 
that  for  small  xQ/T. 

The  responses  of  Figures  A15-A16  are  similar  to  the  ones 
of  Figures  6-10  -  6-11  in  the  main  body  of  the  report,  and 
they  demonstrate  the  effect  of  increasing  the  length  of  the 
delay  line  by  a  factor  of  8  to  T=.8  y sec  for  the  system  pa¬ 
rameters  f  =  60  and  B  =  10  MHz.t  Note  again  that  as  predicted 
by  sampling  theory,  the  Pr  responses  in  Figure  A16  for 
N=64,128  are  at  least  as  good  as  the  ones  in  Figure  A15  for 
N=6,N=16.  Also  note  the  plateauing  effect  in  Figure  A16 
occurring  in  6  dB  steps  for  x0/T  s  1,  1/2,  1/3,  1/4,...,  as 
predicted  by  (6-57). 

The  final  plot  for  the  wideband  jammer  in  Figure  A17 
demonstrates  the  effect  of  adding  a  delay  line  of  length  T/2 
in  the  main  channel.  In  the  absence  of  auxiliary  multipath 
this  would  be  highly  desirable  in  order  to  provide  good 


*  This  is  one  of  the  results  of  Appendix  B. 

t  The  experimental  AOAP  laboratory  implementation  is  being 
designed  for  f  =  60,  B  =  10  MHz  and  T  =  7  Msec. 
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interpolation*  of  the  main  signal  multipath.  However,  for 
wideband  jammer  cancellation  with  auxiliary  multipath,  this 
modification  is  detrimental  since  the  effective  total  auxil¬ 
iary  delay  is  halved,  as  evident  by  comparison  of  Figures  A6 
and  A17. 

Additional  Pr  plots  for  a  narrowband  jammer  are  shown 
in  Figures  A18-A19  for  equal  multipath  component  strength 
(am=aa=1  and  am=aa=0  correspondingly).  The  ripples 
in  Figure  A18  can  be  explained  by  interference  of  the  direct 
and  multipath  sinusoids  at  the  main  and  auxiliary  antennas. 
In  particular  for  t0/T  =  .5  +  m,  where  m  is  an  integer,  the 
main  signal  entering  the  processor  is  zero,  resulting  in  a 
null  in  the  Pr  response.  Similarly  for  the  peaks  at  tq/T  = 
.59  +  m  and  tq/T  =  .61  +  m  are  due  to  auxiliary  signal  can¬ 
cellation. 
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Figure  Al .  Output  Residue  Power  vs.  Normalized  Multipath 
Delay  for  a  Wideband  Jammer  and  am=l •  “a=1 
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Figure  A2.  Output  Residue  Power  vs.  Normalized  Multipath 
Delay  for  a  Wideband  Jammer  and  am=.5,  aa=1 
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Figure  A3.  Output  Residue  Power  vs.  Normalizd  Multipath 
Delay  for  a  Wideband  Jammer  and  am=0,  aa=1 


Figure  A4.  Output  Residue  Power  vs.  Normalized  Multipath 
Delay  for  a  Wideband  Jammer  and  am=1 ,  aa=.5 
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Figure  A5.  Output  Residue  Power  vs.  Normalized  Multipath 
Delay  for  a  Wideband  Jammer  and  am=.5,  aa=.5 
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Figure  A6.  Output  Residue  Power  vs.  Normalized  Multipath 
Delay  for  a  Wideband  Jammer  and  am=0,  aa=.5 


Figure  A7.  Output  Residue  Power  vs.  Normalized  Multipath 
Delay  for  a  Wideband  Jammer  and  am=1 ,  aa=0 
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Figure  A8.  Output  Residue  Power  vs.  Normalized  Multipath 
Delay  for  a  Wideband  Jammer  and  am=.5,  aa=0 
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Figure  All.  Output  Residue  Power  vs.  Normalized  Multipath 
Delay  for  a  Wideband  Jammer  and  ■vl/T=.‘\ 


Figure  A12.  Output  Residue  Power  vs.  Normalized  Multipath 
Delay  for  a  Wideband  Jammer  and  Tj/T^.5 


Figure  A13.  Output  Residue  Power  vs.  Normalized  Multipath 
Delay  for  a  Wideband  Jammer  and  t1/T=1 


Figure  A14.  Output  Residue  Power  vs.  Normalized  Multipath 
Delay  for  a  Wideband  Jammer  and  AR«0 


Figure  A17.  Output  Residue  Power  vs.  Normalized  Multipath 
Delay  for  a  Wideband  Jammer  and  a  T/2  Delay 
in  the  Main  Channel 


Figure  A18.  Output  Residue  Power  vs.  Normalized  Multipath 
Delay  for  a  Narrowband  Jammer  and  am=1,  a ,=1 


Figure  A19.  Output  Residue  Power  vs.  Normalized  Multipath 
Delay  for  a  Narrowband  Jammer  and  am“0,  aa=0 
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APPENDIX  B 


INTERPOLATION  AND  EXTRAPOLATION  OF  A  BANDbIMITED  PROCESS 


B . 1  GENERAL 

Consider  a  random  process  x(t)  defined  on  the  infinite 
interval.  An  estimate  of  the  process  x(t)  can  be  obtained  as 
an  appropriately  weighted  sum  of  N  equally-spaced  samples  of 
the  process  covering  a  total  delay  T,  and  is  expressed  as 

?  ,nT 

x ( t )  =  l  an  x  (— J  .  (B-1) 

n  =  1  N 

An  illustration  of  this  process  is  shown  in  Figure  B1 . 


The  estimation  error  associated  with  (B-l)  is  then  given 


L**- 


as 


N 


e(t)  -  x ( t )  -  x ( t ) 


nT» 


(t)  “  1  anx( * 


n*1 


(B-2 ) 


Minimizing  the  mean-squared  estimation  error, 
e  »  E  {e2(t)},  is  equivalent  to  applying  the  orthogonality 
principle  [ B— 1 ] s 


N 


r  ,  mT.  "  (n-m)T 

E  [  e(t  )x(— )  ]  -  Rx(t-— )  '  l  anRx! - — ]  *  ° 

N  N  n=1 


N 


N 


m  =  1 , . .  .  rN 


(B-3 ) 


where 


R  (T )  =  E[ X(t)  X ( t+T ) ] 


(B-4 ) 


is  the  autocorrelation  function  of  the  process. 


The  optimal  weight  vector  aQ  is  obtained  by  solving 
(B-3)  which  in  vector  form  can  Be  written  as 


R  a 


( B-5 ) 


where  R  is  a  Nxn  covariance  matrix  and  r  is  a  N*1  column 
vector  with  elements  defined  in  (B-3).  Finally,  using  (B-2), 
(B-4),  and  (B-5),  the  minimum  mean-squared  estimation  error 
is  given  simply  as 


emin=  min  E{ e2 ( t )}  =  Px  -  =  Px-  a^  ,  (B-6) 

3 


where  Px  *  E{e2(t)}  is  the  power  of  the  process. 
B.2  BANDPASS  PROCESS 


Now  consider  the  special  case  of  a  band-limited  random 
process  with  a  frequency  spectrum  of  bandwidth  B  centered  at 
f0.  The  autocorrelation  function  of  this  process  is  the 
inverse  Fourier  transform  of  the  power  spectrum,  which  is 
then  given  as 


Rx(t) 


sine  Bt  cos  2it  f Qx 


t  .  ^o  x  . 

sine  —  BT  cos  [2it  —  BT  —  1  , 

T  V  B  T 


( B-7  ) 
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where  the  sine  function  is  defined  as 


.  A  sin  itx 
sine  x  =  — . —  -  , 

IT  X 

and  for  convenience,  the  process  power  has  been  normalized  to 
unity. 

Let  Rfl  denote  the  Nyquist  sampling  rate  which  is  equal 
to  twice  the  largest  significant  frequency  in  the  process 
power  spectrum  and  let  Rs  denote  the  actual  sampling  rate 
which  is  simply  given  by  N/T.  The  relative  sampling  rate  is 
then  defined  as 


A 

P  = 


'N 


N/T 

2fc+B 


N 


( -  +  1 )  BT 

B/f0 


(B-8) 


Calculations  of  em£n  as  a  function  of  p  and  normal¬ 
ized  point  of  interpolation  t/T  for  a  given  number  of  points 
N  and  fractional  bandwidth  B/f0  were  obtained  by  substi¬ 
tuting  ( B— 7 )  into  (B-6). 


The  solution  for  aQ  in  (B-5)  was  obtained  using  the 
Levinson-Robinson  algorithm  [32J.  Double  precision  arith¬ 
metic  was  employed  to  improve  accuracy.  A  small  positive 
term  AR  was  added  to  the  diagonal  of  the  covariance  matrix 
R,  whenever  deemed  necessary,  (for  large  N  or  small  B/f0), 

Tn  order  to  counteract  ill-conditioning.  The  ill-condition¬ 
ing  of  the  covariance  matrix  R  is  related  to  the  ratio  of 
largest  to  smallest  eigenvalue  or  equivalently  to  the  spec¬ 
tral  variation  of  the  sampled  power  spectrum  [33]  [20] . 


s,'(f) 


( B-9  ) 


For  a  bandpass  spectrum  there  is  no  overlap  in  (B-9)  for 
N>BT ,  and  hence  the  eigenvalue  spread  is  unbounded  and  the 
matrix  is  inherently  ill-conditioned. 


V 


Figure  B4.  Minimum  Mean-Squared  Estimation  Error 

vs.  Normalized  Estimation  Point  and  N=8 


Figure  B5.  Minimum  Mean-Squared  Estimation  Error 

vs.  Normalized  Estimation  Point  and  N=16 
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Figure  B8. 


Minimum  Mean-Squared  Estimation  Error 
vs.  Normalized  Estimation  Point  and 


B/f0=.5 


(.11 


Figure  B9.  Minimum  Mean-Squared  Estimation  Error 
vs.  Normalized  Estimation  Point  and 
B/ f 0= 1 
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APPENDIX  C 


SIMULATION  RESULTS  FOR  QUADRATURE  TAP  SPACING 


tn  this  appendix,  the  multipath  cancellation  performance 
of  the  adaptive  array  processor  considered  in  the  main  part 
of  this  report  is  investigated  for  quadrature  tap  spacing  as 
opposed  to  uniform  tap  spacing. 

The  steady  state  value  of  the  processor  adaptive  weight 
vector  Hb  is  obtained  as  a  solution  to 

R  ]£»  -  r  ,  (C-1  ) 

where  R  is  the  N*N  auxiliary  signal  covariance  matrix  and 
r  is  a  N* 1  column  vector.  For  a  uniform  tap 
spacing  iT/N,  i=0,...,N-1,  where  T  is  the  delay  line  length 
and  N  is  the  number  of  adaptive  taps,  the  matrix  R  is 
Toeplitz  and  can  be  solved  using  the  Levinson-RobTnson  algo¬ 
rithm  as  was  used  in  the  main  body  of  this  report.  For 
uniform  tap  spacing,  R  and  r  are  given  as 


R  J  (rij=  RJ 

T 

i)  -  ]}  i,j=0,...,N-1 

N 

(C-2 ) 

r  :  { r  *  Rma  (-j 

T 

-)}  j  »  0  ,  .  .  .  , N- 1 

N 

(C-3  ) 

where  Ra  is  the  auxiliary  signal  autocorrelation  and 
Rma  *s  the  cross-correlation  between  the  main  and  aux¬ 
iliary  signal. 

The  processor  output  power  Pr  can  be  easily  obtained 
as 


Pr  =  1  -  2$»TR  =  1  -  w.Tr  .  (C-4 ) 

Plots  of  Pr  (dB)  versus  x0/T  for  the  nominal  parameters, 
am=0,  aa=.5,  B=20  MHz  and  T=.1  ysec  were  generated. 

From  sampling  theory  it  is  known  that  in  general,  opti¬ 
mum  performance  will  be  obtained  provided  the  effective  tap 
spacing  is  closer  than  required  by  the  Nyquist  sampling 
rate,  that  is, 


Also  from  sampling  theory,  when  the  taps  are  located  on 
the  delay  line  in  quadrature  pairs,  optimum  performance  is 
obtained  provided  the  effective  tap  spacing  is  less  than  one 
over  the  bandwidth,  that  is 


A  plot  of  Pr  (dB)  versus  tq/t  for  N=4 ,8 , 1 6 , 32 ,64  and 
for  the  nominal  system  frequency  fo=100  MHz  is  shown  in 
Figure  Cl.  The  nominal  phase  of  the  auxiliary  tapped  delay 
line  signals  is  given  by 

iT 

=  2irf_  —  i  =  0,1,...,  N-1 .  (C-7) 

For  N=8 ,  it  can  be  seen  that  successive  taps  are  spaced 
apart  90°  in  phase,  and  thus  the  samples  occur  in  quadrature 
pairs  and  (C-6)  is  satisfied.  Therefore  near-optimum  per¬ 
formance  in  the  output  residue  response  is  obtained  for  N>8. 
On  the  other  hand,  for  N=4  the  taps  are  not  in  quadrature 
pairs.  Therefore  (C-6)  is  not  satisfied,  and  a  degradation 
in  the  Pr  response  is  noted. 

A  plot  of  Pr  versus  x0/T  with  fQ  increased  by  a 
factor  of  four  to  fc  =  400  MHz  is  shown  in  Figure  C2.  For 
N=8 ,  N= 1 6  a  performance  degradation  is  observed  which  is  the 
result  of  sampling  at  less  than  the  Nyquist  rate  (C-5).  For 
N=32 ,64  though  the  quadrature  spacing  of  taps  results  in  a 
near-optimum  performance  since  (C-6)  is  satisfied. 

A  plot  of  Pr  (dB)  versus  t 0/T  for  a  system  frequency 
five  times  larger  than  the  original  or  fo=500  MHz  is  shown 
in  Figure  C3.  For  N=8  it  can  be  seen  from  (C-7)  that  the 
taps  are  spaced  apart  90°  in  phase,  that  is  they  occur  in 
quadrature  pairs  and  (C-6)  applies.  Therefore  near-optimum 
performance  is  obtained  for  N>8.  Comparison  of  Figures  C2 
and  C3  shows  the  dramatic  performance  improvement  possible 
with  quadrature  sampling  in  Figure  C3  even  though  the  corre¬ 
sponding  system  frequency  is  larger  than  the  one  in  Figure 
C2  (500  versus  400  MHz).  The  previous  discussion  explains 
the  choice  of  number  of  taps  (N=8 , 1 6 , 32 , 64 )  in  the  simula¬ 
tion  runs  of  the  main  body  of  this  report. 

Another  way  to  ensure  quadrature  tap  spacing  is  by  sam¬ 
pling  the  auxiliary  signal  in  pairs  spaced  1/4f0  apart. 

In  this  case,  since  nonuniform  tap  spacing  is  employed,  the 
covariance  matrix  R  is  not  Toeplitz  and  a  solution  for  the 
steady-state  weight  w*  in  (C-1)  must  be  obtained  by  matrix 
inversion. 
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A  plot  of  Pr  (dB)  versus  tq/T  for  N=4,8  quadrature 
taps  with  single  precision  arithmetic  i$  shown  in  Figure  C4 
Near-optimum  performance  is  obtained  for  N=8  as  predicted 
since  (C-6)  is  satisfied.  The  degradation  for  N=4  can  be 
explained  perhaps  because  the  number  of  samples  N  is  not 
large  enough. 

A  plot  of  Pr  (dB)  versus  tq/T  for  N=4  and  N=8  for 
uniform  tap  spacing  was  generated  using  again  matrix  inver¬ 
sion  and  single  precision  arithmetic  and  is  shown  for  com¬ 
parison  in  Figure  C5.  Note  again  that  near-optimum  perfor¬ 
mance  for  N=8  is  obtained  due  to  quadrature  sampling  and 
that  performance  is  degraded  for  N=4. 

The  advantage  of  quadrature  sampling  as  opposed  to  uni¬ 
form  sampling  has  been  demonstrated.  For  a  real  signal 
processor  it  requires  an  inflexible  system  with  exact  align 
ment  of  the  taps  according  to  system  frequency  to  insure 
that  consecutive  tap  phase  differences  are  90*  apart.  A 
general  processor  on  the  other  hand  is  more  flexible  and 
will  have  good  performance  provided  (C-5)  is  satisfied,  tha 
is  N  is  large  enough. 
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Figure  Cl.  Output  Residue  Power  vs.  Normalized  Multipath 
Delay  for  a  Wideband  Processor  and  fo«100MHz 
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Figure  C3.  Output  Residue  Power  vs.  Normalized  Multipath 
Delay  for  a  Wideband  Processor  and  fo=500MHz 
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Figure  C4.  Output  Residue  Power  vs.  Normalized  Multipath 
Delay  for  a  Wideband  Processor  and  Quadrature 
Tap  Spacing 
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APPENDIX  D 


PLOTS  OF  JAMMER  AND  NOISE  AUTOCORRELATION  FUNCTIONS 


Consider  a  bandlimited  random  process  x(t)  with  a  fre¬ 
quency  spectrum  of  bandwidth  B  and  center  frequency  fQ. 

The  process  autocorrelation  function  Rx(t  )  is  given  as 

Rx(t  )  =  sine  Bt  cos  2irf0t  ,  (D-1  ) 

where  sine  x  is  defined  as 

sin  it  x 

sine  x  =  -  . 

H  X 

Plots  of  Rx(t )  versus  t  for  the  fractional  bandwidth, 

B/ f 0= . 01, .1,. 2, .5, 1,2  are  shown  in  Figures  D1-D6. 

Let  y(t)  represent  a  narrowband  random  process  con¬ 
sisting  of  two  sinusoids  of  frequencies  ft  and  f2  and  corre¬ 
sponding  powers  and  y2  •  The  process  autocorrelation 
function  Ry(x )  is  given  as 

Ry(T)  =  Y-j  COS  2irf1T  +  Y2  COS  2uf2x  .  (D— 2  ) 

A  plot  of  the  normalized  autocorrelation  function 
Ry(x )/{yl  +y2 )  for  the  nominal  parameter  values  fA =90  MHz, 
f ^  =  110  MHz,  y A  =  1 00  (20  dB)  and  y2  =  1000  (30  dB)  is  shown  in 
Figure  D7. 


APPENDIX  E 


COMPUTER  PROGRAM  LISTINGS 


E.1  WIDEBAND  JAMMER 

This  appendix  describes  the  computer  programs  that  were 
employed  in  the  main  body  of  this  report  for  solving  the 
adaptive  array  processor  interference  cancellation  problem 
in  a  single  multipath  environment.  For  the  wideband  jammer, 
output  residue  data  versus  normalized  multipath  delay  was 
generated  by  the  program-file  WIDBJAM  a  listing  of  which  is 
shown  in  Figure  El. 

The  program  employs  double-precision  arithmetic  for 
variable  definitions  and  calculations  in  order  to  obtain 
sufficient  accuracy.  It  starts  by  reading  in  from  file  24 
the  parameter  values  for  the  delay  length  T,  the  system 
frequency  F  and  bandwidth  BW,  the  jammer  signal  to  noise 
power  ratio  G  and  the  multipath  component  strengths  AM  and 
AA.  Since  the  auxiliary  channel  covariance  matrix  is 
Toeplitz,  a  one-dimensional  array  R(N )  is  used  to  store  its 
elements.  The  cross  correlations  of  the  main  with  the  de¬ 
layed  auxiliary  signals  are  stored  in  array  B(N).  Loop  2 
calculates  the  elements  in  R(J),  B(J),  J=1,...,N  by  calling 
functions  RA  and  RMA .  Functions  RA,  RM,  RN  and  RMA  calcu¬ 
late  for  a  given  delay  TAU  the  corresponding  auxiliary,  main 
and  noise  signal  autocorrelations  and  the  crosscorrelation 
between  the  main  and  auxiliary  signals. 

A  small  positive  number  RD=10-3  is  added  to  R( 1 )  to 
counteract  the  ill-conditioning  of  the  auxiliary  covariance 
matrix.  A  solution  for  the  adaptive  weights  W(N)  is  ob¬ 
tained  from  subroutine  SOLVE,  which  is  a  modified  version  of 
a  subroutine  listed  in  Appendix  A  of  [34],  based  on 
Levinson's  method. 

Subroutine  BERROR  calculates  the  average  and  maximum 
errors  in  the  estimates  of  B(N)  due  to  the  finite  accuracy 
in  the  calculation  of  W(J),  J=1,...N  and  the  ill-condition¬ 
ing  of  the  covariance  matrix. 

The  output  residue  PT  is  calculated  in  decibels  by  loop 
4  and  the  five  ensuing  statements.  If  the  output  residue  is 
less  than  or  equal  to  zero,  due  to  finite  computer  accuracy 
or  ill-conditioning,  its  decibel  value  is  set  by  statement  6 
to  -30. 

Finally  the  program  writes  the  output  residue  values  PT 
and  the  corresponding  normalized  multipath  delays 
TR1  =  TO/T  =  0,  .01  ,  .02, . . . , 1  .80  to  file  101. 


E.2  NARROWBAND  JAMMER 

The  computer  program-file  N ARB JAM  generates  output  resi¬ 
due  versus  normalized  multipath  delay  data  for  the  narrow- 
band  jammer  adaptive  cancellation  problem.  The  program  is 
a  variation  of  the  wideband  jammer  program  WIDBJAM.  It 
reads  from  file  3  the  delay  length  T,  the  frequencies  of  the 
two  sinusoids  F(1)  and  F(2),  the  system  bandwidth  BW,  the 
sinusoid  signal  to  noise  power  ratios  G ( 1 )  and  G(2)  and  the 
multipath  parameters  AM  and  AA.  The  program  has  the  addi¬ 
tional  function  RS  which  calculates  the  autocorrelation  of 
the  jammer.  Its  main  difference  to  WIDBJAM  is  that  it  uses 
complex  double  precision  arithmetic  to  calculate  the  output 
residues  due  to  the  sinusoids  in  loop  7.  Finally  the  pro¬ 
gram  writes  to  file  102  the  normalized  delay  values 
TR1=0, .01 , .02, . . . , 1 .8  the  total  output  residue  PT  and  the 
individual  power  residues  PR(K ) ,K=1 , . . . ,3. 

E .3  INTERPOLATION  AND  EXTRAPOLATION  OF  A  BANDLIMITED  PROCESS 

The  program-file  INTERPOL  in  Figure  E3  generates  output 
data  for  the  minimum  mean-squared  estimation  error  (in  deci¬ 
bels)  EMIN  of  a  bandlimited  process  as  a  function  of  normal¬ 
ized  estimation  point  TT.  The  program  is  very  similar  to 
WIDBJAM  and  NARBJAM.  For  a  given  number  of  taps  N,  func¬ 
tional  bandwidth  BF  and  relative  sampling  to  Nyquist  rate 
RT(J),  1=1,... ,13,  it  generates  output  data  of 
TT=0 , .0025, . . . ,4  and  EMIN.  Its  function  RX  calculates  the 
autocorrelation  of  the  bandlimited  process  and  is  essential¬ 
ly  the  same  as  function  RN  of  programs  WIDBJAM  and  NARBJAM. 

E.4  PLOT  ROUTINES 

A  plot  of  the  output  residue  PT  versus  normalized  multi- 
path  delay  TR1=0, .01 , . . . , 1 .8  was  generated  by  the  program- 
file  P LTD AT .  The  program  reads  the  number  of  taps  N  and  the 
corresponding  TR1 ,  PT  values  from  file  IF  and  stores  them  in 
array,  X(180),  Y(180).  For  each  N  it  calls  subroutine  LINE 
to  plot  the  points  X(J),  Y(J),  J=1,...,180.  The  X  and  Y 
axes  are  plotted  by  calling  subroutine  AXIS  twice. 

The  program  PLTDAT  was  also  used  to  generate  plots  for 
the  narrowband  jammer  after  simply  changing  the  format  of 
the  READ  and  WRITE  statements  in  it.  For  the  interpolation 
plots  in  addition  to  changing  the  READ  and  WRITE  statements, 
the  dimension  of  X  and  Y  was  increased  to  1600  and  IN 
changed  to  1600  since  a  total  of  1600  points  were  read  and 
plotted.  A  PLXY  plotter  was  used  to  generate  all  the 
plots . 
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FILE*  WIDBJAM 


1,  THE  JAMMER  CONSISTS  OF  A  WIDEBAND  SIGNAL  (30  DB) 

2.  ALL  VARIABLES  AND  CONSTANTS  ARE  IN  DOUBLE  PRECISION 

3.  A  TYPICAL  ANGLE  IS  EVALUATED  ASt 
ONEGA  a  MOD(2.DO*F(I)*TAU,2.DO)*PI 

4,  AA  AND  AM  ARE  THE  MULTIPATH  GAINS  FOR  THE  AUXILIARY  AND  MAIN 
CHANNELS  correspondingly, 

DOUBLE  PRECISION  R ( 260  )  ,  B(  2o0  )  ,  w  (  2b0 )  ,  Y  1 ,  Y 2 ,  YT 

DOUBLE  PRECISION  PT, XN , PI , T , TO , T1 , TR , TR1 , XX , Y Y 

DOUBLE  PRECISION  BHAT( 260) , BW , F ,G, Y ,P, KMA, HA , RM, AA , AM ,RD 

COMMON  AA, AM,T0,T1 ,G 

COMMON  /B1/F,BW,PI 

PI*3.14l592b535B97932Df0 

RD*l.l>*9 

IF* 101 

READ  (24,*)  T,F ,B«,G,AM,AA 

THE  FOLLOWING  NORM ALIZA T1UNS  ARE  ASSUMED  IN  THE  DEFINITIONS 
OF  THE  PROGRAM  VARIABLES* 

KaFM.E-6  3wsBW*l,E-6 

T0*TU* I , E*6  Tl*Tl*l,EFb  f*T*l,E4b 

Tl*l,UD-4 

TRsTl/T 

WRITE  (IF, 100)  F,Bw,G, AM, AA,T,TR,DD 
100  FORMAT  (/,5X,'F0  * ' , F7 , 2 , 23X, *  NOISE  bAnDwIDTH  *  * ,FS , 1 , 10A, » (MEG 
»  AHEKTZ ) * ,/ ,SX, 'GAMMA  s', F7 , 1 , 20X ,' ALPHA  OF  MAIN  *  ',F8,6,iOX, 

, 'ALPHA  OF  AUXILIARY  *  ♦ , FB, b , / , SA , ' TAP 

, DELAY  LENGTH  ** ,FS.3,tlX, 'Tl/T  s', r 7 , 4 , l 9X ,»( MICROSECONDS) ',/ , 
,SX, 'WHITE  NOISE  ADDED  TO  HA(I)  *',Dl0,2) 

N*B 

DO  1  MSI, 4 
YTSO.DO 
TO*-.OlDO*T 
WRITE  (IF, 200)  N 

200  FORMAT  (//,5X,»#  OF  DELAY  TAPS  s' ,13) 

XNSDBLE(N) 

DO  J  L*l,lBO 

TO*TO+,OIDO*T 

TRlsTO/1 


Figure  El.  Wideband  Jam?.ier  (WIDBJAM)  Program  Listing 
(Sheet  1  of  4  Sheets) 
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DO  2  J=i , n 
Yl*DBLE(J-l)*l/XN 
Y2sUttLE(-J-fl  )*I/Xn 
R(J)sRA(Y1) 

2  B(J)=RMA( X2) 

R(1)sR(1)*RD 

CALL  SOLVE  (R,w,b,N) 

R(t)sR(l)-RD 

CALL  ERROR  (R  ,  *  ,  b ,  BHA T ,  XX ,  Y  Y ,  N  ) 

YTsYl+YY 
PaO.DO 
YaO.DO 
DO  4  K»1,W 
DO  b  J=1#N 
JJ=iABS(K-0)+l 
b  PsP+w(M*«lO)*R(JJ) 

4  Y*Y4W(K)*B(K) 

PT=HM(0.D0)tP-2.D0*Y 
IF  (RT.LE.O.DO)  GO  TO  b 
RT®1 0 , DU 4  0 LOG  10 ( PT) 

GO  TO  3 
6  RTs-30, DO 

3  wWITE  (IF, 300)  TH1,RT,XX,YY 

300  FORMAT  (lA,'TO/T  a  « , F4 , 2 , 1  OX , * RR  a  •  ,  F 1 o , 1 U , 1  OX , 

,'RAX,  NORM,  6  ERROR  a »  ,  01 4 , 7 , 5X , ’ A VE ,  NORM.  b  ERROR  a* ,014.7) 
Y  T  a | 1/180, DU 
WRITE  (IF,4uO)  IT 

40U  FORMAT  (lUX, 'AVERAGE  NORMALIZED  b  ERROR  a  «,U14.7) 
l  N®N42 
STOP 
KtiO 

FUNCTION  RA  (TAJ) 

COMMON  AA,Am, ro,Tl ,G 

DOUBLE  PRECISION  A  A  ,  AM  ,  f  AU  ,  G  ,  l’U  ,  !  I  ,  S  1  ,  S2  ,  S4  ,  S4  ,  R  A  ,  Rn 

SIsrN  U’AU,G) 

S2*RN  ( TAU-TO-Tl ,G) 

S3*RN  ( TAU+I0+T1 ,G) 

S4sRN  ( T  A  0 , 1,JUU) 

RA*(  1  ,D04AA**2)*S1*AAMS2  +  Si)+S4 

RETURN 

END 


Figure  El.  Wideband  Jammer  ( WIDBJAM)  Program  Listing 
(Sheet  2  of  4  Sheets) 


FUNCTION  KM  (TAU) 

COMMON  AA,AN,TQ,T1,G 

DOUBLE  PRECISION  A A , AM , TAU , G , TO , T 1 , S 1 , S2 , S J , RM , RN 
SI *HN  (TAU,G) 

S2*KN  (TAU-T0,G) 

S3«RN  (TAU*T0,G) 

RMa(l,DO+AM**2)*Sl+AM*(S2+Sj) 

RETURN 

END 


FUNCTION  KMA(TAU) 

DOUBLE  PRECISION  A A , A N , G , TAU , TO , T 1 , S 1 , S2 , S3 , S4 , RMA ,  RN 
COMMON  AA , AM  ,  TO , T 1 , G 
S1*RN  (TAU, G) 

S2*RN  (TAU-TO-TI,G) 

53*RN  ( TAU+TO , G ) 

S4*RN  (TAU-T1,G) 

RMAsSlTAA*S2+AM»S3+AA*AM*S4 

RETURN 

END 


FUNCTION  RN  (TAUrGl ) 

DOUBLE  PRECISION  F, G1 ,B*,TAU,X1 ,X2,XJ,RN,PI 

COMMON  /Bl/F,BM,PI 

RN*G1 

IF  (TAU.EU.O.ODO)  RETURN 
XlaBN*TAU 
X2«2.0D0*F*TAU 
X3*DMOD(Xl ,2«OUO)*PI 
XlsXl*PI 

X2aUMOD(X2,2.000)*PI 

RNaGl*DSIN(A3J*DCOS(X2)/Xl 

RETURN 

END 


SUBROUTINE  SOLVE  (R,A,B,M) 

DOUBLE  PRECISION  R ( 260 ) , A ( 260 ) , B ( 2bQ) , F ( 260) ,TF( 260 ) 
DOUBLE  PRECISION  ALPHA , BETA , GAMMA ,C ,U 
F(1 Jal.DO 
ALPHAaR(l) 

BETAaR ( 2 ) 


Figure  El.  Wideband  Jammer  ( WIDBJAM)  Program  Listing 
(Sheet  3  of  4  Sheets) 
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Acn=aci)/HU) 

GAMMA*A(1 )*K(2) 

00  1  N =2,M 
NMlSN-1 
C=-BETA/ALPHA 
IF  (N-2)  2,2,3 

3  DO  4  Js2,MU 
NNsn-J+1 

4  TE(J)=F(J)+C*F(NN) 

DO  t>  J32,MM1 

5>  F(J)*TF(J) 

2  F  (  N  ) *C*F  ( 1  ) 

ALPHAsAlPHA*C*B£TA 

BETAsO.DO 

DO  6  J=i,* 

NNsN-J+2 

t>  BETA*BETA-*F(J)*R(NN) 

Os(B(n)-GA*mA) /ALPHA 
DO  7  J  * 1 , N  M  1 
NfVSN  — J+l 

7  A(JJs:A(J)-KJ*F(mn) 

A(N)sU*F( 1  ) 

GAmmAsO.DO 

DO  8.J*1,N 
NNsn-J+2 

8  GAM.<!A=GAM*A  +  A(JJ*R(NN) 

1  CONTINUE 
RETURN 
END 

SUBROUTINE  ERROR  (  R  ,  # , B , BH AT , X , Y , N) 

DOUBLE  PRECISION  R  ( N  )  ,  V*  ( N )  ,  6( N )  , BHA T C N  J  , SUM , X ,  ¥ 
X=1 ,0-16 
YsU.DO 
DO  1  J=l,» 

SURsO.DO 
DO  2  K*1,n 
KKsl+J-K 

IF  CK.GT.J)  KK*1+k-j 

2  SUMsSUM4R(KK)4m(K) 

IF  (OABSCbCJ) J,LT, 1,0-127  GO  TO  3 
bHAT(J)*DAbS( (SUN-B(J) J/B( J) ) 

3  Y® Y+BHAT (  J  ) 

1  IF  C  BHAT ( J ) ,GT , X )  XsBHAT(J) 

YsY/DBLECN) 

RETURN 

END 


Figure  El.  Wideband  Jammer  (WIDBJAM)  Program  Listing 
(Sheet  4  of  4  Sheets) 
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FILE!  NARttJAM 


ALL  VARIABLES  AND  CONSTANTS  ArE  In  DOUBLE  PRECISION 
A  TYPICAL  ANGLE  IS  EVALUATED  AS: 

OMEGA  a  MOD(2.D0*F(I)*TAU*Pl  ,2.D0*Pi) 

AA  AND  Am  ARE  THE  MULTIPATH  GAINS  FOR  THE  AUXILIARY  AND 
MAIN  CHANNELS  CORRESPONDINGLY 


COMPLEX* 1 b  SUM ,  C  ,  X  (  J  ) 

DOUbLE  PRECISION  H( bOO ) , U(60Q ) , W ( bOU ) , PR l i J , PT , Y 1 , Y 2 , YT 

DOUHLE  PRECISION  RK ( bUU J , Y ( 3 ) , XX , XK , TAU , S4 , P I , T , TO , T 1 , 1R , TH 1 ,Rd 

DOUbLE  PRECISION  BHAT ( bOO ) , Bw , F ( 3 ) , G l 2 ) , I Y , XT , P , K A , RM A , RN , A A , AM 

COMMON  A A,  AM,  TO, T1 

COMMON  /BI/K,P1 

CUMMOw  /B2/G 

COMMON  /Bi/rt* 

COMMON  /B4/S4 

PI  =  J.  I4l592bSiSHS  /932D+0 

RDsl .OD-9 

IFX102 

READ  (3.*)  l,(F(l),lal,2),bw,(G(l),ial,2),AM,AA 


m 


EMM 


THE  FOLLOWING  NORM AL IZA TiUrtS  ARE  ASSUMED  iw  THE  DEFINITIONS 
OF  THE  PROGRAM  VARIABLES! 

FU)=Fl*l.t-b  E(2)sF2*I  .E-b  B*=b**l.E-o 

TOa'lO*  1 .  E+b  TlaTl*l.K+h  T=T*l.E+b 


^  •  .*• . 


Tlsl.U-4 

TRsTl/T 

F(3)a(F(l )*F(2))/2.D0 

WRITE  (IF,loU)  (F(I),Isl,2),Bw,(GU>,isl,2),AM,AA,T,TR,Ol> 
100  FORMAT  ( / ,bA ,  •  FI  a  •  , FS.  1 , 22X ,  •  F2  a  »  ,Fb.  1  ,  lbX,  •  NOISE 
, BANDWIDTH  s  •  ,FS.l,iOX, • (MEGAHERTZ) ’ , 

,/,5A, 'GAMmAI  a  •  , Kb . I , 1 HX , • GAMM A2  a  • 

,Fb.l,  10X,  'ALPHA  OF  MAIN  a  ',FH,b,10A, 

, 'ALPHA  OF  AUXILIARY  *  *,F«,b,/,5X, 

, 'TAP  DELAY  LENGTH  a  ' , Fb . 3 , 1 0 A , »  T 1 / T  =  ',D12.4,3bX, 

,' (MICROSECONDS) ' ,/,5X, • WHITE  NOISE  ADDED  Tu  RA(1)  a',ui0.2) 
Nat) 

DO  I  Mai, 4 

YTau ,uO 

WRITE  (IF, 200)  N 

200  FORMAT  (//,bA,'«  OF  DELAY  TAPS  a  ',13,/) 

XivSi/hLE  (  rt  ) 

TO*-,  01  DO*’) 

WRITE  (IF, 3 00) 

300  FORMAT  (lA,  'TO/I' , 1 2A , • PH • , I 4X  , 

, 'PHI' ,14X,'PH2' ,14X, 'PRO* ,12A, 'MAX.  NORM.  H  ERROR ', ttX ,» AVER , 
,  norm,  w  ERROR',/) 


•  V-  V 
•-■'.‘.NV 


:■  \\.\j 

--  *W  N  J 


•v.  -i 
i-  .*#1 


Figure  E2.  Narrowband  Januner  (NARBJAM)  Program  Listing 
(Sheet  1  of  4  Sheets) 
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DU  i  Lsl, 1  HU 

'l  OsTO-f  ,  0 1  00*'f 

•IklsTO/f 

DU  2  J=1,N 

YlaDBLE(l-J)*i/XN 

Y2sOBLK(-J+1)*T/Xn 

H  (  J } SKA (  Y  1  ) 

kk( J)sS4 

2  b(J)sRtfA(Y2) 

k  (  1  )  SR  (  1  )  +RL) 

Call  solvk  ik,«,b,nj 

HUJSKU  J.Ru 

CALL  fclHnUR  (k,W,H,BHAT,XX, |Y,N) 

Y'fsYT  +  Y  Y 
DU  7  Ks  1  #  ^ 

SUrtsDCMkLA l 0,00, O.DU) 

DU  4  .isl#,N 
X  IS0BL£(  J-l )*f/XW 
Y(U*2.D0*F  (K)*XT 
If  (2)*2,U0*F  (K)*(Xl>l’0  +  Tl) 

Y ( 3 ) *2 ,  DO*!*  (K)*ro 

DU  H  D Js t ,3 

Y(JJ)  snmUD ( i ( JJ ) , 2 . UDO )  *Pi 
«  X(JJ)sDC*PLX(0.l)0,-Y(Ju)  ) 

4  SlJMsSUtf  +  *cU)*(COfcXP(X(  1  )  )+AA*Cl»KXPUU)  )  ) 
C=l,D0+AM4CDKXPCXl JJ J-SUh 
7  Pk(K)sG(K)*(CDABS(C)**2) 

P=0,DU 

DU  b  K  s  1  ,  in, 

DU  b 

JJsiAMSl A-J J+i 
b  psp+»(h )**lJ)*kH(JJ ) 

PH ( 3 ) sp 

PI  sPR ( 1 )+Pk 12)+PH(J) 

PTs10.du*ULUG10(PX) 

DU  h  Ks  1  ,  J 

b  PHir)slO.UDADLUGlUCPK(K)  ) 

*Kllfc  (1F,4U0)  fkl ,Pt, (PH(K) ,Ksi r J),xX, it 

400  FOKriAT  (lX,F4.2,4(bX,FU,o),2UOX,Oi4.7)) 
i  CUtfll>»Ut 
Y  1  =  Y1/1b0,I>0 
wHIIF  t  l  F ,  b  0  0  )  YT 

boo  F  0  k  M  A  '£  (  lux  ,  *  AVERAGE  uOHMAUZED  k  trtKUK  s  ',014,7) 

1  usMY 
STUP 
tOD 

KUNCllUiV  KA(TAU) 

CUWMJN  A  A  ,  A  >1  ,  f  0 ,1 1 
CUh.-sUN  /  B  4  /  v>  4 

DUOoLfc  PRECISION  AX,A«,  rAU,tO,l|  ,51,S2,SJ,S4,HA,HS,Hi» 
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ALPmAsK {  1  ) 
bt TA  =  H(2) 
A(l)=tt(l)/klU 
GAMMAsAU  )»K(2) 

00  1  n  =  2  ,  M 

Ah"1  1  =N«  l 
C=-Ht:TA/ALPnA 
IK  ('•-2)  2,2,3 

3  00  4  .1=2,  Mr  l 
.MN  =  h-J  +  l 

4  TF(0)=KJ)  +  C*b (hn) 

0U  b  0  =  2,  NM 

b  F(0)  =  J+(01 
2  F(M=C4FU> 

AOPHAsAUPriA+ObK  TA 
BfcJ  A  =  0,00 
0U  b  0=1, « 

(vfj  =  fv».)  +  2 

b  bb.  I  4= OFT A +b  l  J  )  *d (J H  J 
0=  (  b  C  **  )  »GAb  iA  1  /  ALPb  A 
00  /  0=1,  NM 

I.)  IJSII  m  0+1 

/  A  (  J )  =A  ( .J  )  fo  +  F  (I'M  ) 

A(Nl=0*f  C  1  ) 

GAi*b4  =  0,L>0 

uu  n  0  =  1, n 

.D  ID  =  iD  -  0  +  2 

«  GAMrl  a=GAMmA  +A  CJ  )  +  K  ( ‘Mu  ) 
1  COm’1  l  wUE 
HFTOM 
fc  U  0 


•V-V-V 


r.| 


■’WWW 


SOBKOUilNt  FHKUK  (K,  «*  ,  h,  UhA  t  ,  A 

OUliHLt  PKtCl.blU  D  K  (  Olio  )  ,  w  (  OU0  ) 

A  =  1 ,00- lb 

Y  =  0  «  D  0 

OU  1  0=1  ,  *D 
SUM=U,D0 
00  2  K= 1 , d 
KKsi+0-k 

lb  (a.GT.O)  Kk=l  +  t\~0 
SUM=SUN  +  R(l’A)*rt(M 
lb  (i)AB.S(bl  J  )  )  .  Ij  r  .  ID- l  2  >  GO  lu 
riMA  l'CO  )=0A»  .>(  )  )/ti(0  )  ) 

Y  =  Y  +  bM A  I C  o  > 

lb  (<3rtA  1  (0  )  ,G  t  ,  A  )  A  =  r)ilAl  (0> 

Y  =  Y  / 1)  M  L  e,  C  '•» ) 

K  K  1  0  H  id 

b.fii> 


(  K  ,  *>  ,  H  ,  i3  M  A  1  ,  A,  Y  ,  N  ) 

KloUO) , w(ooo) ,hloUO) , oHaT (  bOO ) ,burt,X, Y 


* »  \- 
dV-W.n 


v 


--  > 


::-o% 

*  /  .* 
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‘V-Kn 

‘»V<. 


V.» 

•v-v-> 


Si;? 

:i-: 


J  > 


c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


fill:  InTEKROL 


1. 


2. 

3. 


4. 

5. 


THIS  PROGRAM  INVESTIGATES  THE  INTERPOLATION  OF  A  BANDLIMITED 
PROCESS  OF  DURATION  T  AS  A  FUNCTION  MAINLY  OF  THE  LOCATIUN  OF 
THE  POINT  TAU  WHERE  THE  PROCESS  IS  INTERPOLATED. 

THE  PROCESS  IS  SAMPLED  AT  N  POINTS. 

THE  MINIMUM  MEAN  SQUARED  ERROR  (EMIN)  IS  CALCULATED  IN  DB 
AS  A  FUNCTION  OF  THE  BANDWIDTH  TIME  PRODUCT  BT  (B*T),  FRACTIONAL 
BANDWIDTH  BF  (B/FO),  THE  NUMBER  OF  SAMPLES  N,  AND  THE  NORMALIZED 
POINT  OF  INTERPOLATION  TT  (TAU/T), 

ALL  VARIABLES  AND  CONSTANTS  ARE  IN  DOUBLE  PRECISION 
A  TYPICAL  ANGLE  IS  EVALUATED  ASt 
OMEGA  s  MOD(2,O0*F(I)*TAU,2.D0)*PI 


DOUBLE  PRECISION  R ( 260 ) , B{ 260 ) , w ( 260 ) , BNAT ( 260) 

DOUBLE  PRECISION  EM1N.XN.PI.TT.XX.YY.Y1.Y2.YT 

DOUBLE  PRECISION  BF , BT, Y , P , RX ,|*D ,DTT , RT ( 1 3 ) 

COMMON  /Bl/BF.BT.PI 

PI*3, 1415926535897932 D+0 

RDbO.DO 

BF*.  2 

RT( 1 )*1 . 

RT(2)»1.1 
RT(3)«1,2 
RT(4)«1 .5 
RT(5)*2. 

DO  10  Js6,13 
10  RT(J)sRT(J-l)*2. 

IF*101 

N*2 

XNBDBLE(N) 

DTT«.0025D0 

WRITE  (IF. 100)  BF.N.RD 

100  FORMAT  (/.5X. 'FRACTIONAL  BANDWIDTH  ■  '.FJ.1.10X. 
SAMPLES  s' .13. 

,5X, 'WHITE  NOISE  ADDED  TO  RA(1)  *',D10.2) 

DO  1  Mai. 13 

YTsO.DO 

TT*»DTT 

BT*XN/(RT(H)4 ( 1 , DO+2 . DO/BF ) ) 

WRITE  (IF. 200)  RT(M).BT 

200  FORMAT  (/.5X. 'RATIO  OF  SAMPLING  TO  NYUU1ST  RATE  < 
.F8.1.10X. 'TIME  BANDWIDTH  PRODUCT  s  «.F9.6) 

DO  3  L*l. 1600 

TTsTT+DTT 

DO  2  Jsl.N 

Y1*DBLE(J-1)/XN 

Y2*TT*DBLE(-J)/XN 

R(J)sRX(Yl) 

2  B( J)sRX( Y2) 


Figure  E3 


Interpolation  and  Extrapolation  of  a  Bandlimited 
Process  (INTERPOL)  Program  Listing 
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s' ■-v-v-yv'v'v'V 


.  ■*.  - 


'-'••I 


s  •- 


*■  IMl 


•.Os 


■'.VS 

•.v.sy. 
■*  yCy 


fswsewl 


•V\V\. 


'<  ■  .*•! 


V.s.s-’ 


V.-'.'A 


*  I1* 

*  «  *  -  -  • 
v  J 
vV/j 


R(l)sR(l)+RD 

CALL  SOLVE  (R,W,B,N) 

RCl)*R(l)-RO 

CALL  ERROR  (R , M , 8, BHAT , XX , YY , N ) 

YTaYT+YY 

PaQ.DO 

Y»O.DO 

DO  4  Ks1,n 

00  5  J«1,N 

JJ*IABS(K-J)+1 

5  P«P+W (K)*M(J)4R(JJ) 

4  Y«Y+W(K)4B(K) 

EMIN*1.0D0+P-2.D0*Y 
IF  (EMIN.LE.O.DO)  GO  TO  6 
EMI Nx} 0 , DO*DLOG 1 0 ( EM  IN ) 

GO  TO  3 

b  EMIN*-100.DO 

3  WRITE  (IF » 300)  TT,EMIN,XX,YY 
300  FORMAT  (1X,'TAU/T  a  ',F7,5,5X, 

,  'MIN,  MEAN  SO.  ERROR  a  '  ,  F 1 1 , 6 , 5X , ' M AX  NORM.  B  ERR •  a  ' 

, ,D14,7 ,5X, ' AVE.  NORM,  B  ERR.  a  »,Di4.7) 

YTaYT/lbOO ,00 
WRITE  ( IF ,400)  YT 

400  FORMAT  (1QX,'AVE,  NORM.  B  ERROR  a  * ,014.7) 

1  CONTINUE 
STOP 
END 

FUNCTION  RX  (TAU) 

DOUBLE  PRECISION  BF , BT , TAU , X 1 , X2 ,X3 , KX , PI 

COMMON  /B1 /BF , BT #PI 

RXal.ODO 

IF  (TAU.EO.O.ODO)  RETURN 
X1«BT*TAU 

X3*DMOD(X1,2,ODO)*PI 

XlaX14PI 

X2*2 • ODO*BT*TAU/BF 

X2»DMOO(X2 ,2.000)*PI 

RXaOSIN(X3)*DCOS(X2)/Xl 

RETURN 

END 

SUBROUTINE  SOLVE  (R,A,B,M) 

DOUBLE  PRECISION  RC260) , A(260) ,b(2bO) ,F(260) ,TF(26Q) 

DOUBLE  PRECISION  ALPHA , BETA , GAMMA rC , U 
F(1)*1,D0 

Figure  E3.  Interpolation  and  Extrapolation  of  a  Bandlimited 
Process  (INTERPOL)  Program  Listing 
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ALPHAsRd) 

BETA«R(2) 

A ( 1 ) *B ( 1  )  /R  (  1 ) 

GAMMAsACl )*R(2) 

DO  1  N*2 , M 

Cb-BETA/ALPHA 
IF  (N-2)  2,2,3 

3  00  4  Jb2#NM1 
NNaN-d+1 

4  TF(J)*F<J)+C*F(NN) 

DO  5  d=2»NMl 

5  F(d)»TF(d) 

2  F(  N  ) »C*F ( 1 ) 

ALPHASALPHA+C4BETA 
BETA*0 , DO 
00  6  J*1 , N 
NNsN-d+2 

b  BETA«BETA^F( J)4R(NN) 

0=(B(N) -GAMMA) /ALPHA 
DO  7  J«1,NM1 
NN=N-J+1 

7  A(J)*A(d)+0*F(NN) 

A(N)»0*FC1) 

GAMMAbO , DO 

DO  8  J*1 , N 
NKBN-J+2 

8  GAMM AbGAMM A+A (J)*R(NN) 

1  CONTINUE 
RETURN 
END 

SUBROUTINE  ERROR  (R , « , B, BHAT , X , Y , N ) 

DOUBLE  PRECISION  R ( N ) , W ( N) , B( N ) , BHAT( N ) , SUM , X , Y 

Xal.D-1 6 

YbO.DO 

DO  1  J*1,N 

SUmbO  a  DO 

DO  2  Kb1,n 

KKbI ♦ J-K 

IF  (K,GT,d)  KKbI+k-J 

2  SURsSUM+RCKK)*N(K) 

IF  (DA8S(B(d)),LT.i.D-12)  GO  TO  3 
BHAT( J)bDABS( (SUM-B(d) )/B(J) ) 
i  YaY+BHAT(d) 

1  IF  ( BHAT (J ) , GT ,X )  XaBHAT(d) 

YbY/DBLE(N) 

RETURN 

END 


Figure  E3.  Interpolation  and  Extrapolation  of  a  Bandlimited 
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c 

c 

c 


FILES  PLToA'l 


dIMeasJOh  x(  i«2),Y(182) 

DUUul.K  PRECISION  ill) 

1 «  *  1  8  0 

lb  slot 

IGsiO 

Hb.Au  (lb, 100)  K,*M,G,  AM,AA,  1’,  IK, HO 
100  FORMAT  ( /  ,  9  A  ,  F 7 . 2 , 40X ,  F5  ,  1 ,  / ,  1 2X  ,  F7 . 1 ,  JtoA  ,  Kb  .  6 , 3 1 X  , 
,Fh.o,/,23X,F5. 3,  WX,F/,4,/  ,  1 3X  ,G10,2) 
wRiiE  ( 1 G ,  2  0  0  7  F ,B*,G,A3,AA,  T'THfRV 
200  bliKwAJ  (/,5a, 'FO  a'  ,FJ.2,23X,  'NOISE  dAi«D*IUTH  a  •  ,  K*> ,  1 , 10A ,  '  C^EG 
,AHEKI7.)  •  ,/,5A,  'GAMMA  s«,K  7  .  t ,  2<>X  ,  '  ALPH  A  uF  F.A  IN  a  '  ,  Ffc  .  fa ,  1  OX , 
,'AlibnA  UF  AUXILIARY  a  •  ,  F8  ,  t>  ,  /  ,  5X  ,  •  TAP 

, DELAY  LFNbl'ri  a',F5.3,UX,'Tl/T  * ' , F 7 . 4 , 1 9A , » ( *  ICR  OSECOR05 ) • , / , 

,  5X  ,  '  aHITE  \U1.SE  ADDED  TO  HA(1)  a*, 010. 2) 

00  1  J*l,4 
HEAD  (IF, 400)  N 
300  FORMAT  (//,22a, 13) 
isun:  (IG,400)  n 

400  FORMAT  (//,bX,'»  OF  DELAY  IaPS  s',  13) 

DO  2  1*1 f (N 

READ  (IF, 500)  X(I),Y(I) 

IF  l  Yd)  .bT.-JO.O)  ¥(l)a-30.0 

2  CONTINUE 

500  FORMAT  («X,F4.2,15X,F16.10) 

IF  (J.GT.l)  GO  I’O  3 
CAUL.  PU.TST  (1.,2RIN,1) 
i  (  In  +  1  )3- JO . 

Y  ( li**2)ai0. 

X(  liHDaU, 

X ( Ir+2)s,2 

CALL  AXIS  (.5,.5,4riT0/T,-4,9,,0.,XUN41),X(lN  +  2)) 

CAUL,  AXIS  ( .5,  ,5,oHPH(DB)  ,6,6,  ,90,  ,Y  l  INF  1  )  ,¥  t  IN  +  2 )  ) 

CaLL  PLOT  (,5,.5,-3) 

3  CALL  LIrE  (X,  i ,  Irt,  1 ,0, 1 ) 

1  READ  (IF,*) 

CALL.  PLOT  (0,,0,,-3) 

CALL  PLOTrtD 
CALL  EXIT 
STOP 
END 


Figure  E4.  Plotting  Routine  Program  Listing  (PLTDAT) 
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